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Preface 



Many books on mesoscopic systems have been published as progress has con- 
tinued in the fields of nanoscience and nanotechnology. The focus in these 
books is mainly on quantum mechanical behavior in artificial electronic sys- 
tems fabricated by nanometer-scale structuring. Such quantum mechanical 
behavior is projected to macroscopic observers and the quantum nature can 
be utilized in practical devices. Quantum computers, another hot topic nowa- 
days, are characterized by excitation coherence properties among nanostruc- 
tures, and the ability to maintain excitations is very important when using 
the characteristics as information. In that sense, the device is described as 
a microscopic system and some processes occur before being projected to 
macroscopic observers. In this book, the authors try to describe not only 
the techniques for fabricating nanostructures but also new directions as re- 
gards exciting systems and understanding how energy is dissipated through 
observation. 

The idea of ‘nano-electrodynamics’ underlying the book is an analogy 
with the well-established classical electrodynamics. In contrast to the latter, 
‘nano-electrodynamics’ is still in its infancy and far from well established. 
When a structure is miniaturized as a device, it is essential to have control 
over energy excitation and dissipation. Otherwise, when a device is squeezed 
down beyond a certain size and the energy dissipation becomes overwhelmed, 
the device will eventually collapse. It is our aim in this book to provide some 
thoughts on the task of making devices out of small structures. 

The book is partly the result of the 1999 Tsukuba symposium on the 
interaction of electrons and electromagnetic fields. I wish to thank all the 
participants of the symposium for their contributions and fruitful discussion. 
I would especially like to acknowledge the great support from Dr. H. Hori, 
Dr. Z.-C. Dong, Dr. I.S. Osad’ko, Dr. A.G. Vitukhnovsky, Dr. J. Bae, 
Dr. J.-G. Hou, and Dr. N. Yamada. I also thank Dr. K. Amemiya for his 
help in compiling this book in its present form. 

I would also like to thank Dr. Claus von Klitzing and Dr. Claus Ascheron 
for their original suggestion and encouragement in preparing this book. Fi- 
nally I would like to thank Adelheid Duhm at Springer for checking all the 
manuscripts. 



Tsukuba, Japan, 
August 2002 



Hitoshi Nejo 




Contents 



1 Introduction: Electron and Photon Systems 

H. Nejo and H. Hori 1 

I. 1 Size Considerations Concerning Electromagnetic Fields 3 

1.2 Comparison of Plasmon and Exciton 4 

1.3 Coupling of Light with Plasmons 5 

1.4 Photon Emission from STM 5 

1.5 Photon Emission from a Molecule Under an STM Tip 9 

1.6 Single Polariton Probe 11 

1.7 Coulomb Blockade and Connection 

to a Nonequilibrium Reservoir 13 

1.8 Structure of the Book 15 

References 15 

2 STM-Induced Photon Emission from Single Molecules 

Z.-C. Dong, T. Ohgi, D. Fujita, H. Nejo, S. Yokoyama, T. Terui, 

S. Mashiko, T. Okamoto 21 

2.1 Past and Present 23 

2.1.1 Semiconductor Surfaces 23 

2.1.2 Metal Surfaces 24 

2.1.3 Mechanisms 25 

2.1.4 Factors Affecting STM-induced Photon Emission 

and Modes of Measurements 27 

2.1.5 The Debate on Mechanisms 28 

2.1.6 Molecules on Surfaces 29 

2.2 STM-Induced Light Emission from Porphyrin Molecules 32 

2.2.1 Experimental Setup 32 

2.2.2 Porphyrin Molecules 33 

2.2.3 Tunneling Transport Through Molecules 

and STM Imaging 35 

2.2.4 STM-Induced Photon Emission 38 

2.2.5 Intensity and Quantum Efficiency 39 

2.2.6 Optical Spectra 41 

2.3 Conclusion 44 

References 45 




VIII Contents 



3 Photon Counting Methods in STM and SMS 

A.G. Vitukhnovsky, I.S. Osad’ko 49 

3.1 Introduction 49 

3.2 Light Emission from STM 50 

3.3 Role of Photon Counting in Time-Resolved Spectroscopy 53 

3.4 Ensembles and Single-Molecule Spectroscopy Studies 54 

3.5 Importance of Ensembles for Applications. 

Novel Ideas with Nanosize Ensembles 55 

3.6 Nanoensembles 56 

3.6.1 Regular Dendrimers 56 

3.6.2 Hyperbranched Polymers 57 

3.6.3 J-aggregates 58 

3.7 Time-Correlated Single Photon Counting (TCSPC) Method 58 

3.8 Application of Photon Counting to Kinetic Measurements 

on Nanoensembles 60 

3.8.1 Dendrimer Rotation and Viscosity 60 

3.8.2 J-aggregate Radiative Lifetime and Structure 63 

3.9 Photon Counting in SMS 64 

3.10 Intermittency in Single-Molecule Fluorescence 67 

3.11 Spectroscopy of Single Molecule with Triplet Level. 

Photon Bunching 69 

3.12 Two-Photon Correlators 72 

3.13 Full Two-Photon Correlator 

for Single Molecule with Triplet Level 75 

3.14 Tunneling Dynamics of Polymers Probed by SMS 77 

3.15 Relation Between One- and Two-Photon Counting Methods 

in a More Complex Case 80 

3.16 Local Dynamics of Individual Molecules in Polymers 82 

3.17 Conclusions 84 

3.17.1 Fast Photon Counting vs. CW Measurements 84 

3.17.2 Combination of TCSPC and NSOM: New Advantages ... 85 

3.17.3 Kinetics of Energy Transfer Between the NSOM Tip 

and Single Molecules/ Ensembles on Surfaces 85 

References 86 

4 Scanning Tunneling Spectroscopy 

and Electronic Properties of Single Fullerene Molecules 

J.G. Hou, B. Li 91 

4.1 Scanning Tunneling Spectroscopy 92 

4.2 Electronic Structures of Fullerene Molecules 94 

4.3 Scanning Tunneling Spectroscopy of Single Ceo Molecules 

Adsorbed on a Si Surface 97 

4.3.1 Ceo Molecules on an Si(lll)-(7x7) Surface 97 

4.3.2 Ceo Molecules on an Si(100)-(2xl) Surface 105 




Contents 



IX 



4.4 Electronic Transport Properties of a Single Fullerene Molecule . . . 107 

4.4.1 Electrical Resistance of a Single Ceo Molecule 

Determined by STM 108 

4.4.2 Resonant Tunneling of Single Ceo Molecules 

in a Double Barrier Tunneling Junction (DBTJ) Ill 

4.4.3 Molecular Device Using a Single Ceo Molecule 

as the Functional Unit 114 

4.5 Conclusion 117 

References 117 

5 Interaction Between Nonrelativistic Electrons 
and Optical Evanescent Waves 

J. Bae, R. Ishikawa, K. Mizuno 121 

5.1 Microgap Interaction Circuits 122 

5.1.1 Energy and Momentum Conservation 122 

5.1.2 Transition Rates of Electrons 124 

5.2 Metal Microslit 127 

5.2.1 Near-Field Distributions 127 

5.2.2 Optimum Slit Width 129 

5.2.3 Phase Matching Condition 131 

5.2.4 Interaction Space 131 

5.3 Modulation with Infrared Evanescent Waves 133 

5.3.1 Experimental Setup 133 

5.3.2 Electron Energy Spectrum 134 

5.3.3 Laser Field Dependence 136 

5.3.4 Electron Velocity Dependence 138 

5.4 Interaction with Visible Light 139 

5.5 Conclusion 141 

References 141 

6 The Tunneling Time Problem Revisited 

N. Yamada 143 

6.1 What is the Issue? 143 

6.2 Well-known Approaches 146 

6.2.1 Notation 146 

6.2.2 Following the Motion of a Wave Packet 148 

6.2.3 Velocity Under the Barrier 149 

6.2.4 The Larmor Clock 150 

6.2.5 Time-Modulated Barriers 152 

6.2.6 Dwell Time and its Decomposition 154 

6.2.7 Feynman Paths Approach 157 

6.2.8 Systematic Projector Approach 159 

6.2.9 Weak Measurement Approach 161 

6.2.10 Bohm Trajectory Approach 163 




X 



Contents 



6.3 Tunneling Times as Interfering Alternatives: 

A Feynman Paths Approach 165 

6.3.1 Lesson from Two-Slit Experiments 165 

6.3.2 Tunneling Time Probability Distribution is not Definable. 167 

6.3.3 Relation to Gell-Mann and Hartle’s Formulation 

of Quantum Mechanics 168 

6.3.4 Td = PtTt + PrTr Is Doubtful 168 

6.3.5 Range of Tunneling Times 169 

6.3.6 An Example: The Monochromatic Case 171 

6.4 Conclusion 173 

References 174 

Index 177 




List of Contributors 



Jongsuck Bae 

Research Institute 
of Communication 
Tohoku University 
2-1-1 Katahira, Aoba-ku 
Sendai 980-8577, Japan 
baeOriec . tohoku .ac.jp 

Zhen-Chao Dong 

National Institute 
for Materials Science 
1-2-1 Sengen, Tsukuba 
Ibaraki 305-0047, Japan 
DONG . Zhen-ChaoOnims .go.jp 

Daisuke Fujita 

National Institute 
for Materials Science 
1-2-1 Sengen, Tsukuba 
Ibaraki 305-0047, Japan 
Fuj ita . Daisuke@nims .go.jp 

Hirokazu Hori 

Jamanashi University 
4-3-11 Takeda 
Kofu 400-8511, Japan 
hirohori@es . yamanashi .ac.jp 

J.G. Hou 

Structure Research Laboratory 

and Open Laboratory of Bond 

Selective Chemistry 

University of Science and Technology 

of China 

Hefei 230026 



P. R. China 
j ghou0us t c . edu . cn 

Ryo Ishikawa 

Research Institute 
of Communication 
Tohoku University 
2-1-1 Katahira, Aoba-ku 
Sendai 980-8577, Japan 
issi^riec . tohoku .ac.jp 

Bin Li 

Open Laboratory 

of Bond Selective Chemistry 

University of Science and Technology 

of China 

Hefei 230026 

P. R. China 

libin@mail . ustc . edu . cn 

Shinro Mashiko 

Communication Research 
Laboratory 

Kobe, Hyogo 651-2401, Japan 
masahikoOcrl .go.jp 

Koji Mizuno 

Research Institute 
of Communication 
Tohoku University 
2-1-1 Katahira, Aoba-ku 
Sendai 980-8577, Japan 
ko j iOriec . tohoku .ac.jp 




XII List of Contributors 



Hitoshi Nejo 

National Institute 
for Materials Science 
1-2-1 Sengen, Tsukuba 
Ibaraki 305-0047, Japan 
NEJO . HitoshiOnims .go.jp 

Taizo Ohgi 

National Institute 
for Materials Science 
1-2-1 Sengen, Tsukuba 
Ibaraki 305-0047, Japan 
OHGI . Taizo@nims .go.jp 

Takayuki Okamoto 

RIKEN (Institute of Physical 
and Chemical Research) 

Wako, Saitama 351-0198, Japan 
okamot o^postman . riken .go.jp 

I.S. Osad’ko 

Lebedev Physical Institute RAS 
Leninsky Prospect, 53 
Moscow 119991, Russia 
osadko@sci . lebedev . ru 



Toshifumi Terui 

Communication Research 
Laboratory 

Kobe, Hyogo 651-2401, Japan 
terui@crl.go.jp 

Alexei G. Vitukhnovsky 

Lebedev Physical Institute RAS 
Leninsky Prospect, 53 
Moscow 119991, Russia 
alexeiOsci . lebedev . ru 

Shiyoshi Yokoyama 

Communication Research 
Laboratory 

Kobe, Hyogo 651-2401, Japan 
syoko@crl.go.jp 

Norifumi Yamada 

Fukui University 

3-9-1 Bunkyo 

Fukui 910-8507, Japan 

f yamada@ilnwsl . f uis . f ukui-u .ac.jp 




1 Introduction: Electron and Photon Systems 



H. Nejo and H. Hori 



While previous books about nanometer scale phenomena have dealt with 
either electronic states or electromagnetic fields, this book is concerned with 
both as a unified whole. The characteristics of electrons or electromagnetic 
fields apparent at nanometer scales are quite different from those at the 
macroscopic scale. Nowadays these areas are called nanoscience or nanotech- 
nology. They attract the interest of a wide range of people for their applica- 
tions to new branches of technology [1-15]. 

When we deal with electromagnetic fields at the nanometer scale, the 
scale of interest is much shorter than the wavelength of light, at approxi- 
mately 1 pm. The characteristics of such electromagnetic fields are therefore 
quite different from those pertaining to macroscopic sizes. The former are 
characterized as localized waves, while the latter are characterized as propa- 
gating waves. The former nanoscale wave is usually called an evanescent wave 
and its characteristics have now been studied by many researchers [16-50]. 

For electron waves, the situation is similar to that of a propagating light 
wave in that an electron can propagate through crystalline metals and semi- 
conductors, whereas on nanoscales, when a metal block is cut in half, the 
propagating electron wave starts to penetrate into the vacuum from the sur- 
face. The nature of this penetrating wave can be very similar to that of an 
evanescent light wave. The electron wave penetrating the vacuum appears as 
a tunneling phenomenon: the electron can pass through a potential barrier. 
These features have been studied using the scanning tunneling microscope 
(STM) [51]. 

Consider the relationship between an electron and an electromagnetic 
field. Electromagnetism shows that the moving charged particle causes an 
energy flow, called a Poynting vector. Thus the field of the electron and 
the electromagnetic field are closely related. This interaction has been well 
studied at the macroscopic scale but not at the nanoscale. 

One can see the interaction of an electron and the electromagnetic field in 
the case of a Coulomb blockade [52-62] . A Coulomb blockade can be described 
as follows. A double tunnel junction is biased by a battery. When the battery 
does work and the Fermi level of the emitter reaches the chemical potential 
of the central electrode of the double tunnel junction, electrons can tunnel 
through the first tunnel barrier by overcoming the Coulomb energy e^(7/2, 
where C is the macroscopic capacitance of the tunnel junction. When an 
electron tunnels through the potential barrier, the Poynting vector crosses 
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a resistor that is connected to the double tunnel junction, indicating that 
current flows through the junction. 

While the interaction between electrons and electromagnetic flelds is a 
most interesting phenomenon, it has not yet been studied properly at the 
nanometer scale. The reason for this can be entirely attributed to the lack of 
technological development both in realizing nanometer scale structures and 
also in connecting them to excitation sources and detectors. For example, 
the interaction between a penetrating electron wave, i.e., a tunneling elec- 
tron, and an electromagnetic held has been studied. A planar tunnel junc- 
tion was prepared and then irradiated by light. Since a tunneling current was 
detected, the existence of a penetrating or decaying electron wave was de- 
duced. However, as both the electrode and the tunnel barrier were irradiated, 
it was not possible to extract the net interaction effect between the tunneling 
electron and the evanescent held from the interaction between the propagat- 
ing electron wave and the electromagnetic field in the electrodes. Due to the 
experimental set up, the details of the interaction were not revealed. 

A similar situation occurred in an experiment involving the interaction 
between a propagating electron and an evanescent held, although in this case 
not involving a tunneling electron (Chap. 5). Light was introduced into an 
insulating slab through which the propagating electron flowed. Since the na- 
ture of the electromagnetic held in an insulator is evanescent, the propagating 
electron was assumed to interact with the evanescent wave in the slab. Light 
emission was observed after the electron reached a non-fluorescent screen. 
This is thought to be a consequence of the interaction between the propagat- 
ing electron and the evanescent field. However, as the propagating electron 
interacts with the material of the slab itself, one cannot distinguish the inter- 
action of the electron and the evanescent field from the interaction between 
the electron and the matter. 

These examples show that it is no trivial task to set up an experiment that 
distinguishes the interaction between a tunneling electron and an evanescent 
field from the other possibilities. In this book, the reader will discover some 
efforts in this direction. Owing to the recent development of nanoscience and 
nanotechnology, we now possess a variety of tools for realizing experimental 
setups that can probe the interaction between a tunneling electron and an 
evanescent held. 

The purpose of this book is to show the importance of constructing 
nanoscale electric and electromagnetic flelds connected to a reservoir. We 
start with a conventional electric or electromagnetic held treatment, for which 
many theoretical and experimental results have been accumulated, and then 
move on to the importance of the reservoir connection. 
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Fig. 1.1. Relation between the size of matter and the extension of the electromag- 
netic field that can be detected. Cavity QED is described in Sect. 1.2, SNOM in 
Sect. 1.3, photon emission from STM in Sect. 1.4, elementary excitation STM in 
Sect. 1.5, single polariton probes in Sect. 1.6 and mode connection to the macro- 
scopic mode is described in Sect. 1.7 



1.1 Size Considerations 
Concerning Electromagnetic Fields 

Size considerations of electromagnetic fields are illustrated in Fig. 1.1 [63]. 
Cavity quantum electrodynamics (QED) is described in Sect. 1.2, scanning 
near-field optical microscopy (SNOM) in Sect. 1.3, photon emission from STM 
in Sect. 1.4, elementary excitation STM in Sect. 1.5, single polariton probes 
in Sect. 1.6 and mode connection to the macroscopic mode in Sect. 1.7. 

Single-electron tunneling is accompanied by the creation of a quantum 
Uuj of electromagnetic field in accordance with the following energy balance: 
eV — Eo + ^o;. Whenever the microscopic electronic state changes, the energy 
changes by the unit quantum from the overall mode function. When the 
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energy transfer between the electronic state and the electromagnetic state is 
described in this manner, it is assumed that one is treating a mode function 
whose size is three times the light wavelength. Such a picture is applicable 
to a laser cavity, where the size of confinement is of the same order as the 
wavelength. We are concerned with electromagnetic fields of much smaller 
size. When a molecule is placed under an STM tip, the gap size between the 
tip and the substrate is of the order of a few nanometers. The electromagnetic 
field in this nanometer-sized gap is much smaller than the conventional size 
described by the mode function. Such an electromagnetic field is described 
by the combination of the electrons in the metal electrode (tip) and the mode 
function on the larger scale, the so-called dressed photon. On an even smaller 
scale, the effect of radiation on a tunneling electron was studied by Tien and 
Gordon in the early 1960s and described as the side-band effect of the initial 
and final electronic states [64,65]. The total Hamiltonian H is then described 
by H = Hq eV cosa;^, where Hq is the Hamiltonian excluding the effect of 
the electromagnetic field. Johannson and Xiao discussed ways of describing 
the electromagnetic field. There is a further complication for structures in 
an electromagnetic field at nanometer scales. A dressed photon has a large 
k vector which does not follow the dispersion relation u;/c= \k\. Due to this 
large k vector, there is the possibility of transition from the initial k vector 
k{ to the final k vector fef via various states interacting with the k vector of 
the electron. 

1.2 Comparison of Plasmon and Exciton 

The field of cavity QED has been well studied [66-82]. The behavior of a 
quantum- well exciton resonantly coupled to a single-mode microcavity field in 
a linear (low-excitation) regime can be modeled in terms of coupled harmonic 
oscillators between which the energy is coherently recycling back and forth. 
When an exciton is coupled with a polariton, the state is described as an 
exciton-polariton. When a plasmon is coupled with a polariton, it is described 
as a plasmon-polariton. This is well defined when the size under consideration 
is at most half a wavelength. Below half a wavelength, the situation changes 
drastically. At such small scales, the mode is no longer defined. This is the 
situation when we try to see an elementary excitation or a single polariton. 
The difficulty can be brought out using STM as described in the next section. 
Electromagnetic fields in the nanometer region can be measured using light 
emission from a molecule under an STM tip excited by tunneling electrons. 
For any measurement a probe is required to detect signals. STM uses a probe 
that is very close to the sample. The tip-sample distance is usually only a 
few nanometers so that proximity effects always affect the electromagnetic 
field. Light emission from a molecule under an STM tip has been detected 
by two methods and the differences induced by different tip materials have 
been compared using tungsten and ITO-coated fiber tips (Chap. 2). 
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Fig. 1.2. Dispersion relations of an electron {left) and a photon (right). Photon 
states coupled with matter are denoted by horizontal lines. These have a larger k 
vector than those which obey the dispersion relation. SNOM detects these large k 
vector states 



1.3 Coupling of Light with Plasmons 

Figure 1.2 shows the dispersion relation of electrons and photons. SNOM 
uses this large k photon state. The effect between light and an insulator has 
been well studied. Ebbesen has also demonstrated an effect between light 
and metal. When a two-dimensional array on a metal film is irradiated with 
light, transmission efficiency can exceed unity when normalized to the area 
of the holes. This effect has also been studied by SNOM using a metallic tip 
or particle as scatterer [83-90]. 

1.4 Photon Emission from STM 

Photon emission from an STM is explained by dipole moment radiation due 
to a tip-induced surface plasmon (Fig. 1.3) [91-105]. The electron-photon in- 
teraction can be considered by borrowing the idea of virtual photons. When 
an electron with k vector fee tunnels through a potential barrier, it can inter- 
act with the potential surface and change its state to When the tunneling 
electron then experiences inelastic tunneling, a photon with k vector fep is 
generated and the electron leaves the potential barrier. Since the photon is 
a virtual photon, it can interact with matter to borrow k vector and change 
its state to k vector fcp (Fig. 1.4) [106]. 

This electron-photon interaction in a potential barrier in real space is 
illustrated in Fig. 1.5. When we prepare an electron path with a potential 
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Potential barrier 




Fig. 1.3. Schematic diagram showing the tip-induced surface plasmon. When an 
STM tip is biased, a surface plasmon is induced between the STM tip and the 
metal surface. The wave vector k is not conserved, since the periodicity in the 
lateral direction is broken. This surface plasmon is the source of the dipole moment 
radiation 




Potential surface 



Fig. 1.4. Schematic diagram of the interaction between an electron and a photon on 
a potential surface. When an electron with k vector ke interacts with the potential, 
it may borrow virtual k to become k^. This electron interacts with another potential 
and emits a virtual photon with k vector kp. It then returns to a real electron. The 
virtual photon kp interacts with the potential and recovers to a real photon kp 
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Fig. 1.5. Schematic diagram of the interaction between an electron and a photon 
in a potential barrier. When an electron ^ tunnels through a potential barrier it 
may interact with a virtual photon which is the electromagnetic field E coupled 
with matter 



Photon source 





Fig. 1.6. Experimental setup showing the electron-photon interaction based on 
the possibility of interaction in a potential barrier. The evanescent field is intro- 
duced into the tunnel gap through which the electron tunnels. The electron-photon 
interaction is detected as a matter-dressed photon by a near-field photon detector 
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barrier in the middle {x direction), it is possible to introduce a path for the 
photon which is perpendicular to the electron path {y direction) . The process 
described above can occur in such a potential barrier. 

Photon-assisted tunneling was studied by Tien and Gordon in the 1960s. 
At that time they used a planar junction. Later they illuminated an STM tip 
but the state of the photon in the tunnel gap was not well defined. As shown 
in Fig. 1.6, by using a metal-coated fiber tip, if an evanescent wave is intro- 
duced into the tunnel gap, the interaction of the tunneling electron with the 
evanescent wave is well defined. This interaction can be detected as current 
I{uo) and as a matter-dressed photon by a near-held photon detector. Using 
this setup, the lateral size of the interaction is restricted by the diameter of 
the aperture of the hber tip. Again, in this setup, the evanescent wave can 
excite the substrate surface and the net effect must be distinguished from 
this photon-substrate interaction. 




Fig. 1.7. An experimental setup showing the interaction of a tunneling electron 
with a photon. Two metal-coated hber tips face each other. When they come close 
enough, an electron tunnels through the potential barrier. A dressed photon tunnels 
across the gap. The interaction of the electron states, such as momentum p, angular 
momentum I and spin s, with the photon can be detected by a photon detector 
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It should be possible to detect tunneling of momentum, angular momen- 
tum and electron spin using two metal-coated fiber tips (Fig. 1.7). When two 
such tips are placed face to face in close proximity, electrons can tunnel be- 
tween the two metallic electrodes. If light is introduced at one fiber tip, an 
evanescent field is generated between the two fiber tips. This heavy photon, 
dressed with electron momentum, angular momentum and spin, can be de- 
tected by another fiber tip connected to a photon detector. Gimzewski et 
al. pioneered the study of photon emission from an STM tip and sample in 
1988 [107-116]. When a bias voltage is applied across a tunneling barrier 
so that an electron tunnels, the electron in the counter electrode has an ex- 
cess potential energy equal to eV where e is the electronic charge and V the 
bias voltage. The electron will relax by releasing this excess energy. Most of 
the excess energy is converted into heat by collision with phonons, but some 
fraction of the excess energy will be emitted as a photon under appropriate 
conditions. If the bias voltage is a few volts, the emitted photon has energy 
corresponding to the visible region. When one observes such photon emission 
from a pristine metal surface, the photon emission from an STM is consid- 
ered as dipole moment radiation due to the tip-induced surface plasmon. In 
the surface plasmon view, photon emission from an STM can be described 
as follows [117-120]. When the tip-sample separation is no more than 0.5- 
1 nm, the surface plasmons on the two surfaces interact strongly and form 
an interface-plasmon mode. The resonance in the light emission spectrum 
occurs at the frequency for which half a wavelength of the interface mode 
fits into the cavity between tip and sample (although later, the reader will 
discover that in reality this does not apply to such a small system). The 
charge oscillations associated with the interface plasmon have opposite signs 
to the two electrodes. Photon emission from an STM as described via the 
surface-plasmon-polariton mode will be described in the next section. 

1.5 Photon Emission from a Molecule 
Under an STM Tip 

To observe photon emission, the object of interest should be connected to 
some macroscopic mode [121-131]. Here the object of interest is a single 
molecule, a typical quantum mechanical object [132]. When the molecule is 
excited by a tunneling electron, the dynamic dipole moment couples with 
the surface plasmon between the tip and the substrate. This surface plasmon 
is localized between the tip apex and the counter substrate surface. It is 
a coupled mode of the electromagnetic field in the cavity (formed by the 
tip and the substrate surface) and the electrical system of the tip and the 
substrate (and the molecule itself) [133,134]. This surface plasmon propagates 
along the surface of the tip and/or the substrate surface and finally connects 
with the macroscopic mode. The size of the SNOM electromagnetic field is 
a few tens of nanometers between the tip and sample, since the tip height is 
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Fig. 1.8. Schematic diagram showing the evanescent field around a molecule. When 
a molecule is excited by a tunneling electron, it subsequently generates an evanes- 
cent field as it relaxes. This evanescent field exists only around the molecule whereas 
the electromagnetic field extends over the entire gap between the tip and the sub- 
strate. The diagram also compares evanescent field detection by a fiber tip and 
dipole moment radiation 



regulated using force between the tip and sample. While the lateral extent of 
the electromagnetic field of photon emission from STM is also some tens of 
nanometers, as the extension of the plasmon is great along a metal surface, it 
is possible to considerably reduce the vertical size of the electromagnetic field 
by putting a molecule in the tunnel gap (Fig. 1.8) [135-142]. The molecule 
can be excited by a tunneling electron by holding a metal-coated fiber tip 
above it. 

The eflFect of scattering a photon with a large k vector can be verified 
using two methods. In the first method, a molecule is placed on a metal 
substrate and a standard W STM tip is placed above the molecule (Fig. 1.9a). 
The molecule is excited by a tunneling electron and emits light when the 
molecule relaxes. A photon with a large k vector is scattered by the tip apex 
and turns into a photon with a smaller k vector that subsequently reaches a 
photon detector. 

The second method uses an indium-tin-oxide (ITO) coated fiber tip [143] 
supported above a molecule which is excited by a tunneling electron. When 
this molecule emits light with a large k vector, in the ITO-coated fiber tip 
case, a photon with a large k vector is scattered and introduced into the fiber. 
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Fig. 1.9. A comparison of near-held detection using (a) a W tip, (b) a hber tip. 
The evanescent held with a large k vector is scattered by the metal tip apex and 
evolves into a small k vector which subsequently reaches a photon detector (a). 
The evanescent held induces dipoles in the hber tip and these dipoles propagate in 
the hber tip to reach a photon detector (b) 



In other words, a dipole is induced in the hber and these dipoles propagate 
along the hber. This dipole is also induced in the ITO layer, since ITO acts as 
an insulator for higher frequencies. A suitable experimental setup for studying 
an electromagnetic held at the nanometer scale is to study light emission from 
a molecule placed under an STM tip. In this system, the existence of a dipole 
moment is well dehned. The dipole moment is not static but dynamic, and is 
responsible for the light emission. The direction of the dipole moment may 
then be manipulated parallel or perpendicular to the substrate surface. When 
the dipole moment is perpendicular to the surface, i.e., parallel to the tip- 
substrate direction, strong coupling with the plasmon results between the tip 
and the substrate surface and this plasmon turns into a radiative held. In 
contrast, when the dipole is parallel to the substrate surface, the coupling 
with the plasmon mode is weak so that the excited dipole is transferred to a 
neighboring surface molecule if there are many such molecules on the surface. 
This is a realistic situation when molecules are evaporated onto the surface. 



1.6 Single Polariton Probe 

As an example of a single polariton probe, consider the following experimen- 
tal setup (Fig. 1.10) [144]. Two molecules are placed in nano-gaps between 
electrodes separated by a distance much smaller than an optical wavelength. 
Let us assume that these molecules are excited by a coherent electron wave 
with different phase. Then the phase of the electromagnetic field radiated 
from each molecule will be different so that the electromagnetic field cancels 
out at the far field. 
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Cancel at far field 

Fig. 1.10. Schematic view of prohibited photon emission. Two molecules are put 
in the gaps of the electric circuit and excited by a coherent electron. If the phase 
of the dipole moment of each molecule in the gap differs by tt , the electromagnetic 
field cancels at the far field 
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Fig. 1.11. Schematic view of allowed radiation. When a near-field probe is brought 
to one of the molecules so that it absorbs radiation from one side of the molecule, 
the electromagnetic field no longer cancels at the far field and the system emits 
light 
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Phase shifter 



Fig. 1.12. Schematic view of a setup for observing photon emission. When the 
magnetic field is applied perpendicularly to this system, the phase of the electron 
is changed along a path shown in the drawing as a function of the magnetic field. 
The photon emission is thus modulated as a function of the magnetic field 




Let us now imagine that we can place a probe above one of the molecules 
and thereby pick up a polariton (Fig. 1.11). As a consequence, the electro- 
magnetic field radiated from the other molecule can reach the far field and 
this electromagnetic field can be observed. 

This gedanken experiment can be extended. We consider a magnetic field 
applied vertically to the system (Fig. 1.12). By considering a loop along the 
electrode, the phase of the electron wave can be changed using the Aharanov- 
Bohm effect. The phase of two electromagnetic fields generated from each 
molecule can then be exchanged. The amplitude of the total electromagnetic 
field at the far field can thus be varied as a function of the magnetic field. 

1.7 Coulomb Blockade 

and Connection to a Nonequilibrium Reservoir 

In this section, we consider the meaning of the term ‘observation’, whereby 
the mode of an electromagnetic field as discussed above is connected to the 
macroscopic mode. 

A Coulomb blockade should be understood in terms of the connection be- 
tween a nanostructure and a macroscopic reservoir. Consider a double tunnel 
junction connected to a battery and an Ampermeter (Fig. 1.13). The system 
consists of two loops. The first loop includes the emitter, the central electrode 



14 



H. Nejo and H. Hori 



Electron exchanges in the central 
electrode 



Central 

electrode 




When an electron goes through the first junction emitting 
a photon, another electron goes through a resistor. 

Which Is faster? 

Fig. 1.13. Schematic view of a Coulomb blockade and electromagnetic field. When 
an electron tunnels through the first tunnel junction, the Poynting vector crosses 
the first resistor and this electron cannot be detected. When the first electron tun- 
nels, another electron tunnels simultaneously through the second tunnel junction. 
Assuming that a photon is emitted during the first electron tunneling, the photon 
is detected by a photon detector 



and a resistor connected to the central electrode via stray capacitance. The 
second loop includes the collector, the central electrode and another resis- 
tor connected to the central electrode. A Coulomb blockade arises when an 
initial electron is localized in the central electrode. A second electron can- 
not tunnel through the first potential barrier until the first electron leaves 
the central electrode. By observing the first electron, the second electron is 
blocked from tunneling. We consider the voltage range up to e/2(7 where 
the charge of the central electrode is zero (Fig. 1.14). This suggests that the 
first electron resides virtually in the central electrode and cannot be observed 
from outside. This in turn leads to the consideration that the resistor in the 
first loop is connected to the reservoir but that this resistor is not connected 
to the detector. There is a correlation between the first and second electrons. 
When the first electron tunnels through the first potential barrier and reaches 
the central electrode, the Poynting vector crosses the resistor. This Poynting 
vector crossing is equated with observation of the electron. Put another way, 
when the first electron tunnels through the first potential barrier, the elec- 
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Fig. 1.14. Schematic view of a Coulomb staircase where the first electron resides 
virtually in the central electrode. The Poynting vector crosses a resistor that is 
connected to a reservoir but not to a detector so that we cannot observe the first 
electron in the Coulomb gap 



tron which stayed in the second electrode should leave the central electrode 
simultaneously and this first electron passes through the resistor. According 
to this scheme, a Coulomb blockade is closely connected to the process of 
observation. 

1.8 Structure of the Book 

Chapter 2 describes photon emission experiments from single molecules using 
STM. The mechanism of single molecule photon emission is further described 
in Chap. 3. Fullerene molecules are commonly used for realizing nanoscale 
structures, and a detailed study of fullerene molecules is reported in Chap. 4. 
Chapter 5 describes the interaction between electrons and photons as demon- 
strated experimentally, whilst Chap. 6 provides the theoretical treatment of 
electron tunneling. 
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The investigation of single quantum structures lies at the frontier of 
modern chemistry, physics, and nanoelectronics. A single molecule is a sim- 
ple quantum system for studying electron-electron and electron-photon in- 
teractions and other quantum mechanical effects. The quantum nature of 
molecular systems makes single-molecule studies important not only for a 
fundamental understanding of basic phenomena underlying bulk molecular 
sample behavior, such as in organic electronics, but also for shaping emerg- 
ing concepts and principles of molecular nanotechnology. For instance, over 
the last decade, organic thin film science has been developing at an un- 
precedented pace with attractive technological applications in organic light- 
emitting diodes, dye lasers, organic transistors, liquid-crystal displays, and 
photovoltaic cells. However, even with these technological breakthroughs, 
there remain deep questions about the underlying physics of charge and ex- 
citon transport, energy absorption-transfer-decay processes, optical interac- 
tions, and the nature of electronic states within these bulk molecular materi- 
als. Molecular-scale investigations will clarify the fundamental mechanisms of 
these phenomena and processes and thus offer guidelines for designing better 
molecular materials and devices. 

Single molecular studies are also driven by the emerging field of molecular 
scale electronics, which aims to develop electronic components such as wires, 
switches, amplifiers, and transistors via the use of just a few molecules or 
even a single supramolecular assembly, and to integrate these into a device. 
The field involves many challenging unsolved issues such as the electrical ad- 
dressing of molecules and circuit design (e.g., input/output interconnects and 
clock frequency). We are still at the materials stage, seeking to understand at 
a fundamental level what is going on in molecules by measuring their trans- 
port and optoelectronic characteristics. To date, most transport phenomena 
measured seem to be of contact-type rather than a true conductivity of the 
molecule itself, except for those related to internal motion (e.g., vibronic and 
rotational transitions). Scanning tunneling microscopy provides controlled 
two-terminal measurements, and its development together with other tech- 
niques should offer new experimental approaches for probing electron-photon 
interactions in, and electron transport through individual molecules. Even- 
tually, this may lead to the determination of true molecular conductance, a 
measure of the electronic transparency of a single molecule. 
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Single-molecule studies comprise three aspects; microscopy in the space 
domain, spectroscopy in the frequency domain, and fluorescence decay in the 
time domain. To understand single-molecule behavior in a nano-environment, 
it is necessary to describe the electromagnetic radiation in the nanoscopic sys- 
tem. Scanning probe techniques (e.g., scanning tunneling microscopy) allow 
study of individual molecules on surfaces, while optical techniques (e.g., ultra- 
sensitive fluorescence spectroscopy) enable characterization of single-molecule 
behavior in complex condensed environments. Single- molecule measurements 
make it possible to unravel the normal ensemble averaging that occurs when 
a large number of molecules are probed simultaneously. This allows direct ob- 
servation of the true distribution and time-resolved information about a given 
property (structure and functions) of a molecule, or many single molecules 
but on a one-by-one basis [1]. 

The invention of the scanning tunneling microscope (STM) has revolution- 
ized our way of dealing with individual atoms and molecules. It can provide 
atomically resolved information on the conformation and adsorption nature 
of individual molecules on a surface. Yet it can do more than just observe the 
nanoworld. It can also be used to manipulate atoms and molecules and to 
construct supramolecular assemblies and functional molecular machines. Fur- 
thermore, combined STM and optical techniques are particularly attractive 
for single- molecule studies. It becomes possible not only to image individ- 
ual molecules on a surface, but also, through probing photons emitted via 
inelastic processes, to gain insights into the tunneling regime, e.g., specific 
characteristics of diverse molecule and substrate excitations as well as modi- 
fied optical properties due to quantum confinement [2-4]. 

Inelastic processes excited by tunneling electrons carry chemically specific 
information on molecular systems, but are very difficult to detect by conven- 
tional scanning tunneling spectroscopy alone due to their small contribution 
to the total tunneling current [5]. Tunneling of electrons or holes from an 
STM tip provides an extremely localized source of low-energy carriers (only 
a few eV), which offers resolution down to the subnanometer scale for photon 
emission mapping. Spatial mapping of optical signals may enable us to add 
‘color’ to STM images and perform chemical mapping on surfaces. Such high 
resolution spectral photon maps will give access to the electromagnetic cou- 
pling between the tip and sample and its variation even on the atomic scale. 
Another advantage of STM is its freedom to vary experimental parameters 
on well-characterized surfaces, thus permitting detailed investigation of the 
processes involved. 

The atomic resolution of a scanning tunneling microscope arises from 
the extremely high sensitivity of tunneling currents to the small tail overlap 
between the wavefunctions of both tip and substrate across the tunneling 
barrier. How does the molecule couple with the electrodes and what happens 
optically when a single molecule is inserted into this subnanometer gap and 
the molecular junction is biased (Fig. 2.1)? This is the issue we shall address 
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Fig. 2.1. Schematic view of a molecule sandwiched between tip and substrate. Light 
is emitted due to excitation by tunneling electrons 



in the present chapter. Following a brief review of STM-induced photon emis- 
sion from metal and semiconductor surfaces, we shall report our research on 
tunneling-electron-induced light emission from porphyrin molecules. General 
concepts and principles involved in STM-induced photon emission phenom- 
ena, experimental setup, factors affecting photon emission, light emission 
mechanisms, and current challenges in the field will be discussed. 

2.1 Past and Present 

Before we proceed to the discussion of STM-induced photon emission from 
single molecules, we shall briefly review the field of light emission induced by 
tunneling electrons on semiconductor and metal surfaces [2,4,6]. 

2.1.1 Semiconductor Surfaces 

Photon emission stimulated by tunneling electrons can be traced back to 
Lambe and McCarthy in 1976 on metal-oxide-metal solid-state junctions 
[7]. The first experimental evidence of STM-induced light emission was re- 
ported by Gimzewski et al. from Ta and Si(lll) surfaces [8] and soon after- 
wards, enhanced light emission from Ag films was also observed [9]. These 
results showed similarities to conventional inverse photoemission [10] with 
additional spectral features assigned to optical transitions from field-induced 
states. STM is also used for the investigation of local cathodoluminescence 
(CL) in GaAs/Gai_a:Ala;As heterostructures [11-13] and CdS surfaces [14] 
in which emission arises from optical transitions between band edges. The 
polarization effects of STM-induced photon emission have also been reported 
for GaAs [15] and Si(lll) [16]. There are also reports using a conductive and 
transparent indium-tin-oxide (ITO) tip for both carrier injection and photon 
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collection [17], the so-called near-field detection (to be contrasted with far- 
field detection using a lens system located far away from the light source). 
Recent progress in the field has provided photon maps with spatial resolution 
down to the atomic scale for Si(lll) [18] and hydrogen-passivated Si(lOO) [19] 
surfaces. 

2.1.2 Metal Surfaces 

On the other hand, research on STM-induced photon emission from metal 
surfaces has also received a great deal of attention because the mechanism in- 
volved in the emission process is new and specific to the nanometer proximity 
and electromagnetic coupling of the tip to the surface [2,4,6]. An energetic 
electron impinging on a metal surface can generate electromagnetic radiation 
via dipolar resonance mechanisms such as transition radiation and surface 
plasmon excitation. An STM tip can be used as an electron source for gen- 
erating fluorescence on metal surfaces in both the far-field emission regime 
and regimes of proximity field emission and tunneling. At high bias volt- 
ages (> 100 V), emission phenomena were found to be similar to those due 
to conventional electron guns [20]. However, in the proximity field-emission 
regime (< 40 V) and in the tunneling regime (< 5 V), a new and interesting 
phenomenon was observed [21]. Instead of the absence of photon emission 
by conventional electron guns for such low electron energies, intense photon 
emission occurs with redshifted optical spectra and an enhanced quantum 
efficiency up to ~ 10"^ photons/electron for a W tip on Ag. A universal in- 
tensity maximum for bias voltages between 3-4 V is observed in the tunneling 
regime on all metallic surfaces. This has been attributed to competition be- 
tween the initial increase in allowed channels for photon emisison and the 
subsequent field-enhancement weakening due to the vertical tip retraction in 
that region [22]. In an STM configuration, the tip-sample distance is in the 
nanometer region, much shorter than the wavelength of the emitted light. 
This results in localized electromagnetic coupling between tip and sample 
and hence new emission features. 

Different single-crystal metal surfaces have been shown to exhibit different 
emission characteristics [20] , but the spatial resolution of chemical signals has 
not been fully investigated. The localization of photon emission to a small 
area beneath the tip promises a form of chemical microscopy, because the 
photon signal is related to the local chemical environment, and chemical 
information can be extracted from the complicated photon emission signals 
[18,23]. Strictly speaking, the spectroscopy of emitted light is the only way 
to identify unknown metallic structures. Nevertheless, since different metals 
have different optical properties, their emission efficiencies can vary over two 
orders of magnitude and photon mapping can thus be effectively used to 
distinguish known regions on flat surfaces. Analysis of photon maps is based 
on simultaneous topographic, bias-dependent imaging in combination with 
spectroscopic data. A contrast mechanism thus formulated can help to reveal 
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various factors that determine the local efficiency of photon emission and 
to gain further insight into the physical origin of factors governing inelastic 
tunneling excitation of photons in mesoscopic systems. 

2.1.3 Mechanisms 

Theoretically, photon emission from metal surfaces was attributed to radia- 
tive decay of plasmons for solid state tunneling junctions and STM-induced 
photon emission. These plasmons are excited either by inelastic tunneling 
(lET) [24,25] or by hot-electron thermalization (HET) [26]. 

For lET excitation, a tunneling electron gives away a fraction of its en- 
ergy to a plasmon that is subsequently converted to a photon, as shown in 
Fig. 2.2a and b. Tunneling electrons directly excite the charge density oscil- 
lation between the tip and surface, i.e., the local (interface) plasmon modes. 
The local plasmons transfer their energy to photons that are then emitted. 

In HET, an electron tunnels through a barrier elastically and then loses 
its energy to the sample by thermalization (or hot-electron decay) via plas- 
mon generation. The inelastic event, which leads to photon emission, occurs 
during tunneling for lET but after tunneling within the sample for HET. 
Hot-electron injection is still debated as one of the mechanisms for photon 
emission from metal-oxide-metal junctions, whereas in STM configurations, 
inelastic tunneling excitation is believed to be dominant in the photon emis- 
sion process from metal and indirect semiconductor surfaces. Hot-electron 
excitation is negligible compared with lET excitation in terms of high inten- 
sities and lack of significant bias dependence (suggestive of a resonant process 




Fig. 2.2. Schematics of mechanisms for STM-induced photon emission: (a) excita- 
tion of localized surface plasmons (ESP) on metal surfaces via inelastic tunneling 
(lET), the so-called tip-induced plasmon (TIP) mode, (b) radiative decay via inelas- 
tic tunneling into surface states on metal and semiconductor surfaces, (c) conven- 
tional luminescence via electron-hole pair recombination through elastic tunneling 
on direct band gap semiconductor surfaces 
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during tunneling) as well as consistent fluorescence spectra between theory 
and experiment. 

Barrier height measurements show almost the same value for both elastic 
and inelastic tunneling [27]. This is somewhat unexpected because one would 
expect the inelastic barrier height to be larger than the elastic one by ~ /iz//2. 
An electron that tunnels inelastically is therefore thought to traverse most of 
the barrier elastically and undergo an inelastic event only toward the end of 
tunneling. On the other hand, for STM-induced light emission from semicon- 
ductor surfaces, in particular direct band gap materials such as p-type GaAs 
[11], hot-electron decay within the semiconductor is the primary mechanism 
for luminescence. Photon emission is due to elastic tunneling of electrons from 
the tip into the GaAs conduction band and subsequent recombination with 
holes in the valence band, i.e., electron-hole pair recombination, as shown in 
Fig. 2.2c. 

There is increasing interest in the investigation of light emission from 
surface states associated with either intrinsic surfaces or adsorbed atoms or 
molecules. This is due on the one hand to an ongoing effort to understand the 
electronic properties of a surface with or without adsorbates on a very local 
scale, and on the other hand, the lack of a well-established physical picture 
for the light emission process. The mechanism shown in Fig. 2.2b, though 
favored, is not the only one available for interpreting light emission [19]. An 
HET mechanism via surface states is also possible. 

Two models have been proposed to explain STM-induced photon emission 
from metal surfaces [4] . The model due to Persson and Baratoff [28] assumes 
that tunneling occurs from an s-like orbital at the tip apex to a spherical 
free-electron-like metallic particle. This model predicts that enhanced pho- 
ton emission occurs if tunneling electrons excite the dipole plasmon resonance 
of the metal particle, while direct transition radiation is found to be low. 
Furthermore, the dipolar plasmon excitation is predominantly caused by the 
inelastic tunneling process (with a quantum efficiency up to 10~^) in com- 
parison with hot-electron decay (10~^). Intensive photon emission has been 
interpreted as arising from the radiative decay of dipolar plasmons created 
by the close proximity of a tip to a surface. 

Another model (a fashionable electrostatic one) due to Johansson et al. 
[29] takes into account the role of a tip, viewing tunneling as occurring be- 
tween a metallic nanosphere and a planar metal surface. The classical electro- 
magnetic response of the system is calculated to obtain the strength of field 
fluctuations that are responsible for photon emission. The emission intensity 
is determined from the product of the power spectrum of the tunneling cur- 
rent and that of the enhanced electric held induced by the tip. Light emission 
is resonantly enhanced due to the formation of local coupled plasmons, the 
so-called tip-induced plasmon (TIP) modes in the cavity formed between the 
tip and sample. All modes of the electromagnetic field inside the tunneling 
gap are strongly coupled to the collective motion of the electron gas inside 
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the tip and sample, and hence, all of them should be called plasmons. These 
plasmons decay into photons outside the gap, in other words, the decoupling 
occurs far away from the metal surface in the free space where the modes 
of the free electromagnetic field can be called photons. The resonance occurs 
at the frequency at which half a wavelength of the TIP modes fits into the 
cavity between the tip and sample. Compared with light emission in vacuum, 
the enhancement could be as large as 4-5 orders of magnitude. The quantum 
efhciency is typically 10“"^ photons/electron and is as high as 10'^ for an Ag 
tip on an Ag surface [14]. 

These TIP modes are similar to the collective charge density oscillations 
observed at two metal-dielectric interfaces in close proximity [30]. The lowest- 
order dipolar mode gives rise to a high probability of inelastic tunneling and 
photon emission. For an enhanced electromagnetic field under the tip, tunnel- 
ing electrons that have a wave vector parallel to the dipole mode can couple 
efficiently via the lET mechanism. A recent study on Ag/Si(lll) [31] sug- 
gested that higher-order multipole modes are responsible for higher-energy 
spectral peaks. If sample and tip are made of the same material, the reso- 
nance frequency of localized plasmons between tip and sample is given ap- 
proximately by [25,31] 
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where n is the order of the multipole mode (roughly speaking, n = 1,2,3 
correspond to the dipole, quadrupole, and hexapole mode, respectively), Up 
the bulk plasma frequency of the free electron gas, qn the wave number of 
the resonance multipole mode [qn = nqi ^ n(2Rd)~^/^], d the tip-sample 
distance, and R the tip radius. Using the lET-TIP mechanism, Berndt et 
al. have further demonstrated that photon spectra can be modeled fairly 
well for Cu, Au, and (to a lesser extent) Ag surfaces regarding the positions 
of maxima, cutoff wavelength, and signal intensities [32]. The intensity of 
plasmon excitation was also evaluated by Tsukada et al. [33] using the first- 
principles calculation within the framework of surface plasmon excitation via 
inelastic tunneling. 



2.1.4 Factors Affecting STM-induced Photon Emission 
and Modes of Measurements 

STM-induced photon emission is determined by a variety of factors. In addi- 
tion to the tunneling current and bias voltage, the lET-TIP model implies 
that the tip-sample distance, tip radius, and work functions and dielectric 
properties of the tip and sample all affect the local TIP modes. In order to 
separate the factors that affect photon emission, STM-induced light emission 
experiments are mainly carried out in three different modes: photon mapping, 
optical emission spectroscopy, and isochromat spectroscopy [4]. The lateral 
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extent of these lET-TIP modes on fiat surfaces is theoretically determined 
to be below ~ 10 nm [34] (the localization length of localized plasmons is 
equal to [35]), and hence we can expect contrast in photon maps for 

surface features with geometric dimensions on that scale. On the other hand, 
electronic properties determine the initial and final density of states avail- 
able for inelastic processes, and they will modify the branching ratio between 
elastic and inelastic channels and hence photon emission. Consequently, one 
expects local variations in such factors to give rise to contrast in maps of the 
integral photon intensity (i.e., photon maps). Analysis of photon maps has to 
take into account both geometric and electronic aspects via bias-dependent 
topography and spectroscopy. Photon maps [36] and emission spectra [42] 
with atomic spatial resolution were recently reported for the Au(llO) surface 
with different features for over-row and between-row atoms. Polarization ef- 
fects of photon emission have also been investigated for Co films and were 
found to be related to magneto-optical effects [38,39]. 

2.1.5 The Debate on Mechanisms 

Although the theory of the tip-induced plasmon mode via inelastic tunneling 
(lET-TIP) is generally accepted for the modeling of STM-induced photon 
emission on metal surfaces, there is still some debate about the retardation 
effect and even the role of localized surface plasmons. A recent fully-retarded 
calculation by Johansson [40] indicates that the inclusion of retardation ef- 
fects does not change the basic picture of STM-induced photon emission 
described above but does change spectral positions at a quantitative level. 
In general, retardation effects become important when the sample or tip is a 
good conductor (e.g., Ag). The energy dissipation is then rather small in the 
near-field zone and a localized mode will instead be damped because some 
energy propagates away from the near-field zone. Retardation therefore leads 
to an additional redshift of the resonance frequency. For Au and Cu samples, 
the retardation effect is negligible at the resonance frequency 2.1 eV for 
Au and 2.0 eV for Cu) and the emission intensity is essentially unchanged. 

Other theories have also been proposed to explain STM-induced photon 
emission. Since the tunneling distance, the localization length of localized 
plasmons, and the characteristic plasmon length of the electron gas in the tip 
and sample all lie in the nanometer (or subnanometer) range and are quite 
comparable to each other, a precise theoretical model of inelastic electron 
tunneling should take into account the spatial dispersion of the dielectric 
function of the electron gas, i.e., the non-local effects. Ushioda and coworkers 
[41,37] developed a theory for visible light emission from STM on the basis of 
the theory designed for solid-state tunnel junctions. Their theory concerns the 
differences between the branching ratios of surface plasmon polaritons that 
propagate along the surface (in other words, are delocalized) and those that 
are localized under the tip (TIP modes) and emit light directly. A radiation 
mechanism through the Cherenkov effect (or ‘inverse Landau damping’) was 
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also proposed to explain the resonant excitation of interface plasmons by 
tunneling electrons via phase matching between the plasmon phase velocity 
and the Fermi velocity of electrons [43]. 

Recently, Xiao opened a debate about the importance of surface plasmons 
in the process of photon emission [44-47]. Specifically, does light emission 
require the excitation and de-excitation of surface plasmons? Using the point 
dipole model, he proposed the following simplified Hamiltonian to explain 
light emission from the tip-sample junction: 

H — Hq + eV cosujt , (2.2) 

where Hq is the Hamiltonian of the electron moving through the potential 
barrier between tip and sample, V is the bias voltage applied, and uj is the 
characteristic frequency of the tunneling current. According to Xiao, a system 
with such a simplified Hamiltonian is able to demonstrate that the coupling 
modes inside the tip-sample cavity can be excited by the fluctuation of tun- 
neling currents in analogy with the situation in cavity quantum electrody- 
namics [48]. These excited coupling modes are responsible for light emission 
in STM; the use of local surface plasmons is not needed to explain the ob- 
served photon emission. It is worth noting that the interaction term in (2.2) 
is constant in space and is not very suitable for inducing transitions due to 
the fluctuation of electric fields [45]. Nevertheless, these arguments indicate 
that the theory of photon emission induced by tunneling electrons is not yet 
well established. 

2.1.6 Molecules on Surfaces 

In contrast with the extensive investigations of STM-induced light emission 
on semiconductor and metal surfaces, only a few scattered reports have been 
found on photon emission from molecules adsorbed on metal surfaces. The 
scarcity of data on molecules is mainly due to the difficulty in obtaining 
reliable photon maps and optical spectra without modifying or damaging the 
molecules. 

The first (and probably by far the best defined) photon map with molec- 
ular resolution was obtained by Berndt et al. [49] with a simultaneous topo- 
graph [Vs — 2.8 V, It = 4.4 nA) for a monolayer of Ceo on Au(llO). The Ceo 
molecules appear as a hexagonal array of individual photon emitters with 
diameter ~ 4A and separation ~10A. The high symmetry and stability of 
Ceo molecules are the key to the success of this experiment. The authors 
observed more intense emission when the STM tip was positioned above a 
Ceo molecule and went on to claim that the molecule itself was the primary 
photon source. However, due to the lack of optical spectra, it is still not con- 
clusively demonstrated that the observed molecular contrast on photon maps 
involves molecular fluorescence from Ceo- The Ceo molecular layers could act 
like some kind of grating which modifies the spatial distribution of light 
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emission from the Au(llO) surface. In any case, it is clear that the molecules 
couple strongly with the electromagnetic modes of the cavity between the tip 
and sample and thus play an active role in the emission process. 

Excited molecules are known to couple with the plasmons of a metal 
surface [50]. When a molecule is sandwiched between tip and substrate, 
molecular fluorescence may be enhanced via a similar mechanism to that 
in surface-enhanced Raman spectroscopy (SERS) [51]. Electromagnetic con- 
flnement may mediate and amplify emission from the coupled tip-molecule- 
substrate system. To see how emission is modifled by such conflnement, it 
would be desirable to obtain the optical spectra for such a configuration and 
compare them with photoluminescence and electroluminescence studies. 

Berndt et al. also reported molecularly resolved photon maps of individ- 
ual anthracene molecules on Ag(llO) [4] {Vs = —2.6V, A = 30 nA), but 
there was no simultaneously acquired topography image to compare with 
and the characteristics of molecules after photon emission remain unidenti- 
fied. According to Berndt, the molecular emission is more intense than that 
from the Ag substrate. The emission intensity from anthracene molecules is 
strongly dependent on bias polarity whereas emission from the Ag substrate 
is relatively insensitive to it. 

The nonradiative damping of molecular excited states by metal surfaces 
is thought to drastically quench fluorescence when molecules are in close 
proximity to metallic surfaces [52]. To suppress the quenching of molecu- 
lar emission on metal surfaces, Flaxer et al. [53] used the transparent and 
conducting thin film of indium-tin-oxides (ITO) as the substrate for the de- 
position of submonolayers of dichloroanthracene or coumarin molecules. No 
STM images, photon maps, and optical spectra were reported there, so it 
remained unclear what the surface looked like and what happened to the 
molecules under a tunneling current up to 50 nA and a bias voltage above 
2.6 V. Nevertheless, the bias and coverage dependence of photon yields from 
the squeezable tunnel junction does strongly suggest the presence of light 
emission associated with molecular fluorescence excited by inelastic electron 
tunneling, although the 10“^% quantum efficiency estimated is still very low. 
The poor emission efficiency from molecules on metal and semiconductor sur- 
faces is due to the fast deexcitation processes governed by energy transfer to 
nonradiative surface excitations (Auger processes) [54]. 

Another STM-induced light emission from molecules in air was carried out 
for monolayers of copper phthalocyanine (CuPc) molecules on Au(lll) [43]. 
Pleated structures of CuPc molecules were resolved on both STM topographs 
and photon maps but with a poor resolution, partly due to the large current 
used (lOnA). On the basis of photon maps and bias dependence of emission 
intensities, the authors viewed each CuPc molecule as a source of secondary 
plasmons upon interaction with the interface plasmons excited by tunneling 
electrons. The periodicity of the molecular structure should be preserved in 
the photon map, but the positions of emission maxima do not necessarily 
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correspond to the positions of individual molecules. It is worth noting that 
measurements in air may suffer from complications caused by contamination 
and the resultant altered tunneling behavior as well as electrical breakdown 
and surface modifications at high bias voltages [4,55]. 

Recently, STM-induced photon emission has been reported for submono- 
layers of Cu-TBP porphyrin molecules on Cu(lOO) in ultrahigh vacuum via 
near-field detection from a conductive fiber tip [56]. However, it is difficult 
there to discern how the molecules are related to the correlation between the 
photon map and simultaneous topograph acquired under tunneling conditions 
of 6.0 V and 5.0 nA, since the molecules are likely to be damaged in such a 
high-current, field-emission regime. The authors estimated the quantum ef- 
ficiency to be about 2 x 10~^ photons/electron and claimed that Cu-TBP 
porphyrin molecules contribute to light emission on the basis of the enhance- 
ment with respect to the noble metal surface and a faint peak at 2.2 eV in 
optical spectra. Light emission with a higher quantum efficiency appears quite 
possible if a designed molecular chromophore such as hexa-butyl decacyclene 
(HB-DC) is used [57]. 

As discussed above, the theory of photon emission induced by tunneling 
electrons on metal surfaces is still not yet well established. The situation 
becomes even more complicated and more difficult to study when molecules 
are sandwiched between the tip and surface. The acquisition, distinction, and 
separation of molecular fluorescence with respect to photon emission from the 
tunnel junction or metal surfaces is not the only experimental challenge. In 
addition, the mechanism whereby electrical energy is transformed into light 
energy via molecules remains a puzzling issue. For tip-molecule-substrate 
structures, the details of this electroluminescent process - energy absorption 
or excitation, energy transfer, and emission - are still far from clarified. The 
focus is on the coupling of molecular electronic states with the local elec- 
tromagnetic field, in particular with regard to the precise way in which a 
molecule couples with tunneling electrons, gets excited, and gives out light. 

We have recently initiated an investigation of organic molecules such as 
porphyrins bound onto a metal surface. This work is conducted with the 
following objectives: 

• to understand the mechanism and underlying physics involved in the 
emission process, in particular the electron transport through a single 
molecule and electron-photon conversion, 

• to explore the ultimate spatial and spectral resolution of photon emission 
for chemical mapping and molecular recognition on surfaces, 

• to provide guidelines for improving the quantum efficiency of electrolu- 
minescence (EL) via molecules for future molecular electronics. 



We report here our preliminary results on STM-induced photon emission 
from Cu(lOO) covered by porphyrin molecules, and present a wavelength- 
resolved optical spectrum that appears consistent with molecular vibrational 
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transitions. Some general issues on experiments, molecules, and tunneling are 
also discussed throughout the paper. 

2.2 STM-Induced Light Emission 
from Porphyrin Molecules 

2.2.1 Experimental Setup 

It is no easy task to perform a successful STM-induced light emission ex- 
periment. It requires a well-defined sample, good tip, efficient photon collec- 
tion system, and highly sensitive photon detection system. Photomultipliers 
(PMT), grating spectrometers, intensified diode array detectors, or intensified 
CCDs are usually used to detect and spectrally resolve the photon emission. 
Since the photon intensity is very low due to the low quantum efficiency 
from the STM tunnel junction, a large collected solid angle becomes one of 
the most desirable features in designing a photon collection system. These 
systems can be divided into three main categories: mirror, lens, and fiber 
setups. The ellipsoidal or parabolic mirror system focuses light from the tip 
region onto the focal point outside the UHV close to a viewport [31,58]. This 
setup covers a large solid angle up to 61% of 27 t, but is relatively difficult to 
make and align. The fiber bunch system is relatively easy to make and can 
also offer a large collection efficiency but it is difficult to make alignments 
[19]. The optical fiber tip system can be classified in this category with its 
unique feature of collecting photons at the near field [17]. Lens systems are 
convenient to use and are the most widely used setup nowadays in STM- 
induced photon emission studies [9,18,42]. However, they suffer from a small 
solid angle of collection. This is also the setup we use for photon collection 
(Fig. 2.3), the so-called far-field detection, since the distance between lens 
and light source is much greater than the wavelength of light emitted. 
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Fig. 2.3. Schematic diagram of our photon detection system through a double lens 
setup in combination with UHV-STM 
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Experiments were performed in an ultrahigh vacuum system (UHV) with 
facilities for sample heating and cooling, molecule deposition, and characteri- 
zation by a scanning tunneling microscope (JEOL JSTM-4500XT). The base 
pressure of the system is less than 2 x 10“^^ torr. Experiments were conducted 
at 80 and 295 K. Cu(lOO) substrates 3 x 7 x 0.5 mm^ were mechanically pol- 
ished using 0.3 pm alumina paste followed by electrochemical polishing in a 
solution of 66% H3PO4 and 34% methanol (20 mA, 2.8 V, 10 min, 4°C). These 
were thoroughly rinsed in methanol before being loaded into the UHV cham- 
ber. Atomically clean Cu(lOO) surfaces were obtained after cycles of sputter- 
ing (Ar~^, IkeV, ^ IpA, ~ 800°C, 30 min) and annealing (500°C, 30 min). 
Cu-tetra-(3,5 di-tertiary-butyl-phenyl)-porphyrin (Cu-TBPP) molecules were 
deposited onto Cu(lOO) by sublimation from a Knudsen cell at ^ 300° C 
or higher. The sample was then either used as is or annealed at ~ 200° C 
for 5 min. All images were taken in the constant-current topographic (CCT) 
mode with the sample biased. Tips used for imaging and photon emission 
were prepared by electrochemical etching of a W wire in a 2N NaOH solution 
followed by in situ cleaning in vacuum by Ar+ ion sputtering or heating the 
tip apex up to ~ 1 000° C. 

A schematic view of our light collection and detection system is shown 
in Eig. 2.3. A condenser lens, 15 mm in diameter, is placed in UHV close to 
the tip-sample region (27.4 mm away) at an angle of 45° with respect to 
the surface. The collected solid angle is 0.22 sterad. Photons emitted from 
the tunneling gap are collected by this lens, transmitted through a view- 
port, and then refocused onto an optical fiber using another condenser lens 
located outside UHV. The fiber cable is then connected to either a PMT 
or a spectrophotometer (Hamamatsu PMA-100). The spectrometer includes 
a grating monochromator (C5094) and an intensified charge-coupled device 
(ICCD) with Peltier-cooled multialkali (S-20) and multichannel plate (MCP). 
The ICCD operates in a pulse counting mode at a dark count rate of ca. 20 
counts per second (cps) at — 15°C for the wavelength range 350-850 nm. To 
ensure that emitted photons are counted properly, the intensity was also mea- 
sured by a photomultiplier (Hamamatsu H6180-01, dark counts ^ 10 cps at 
25°C, 300-650 nm), and a cooled digital CCD camera (Hamamatsu C4880-40, 
average dark counts 0.02 electron/pixel/sec at -50°C, 300-1100 nm). 

2.2.2 Porphyrin Molecules 

The electronic characteristic of metal-free porphyrin molecules is the 18 
electrons delocalized along the inner 16-membered tetrapyrole ring of the 
molecule, forming a closed, completely conjugated aromatic system. This elec- 
tronic structure is believed to be responsible for its optical behavior. Indeed, 
the optical absorption spectrum of porphyrins is successfully rationalized by 
Gouterman’s four-orbital model [59]. This simple model gives rise to a rel- 
atively weak (pseudoparity-forbidden) Q-band in the visible region and an 
intense (optically allowed) Soret or B-band in the near-UV. These features 
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are nicely reproduced in the optical absorption spectra of H 2 -TBP porphyrin 
molecules with a fine- structured Q-band around 519.3 nm and a sharp B- 
band at 422.8 nm. The free-base TBP porphyrin molecules fluoresce strongly 
at 649 nm. Similarly for Cu-TBP porphyrins, absorption spectra show a small 
Q-band at 539.6 nm (2.3 eV) and a sharp B-band at 417.0 nm (3.0 eV), but 
they do not fluoresce at all. The ability of porphyrins to achieve various states 
of oxidation and the ease with which they bond to various metals has made 
them a major functional constituent of biological and molecular electronic 
systems. The introduction of metal atoms into the center of molecules enables 
us to tune not only oxidation states but also photophysical and photochem- 
ical properties of the molecules. A critical issue is how to understand their 
functions on the molecular level. To be more specific, what are the kinetic 
behavior and the mechanisms by which porphyrins oxidize, transfer energy, 
and relax? What is the effect of the central metal atom on the relaxation 
mechanism of electronically excited porphyrins [60]? 

It is generally observed that metalloporphyrins with closed shell metal 
ions of Mg, Ca, and Zn exhibit fluorescence and phosphorescence, while open 
shell diamagnetic metal complexes of Ni do not emit at all; Pd complexes 
emit weak fluorescence and strong phosphorescence, while Pt complexes do 
not fluoresce, but do phosphoresce strongly. Paramagnetic metal complexes 
such as Cu, Ag, and Au either do not show any emission at all, or, as in 
copper complexes, exhibit short-lived phosphorescence [61]. Fe(II)-centered 
porphyrins relax to the ground state via nonradiative decay within 6psec 
after excitation. In contrast, the metal- free porphyrin fluoresces strongly and 
its relaxation time to the ground state is longer than 1 nsec. Evidently, the 
central Fe(II) metal plays a predominant role in the electronic energy excita- 
tion of these porphyrins which proceeds via radiative or nonradiative decay 
channels. 

Specifically for Cu(II)-porphyrins (oxidation states can be verified by 
ESR) [60], the interaction between the porphyrin electron system and the 
unpaired d electron of Cu(II) (d^, 25 H- 1 = 2) means that the ground state 
becomes a doublet ^5 q (25 -fl = 2). The lowest excited singlet state becomes 
the lowest singdoublet ^5i that is closely related to the lowest excited singlet 
state of porphyrin. The lowest excited triplet states are tripdoublet and 
quartet ^Ti (25 + 1 = 4) originating from the lowest excited state of por- 
phyrin. The singlet states of the metal-free porphyrin couple with the doublet 
levels of the central metal atom to form ‘singdoublet’ states (^5i). The triplet 
states couple with the doublets to form ‘tripdoublets’ (^Ti) and quartet states 
C^Ti). The coupling of singlet, doublet, and triplet states in Cu-centered TBP 
porphyrin molecules is responsible for the absence of fluorescence. A detailed 
kinetic investigation will help to clarify the decay pathways and the nature 
of transitions between different states. 

In general, a good photoluminescence material is not necessarily a good 
electroluminescence material; furthermore, a bad photoluminescence material 
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is usually not good for electroluminescence unless triplet states are used for 
photon emission. 

2.2.3 Tunneling Transport Through Molecules and STM Imaging 

An STM is a typical two-terminal device. Tunneling between two electrodes 
consists of elastic and inelastic portions with the former predominant. The 
inelastic tunneling channels contribute only a small fraction (less than 0.1%) 
to the total tunneling current and are affected by factors such as geometry, 
local dielectric properties, and density of states of both tip and surface [6]. 
The molecular image observed by STM is a combined result of topography 
and electronic structures of the molecule, tip, and substrate. To understand 
STM imaging from elastic tunneling and photon emission induced by inelastic 
processes, it is important to have a general picture of the transport behavior 
through individual molecules between two electrodes [3,62]. 

At low voltages, two regimes of elastic transport explain the large scat- 
ter of transparency values measured on single molecules. For a transmission 
coefficient larger than one (T > 1), the transport regime is ballistic: the 
molecular levels are in resonance with the Fermi level of the electrodes, and 
there is complete electronic delocalization along the molecule. For example, 
in the case of single-wall carbon nanotubes at low temperatures, the trans- 
port coherence length can extend to as much as 500 nm. Values of T < 1 
indicate that there are some scattering events during electron tunneling: ei- 
ther the molecule is differently bound to the electrodes or the Fermi level is 
located within the gap between the highest occupied (HOMO) and the low- 
est unoccupied (LUMO) molecular orbitals. In the first case, the tunneling 
regime is determined by the metal-molecule contact, which can be optimized 
by controlling the surface chemistry. To overcome the second situation, the 
molecular orbitals need to be electronically coupled to the electrodes to sta- 
bilize a tunnel pathway for the carriers that is more efficient than vacuum. 
This pathway builds up from the constructive or destructive superposition 
of tunnel channels. For systems with T < 1 and involving a molecule sym- 
metrically chemisorbed on the electrodes, such as Ceo bound to two similar 
electrodes [63], two-terminal measurements exhibit linear I-V characteris- 
tics when V is much lower than the molecule’s effective barrier height. This 
Simmons tunneling regime is analogous to electron tunneling between metals 
at low voltages, resulting in linear I-V curves. Nonlinear I-V curves usu- 
ally occur at higher voltages in systems with T > 1 involving a physisorbed 
molecule, owing to the dominant effect from the contact tunnel barriers at 
the electrodes. 

Tunneling through atoms and molecules whose filled and empty electronic 
states do not lie around the Fermi level occurs by a process known as vir- 
tual or off-resonance tunneling [3,64,65] at low voltages. The electrical re- 
sistance of the junction at contact is defined by the electronic transparency 
of the molecnle. Transparency is a measure of the molecule’s efficiency in 
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extending the metallic wave functions of the electrodes. It is approximately 
proportional to the square of the inter-electrode electronic coupling intro- 
duced by the molecule with respect to the corresponding vacuum gap. Ohmic 
dissipation in the electrodes is one way to evaluate electronic transparency 
(T) from the macroscopically measurable quantity of tunnel resistance i?, as 
T = (/i/2e^)i?~^, where e is the elementary charge and h is Planck’s con- 
stant. However, owing to R being the resistance of the metal-molecule-metal 
tunnel junction, rather than that of the molecule itself, it cannot be used to 
define molecular conductivity. 

The coupling of molecular electronic states with those of electrodes mainly 
occurs through the tails of the molecular orbitals, in particular those of the 
highest occupied and lowest unoccupied orbitals. These frontier orbitals pro- 
vide such an extension through the induced shift of molecular orbital levels 
[66] and are sensitive to changes in hybridization and degeneracy induced 
by external distortion. For example, a mechanically induced distortion of 
0.1 nm is sufficient to modulate the current flow through Ceo by a factor of 
100 [63]. Increasing this distortion leads to unit transparency, equivalent to 
a quantum unit of conductance 2e^//i. Molecular extension of the metallic 
wave functions is particularly relevant to devices operating at low voltages 
and where linear current-voltage characteristics are desirable. On the other 
hand, single-electron tunneling associated with Coulomb blockade effects has 
been observed for Ceo molecules when they are deposited on an insulating 
layer at 4.2 K [67]. The wave function decoupling of molecules is suggestive 
of a single-electron device using a molecule as an ideal quantum dot [68] . 

The above description of transport through individual molecules suggests 
that access to the true molecular conductivity is a challenging issue. But by 
knowing a variety of factors involved in the transport process for a molecule 
sandwiched between two electrodes, one can have a better understanding of 
STM imaging, coupling of molecules with electrodes, molecular conformation 
on surfaces, and even the molecular manipulation process [69,70]. 

Figure 2.4 shows the structure and conformation of a Cu-TBPP molecule. 
The steric repulsion among the bulky substituents (legs) on the phenyl drives 
the phenyl rings roughly perpendicular to the porphyrin ring. STM images 
of a Cu-TBPP molecule on Cu(lOO) consist of four bright lobes, as shown 
in Fig. 2.5a for a very low coverage and Fig. 2.5b for a coverage close to one 
monolayer. The distance between the trans-lobes is ~ 1.8 nm, slightly less 
than the molecular dimension of ~ 2 nm. The bright lobes are thus assigned to 
the di-t-butylphenyl (DTP) side-groups, in accordance with previous reports 
[69,70] and theoretical simulation [71]. The characteristic four-lobe pattern 
is a registry of Cu-TBPP and is used for subsequent molecular recognition 
and controlled tip-positioning above a molecule. 

The molecules are found to be relatively mobile on the surface at room 
temperature. As regards the stability of molecules against bias voltages and 
currents, our STM studies show that the molecules are stable below 3.8 V 
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for a small tunneling current (e.g., 0.1 nA), but are easily broken into frag- 
ments when bias voltages are above 4.5 V. The application of a large current 
tends to drag the molecules across the surface during scanning. Porphyrin 
molecules are chosen because they are fluorescent molecules with aromaticity 
and are known to have high photoluminescence efficiency for converting tt-tt* 
electronic excitations into visible and near-UV fight in their free state or in 
solution [72]. However, in close proximity to metallic surfaces, nonradiative 
damping of molecular excited states by metal surfaces is thought to dras- 
tically quench fluorescence [52]. Our experiments on planar TPP-porphyrin 
molecules [with the C(CH 3 ) replaced by H] appear to support the above 
claim. In order to suppress the fast electron energy dissipation to metal sur- 
faces, a porphyrin architecture that is electronically decoupled from the sub- 
strate appears necessary. The use of bulky t-butyl groups makes the plane 
of phenyl rings rotate out of the plane of the porphyrin ring, and this dra- 
matically weakens the tt electron delocalization over the porphyrin-phenyl 




Fig. 2.4. (a) Structure and (b) conformation of Cu-TBP porphyrin molecules 
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(a) 










Fig. 2.5. STM images of Cu-TBPP on Cu(lOO) at 295 K: (a) Small coverage, 9 x 
9nm^, +0.5 V, 0.2 nA, no annealing applied; (b) monolayer coverage, 25 x 25 nm^, 
+1.5 V. 2nA, sample annealed at ~ 200°C for 10 min. The inset is the top view 
of the electronic density distribution of the molecule simulated by a semi-empirical 
INDO calculation. The features observed are related to the four TBP legs 



bond. The electronic states of the porphyrin core are thus electronically de- 
coupled from the surface since it hangs ~ 5 A above it. Intense STM-induced 
photon emission is indeed observed from this molecule (see below). In ad- 
dition, the crown-like cap configuration of Cu-TBPP together with its weak 
chemisorbed character on Cu(lOO) might serve as a squeezable tunneling 
junction and might allow the molecule to sustain a relatively large current 
for photon emission. 

2.2.4 STM-Induced Photon Emission 

Photon emission induced by tunneling electrons for the present system was 
initially detected by a CCD detector when the tip was positioned above 
molecules and a moderate bias voltage (e.g., 4 V) was applied. The emission 
event can occur for both ‘on’ and ‘off’ states of the feedback loop. The ‘on’ 
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state gives reproducible and long-lasting photon emission (hours, but ‘blink- 
ing’) while ‘off’ state light emission lasts for only a few seconds, presumably 
due to the lack of a feedback loop to maintain a stable tunneling gap. Our 
intensity measurements indicate that light emission starts when the bias volt- 
age is above 1.8 V in both polarities, and appears to reach a local maximum 
around 3.5-4 V. Further increase of the bias voltage up to 7 V (field-emission 
regime) often leads to short but incandescent emission (‘bleaching’), which 
suggests that the molecule is modified or damaged, as revealed by subsequent 
STM imaging. On the other hand, the emission intensity increases almost 
linearly with current, as expected, and becomes visible to the naked eye at 
a current of around 50 nA. For a surface completely covered by Cu-TBPP 
molecules, since the tip is almost always above a molecule, light emission can 
be realized either when the tip is kept stationary above a molecule or during 
scanning. Variation of the bias from 2 to 5 V did not change the color of the 
emitted light, which appeared orange-red to the naked eye when the environ- 
ment was relatively dark. It is worth pointing out that light is emitted only 
when the current fluctuates to some extent, indicating a relationship between 
photon emission and tunneling current behavior, although a quantitative re- 
lation is still hard to discern. 

The first indication of enhanced light emission due to the molecules comes 
from our experiments performed on the surface in Fig. 2.5a with a very low 
coverage of Cu-TBPP molecules. Intense photon emission was readily de- 
tected by a CCD when the tip was positioned above a molecule, but little 
light was detected by CCD when the tip was located above the bare Cu 
surface under the same excitation conditions (e.g., 3.5 V, 10 nA); the latter 
was consistent with experiments performed on a clean Cu(lOO) surface. This 
observation agrees with reports concerning Ceo on Au(llO) [49] in which the 
detected photon intensity is greatest when the tip is placed above an indi- 
vidual Ceo molecule. A molecular origin for light emission was also claimed 
previously by Flaxer et al. [53] in their study of 9,10-dichloroanthracene on 
an ITO substrate. A similar argument is proposed here based on the above 
observation and the optical spectral data given below. Since an STM tip is 
very sharp and generally believed to have only one atom at the tip apex, 
the enhanced emission above a single molecule also suggests that even single 
molecules can have a significant influence on the photon emission process. 



2.2.5 Intensity and Quantum Efficiency 

We are still pursuing efforts to obtain well-defined photon maps for molecules 
on surfaces at low tunneling currents through improvements in collection 
and detection systems. In this chapter, we report our preliminary results on 
emission intensity and optical spectra. Photon intensity and optical spectral 
data were acquired from a monolayer-covered surface in Fig. 2.5b instead of 
the one in Fig. 2.5a. This is because spectral data with a decent S/N ratio 
require a long integration time, but a low-coverage surface makes it difficult 
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to hold the tip above the same molecule for a sufficient time due to the 
thermal drift problem. Although the intensity of photons emitted from an 
STM is generally low, the signal-to-noise ratio is actually quite high, above 
50 in the present setup. Count rates of up to 10^ cps per nA over 350-850 nm 
were recorded even for the small detection solid angle of 0.22 sterad. The 
total intensity was of the order of 10^ cps for a tunneling current of around 
10 nA after calibrating the quantum efficiency of the multialkali plate, MCP, 
solid angle, and attenuation or loss through transmission. (One can rule out 
the hypothesis that the variation in photon intensity arises from variations 
in tunneling currents by establishing that no direct structure correlation is 
observed in simultaneously recorded maps of the tunneling current.) The 
quantum efficiency of photon emission induced by STM for the present system 
was thus estimated to be around photons/electron. This value 

is about ten times larger than photon emission from pristine noble metal 
surfaces [6,28,58]. The enhancement of photon emission due to molecules is 
evident. 

The so-called tip-induced plasmon (TIP) mechanism [20] alone, i.e., pho- 
ton emission through inelastic tunneling processes with radiative decay of lo- 
calized dipolar surface plasmons, cannot explain the enhanced emission above 
molecules. The presence of a molecule would lead to an increase in the dis- 
tance between the tip and metal surface and thus to a smaller field enhance- 
ment (or TIP). This would instead result in intensity minima above molecules, 
contradicting the experimental observations. Conversely, the molecule plays 
an active role in the photon emission process and may even be an important 
photon source itself. Enhanced photon emission above molecules implies a 
strong coupling between the molecule and the electromagnetic modes of the 
cavity between the tip and metal surface [50]. In addition, the fiuorescence 
of molecules may be enhanced in a way that is similar to that invoked in the 
classical explanation of surface-enhanced Raman scattering [51]. Another fac- 
tor for the enhancement may pertain to the quantum confinement of electrons 
in the molecule [73] and the associated resonant tunneling process. 

In principle, a direct electron-photon conversion is a forbidden process in 
the absence of a potential field because of their diflFerent dispersion relations. 
Surface plasmons on a perfectly smooth metal surface cannot radiate since 
such a process could not simultaneously conserve energy and momentum. 
However, when a tip is brought near the surface or when the surface is rough- 
ened, translational invariance along the surface is broken and photon emission 
can occur upon applying appropriate bias voltages. Although it has not yet 
been well understood how a molecule couples with the local electromagnetic 
field and enhances photon emission, it is reasonable to consider the molecule 
as a quantum system where tunneling electrons are more or less confined in 
the molecule during their transport. The confinement-induced enhancement 
in the radiative rate and reduction in the nonradiative rate would result in a 
higher efficiency of luminescence from the semiconducting molecule [73]. In a 
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sense, each molecule on the surface can be viewed as a protrusion, which acts 
as a scattering center facilitating the conversion of plasmons into photons 
[74]. The unique crown-like conformation with bulky legs also helps to sup- 
press the otherwise fast nonradiative decay [54] by weakening the interaction 
of the molecule with the substrate. 

The model developed by Johansson et al. [29] takes the role of the tip into 
account by approximating the tunneling from a metallic sphere to a planar 
metal surface. The intensity of emitted light from tunneling electrons was 
found to split into two separate factors. One is the strength of the electro- 
magnetic field in the vicinity of the tip giving us a measure of the coupling 
strength between tunneling electrons and electromagnetic field. The other 
is a fluctuation in the tunneling current that acts as a source for radiation. 
The presence of a tip close to a sample surface creates a localized interface 
plasmon built up from surface charges of opposite polarity on the tip and 
sample surface, respectively. This plasmon resonance leads to a considerable 
enhancement of the light emission from tunneling electrons. 

Similar arguments involving current fluctuations may explain the en- 
hancement effect for a surface with absorbed molecules. An electronically 
excited molecule can be viewed as an oscillating charge distribution or a dy- 
namic dipole moment surrounded by an oscillating electric held. Molecules 
are very likely subject to structural deformation during tunneling excitation 
at high bias and current. This could cause a small change in the molecular or 
tunneling resistance and thus result in small fluctuations of tunneling currents 
even though a feedback circuit is used to maintain a constant current. 

Two models have been proposed to address the excitation mechanism 
of molecules [57]. The first model views the tip-molecule-substrate config- 
uration as a double barrier structure. When the molecular electronic states 
are aligned with the energy bands of the electrodes, molecules are directly 
excited by inelastic tunneling. This is similar to the ballistic regime men- 
tioned above for transport through molecules [3,62]. The second model, a 
more favored one, considers the tip-molecule-substrate configuration as a 
single barrier structure. Inelastic tunneling excites localized electromagnetic 
modes between the tip and substrate, and then these lET-TIP modes excite 
the molecule indirectly via coupling with molecular electronic states. Since 
the photon emission is related to the dielectric functions of both tip and sam- 
ple, the dielectric properties of the molecules inside the tunneling gap will 
probably affect both the lET-TIP modes and lET probability. The other 
factor for the enhancement is related to the molecular luminescence, which 
contributes directly to light emission (see below). 

2.2.6 Optical Spectra 

Further evidence for the enhancement effect and even fluorescence contri- 
butions to light emission from molecules is revealed in the optical spectra 
acquired over a monolayer of porphyrins on Cu(lOO). Figure 2.6 shows the 




42 



Z.-C. Dong et al. 



14,000 



I'lJ.SM 

s 

S' 

« 

V 13,000 



12,500 

4S0 500 SSO eoo $50 700 750 SOO 050 

Wav«l«ngth (nm) 

Fig. 2.6. Optical spectra induced by STM on Cu(lOO) at 295 K for pristine Cu(lOO) 
and Cu-TBPP molecules 




typical optical spectra of Cu-TBP porphyrins induced by STM compared 
with those on the bare Cu(lOO) surface. The main peak occurs at ~ 610 nm 
for a surface with porphyrin molecules, very different from the broad peak 
around 700 nm for the noble metal surface. Since the photon emission is re- 
lated to the geometry of both tip and sample, peak positions may shift a bit in 
different experiments. A shifted optical spectrum is shown in Fig. 2.7, which 
was measured over 200 s from a different sample (monolayer coverage) and 
using a different tip. The excitation condition was set to +4 V and 100 nA. (It 
is worth noting that a large current was used in these experiments in order 
to get well-defined spectra, and this may complicate spectral analysis. The 
photon intensity in Figs. 2.6 and 2.7 is not corrected for the device response.) 

In general, for bias voltages less than 2 V, no meaningful peaks were ob- 
served on the spectra even for a long integration time, suggesting a threshold 
around 2V for light emission. The threshold at Vb ~ 2V is in agreement 
with the HOMO-LUMO gap of Cu-TBPP molecules {Eg ^ 2eV). Visible- 
light emission from molecules can occur only when Vb ^ no matter what 
kind of mechanism is involved, whether it be direct injection of electrons into 
the molecule or excitation of the molecule by coupling with local electromag- 
netic or plasmon modes. It is worth noting that adsorption of Cu-TBPP on 
Cu(lOO) results in a broadening of the HOMO and LUMO states, in partic- 
ular the LUMO [70]. Accordingly, the tip can also contribute to enhancing 
the tunneling current by pressure-induced distortion of the molecular orbitals 
and resultant broadening when in contact with the molecule. 

Figure 2.7 shows two broad peaks at 660 and 715 nm, respectively, both 
having contributions from the localized surface plasmons on Cu(lOO). Consis- 
tent with the spectra in Fig. 2.6, the peak at ~ 715 nm is almost invariant with 
respect to the excitation voltage and sample location, suggesting an intrinsic 
feature associated with the Cu surface. However, the broad peak from 550 to 
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Fig. 2.7. Optical Spectra of Cu-TBPP on Cu(lOO) at 295 K at a preset excitation 
condition of +4V and 100 nA integrated over 200 s 




680 nm presents new features that are not observed for the pristine surface in 
terms of intensity enhancement and peak shape. This peak is therefore partly 
attributed to the molecules. Since Cu-TBP porphyrins do not fluoresce, the 
observed high-energy broad peak is compared with the fluorescence peak of 
the free-base H 2 -TBP porphyrin molecules at 649 nm. The good agreement 
comes as no surprise because in any case fluorescence has the same origin as 
TT-TT* electronic excitations from the porphyrin cores. There is nevertheless 
a small possibility that the Cu-centered porphyrins are modifled to free-base 
porphyrins under the excitation condition for photon emission. The shape 
and intensity for this broad peak appear to be site-dependent, some with a 
shoulder-peak structure (Fig. 2.7), some with only a single peak (Fig. 2.6), 
implying that the spectra are related to the local molecular arrangements. 
The cutoff at 550 nm is attributed to the electronic interband transition in 
copper [20,75]. 

Another interesting feature in Fig. 2.7 is the presence of three shoulder 
peaks, at approximately 600, 627, and 662 nm, respectively, superimposed on 
the broad peak around 630 nm (^ 2eV). It is possible that the fine structure 
is caused by the different cavity modes between the tip and sample. How- 
ever, the peak spacing of the fine structure is ^ O.leV 800 cm“^, which 
falls in the range of molecular vibrational transitions. This suggests another 
possible cause for the fine structure: the photon emission may arise from the 
coupling of molecular electronic states with the local electromagnetic modes. 
Excited molecules then undergo vibrational relaxation of excited states be- 
fore decaying radiatively to the ground state. This speculation may be partly 
supported by the IR spectrum of the molecules, which does show a spectral 
peak around 800 cm“^. 

The attribution of this vibrational peak has not yet been clarified. How- 
ever, referred to the range of vibrational frequencies reported for the chemical 
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bonds of organic molecules, there are two possibilities. One is C-C stretching 
or bending between the phenyl-porphyrin single bonds (700-1250 cm~^), and 
the other is the C-H deformation vibration of ring hydrogens on the phenyl 
groups or the ring deformation vibration itself (700-900 cm~^). We speculate 
that the observed frequency at 800 cm“^ is related to one of the above vibra- 
tional modes, with the ring C-H deformation vibration being the most likely. 
A Franck-Condon-like (vertical transition) mechanism may be involved dur- 
ing the transition. The current or field fiuctuation is likely to cause molecular 
deformation followed by the vibrational relaxation of molecules. 

The observation of spectral peaks associated with the Cu-TBPP molecules 
is very exciting but puzzling, since Cu-centered porphyrin molecules are 
known not to fiuoresce under photoexcitation. Will the electroexcitation in- 
side a nanogap defined by STM bring out new physics about energy absorp- 
tion, transfer, or decay processes? Does the mixing or coupling of singlet and 
triplet states for molecules behave differently in localized electromagnetic cav- 
ities? Or is it possible that the Cu-TBPP molecules lose their central metal 
atoms and are changed to free-base molecules when a high electrical field 
and large current are applied across the molecular junction? Further research 
is required to clarify these issues. It is worth noting that the broad peak 
associated with the molecules appears dramatically enhanced when the bias 
voltage is above 3 V, in good agreement with the high absorption efficiency of 
the B-band at 3.0 V for the optical excitation of Cu-TBP porphyrins. These 
observations suggest that, although the excitation or decay mechanism of 
molecules by electrons might be different from that by photons, the energy 
absorption behavior appears similar on the basis of the energy dependence 
of optical spectra. 



2.3 Conclusion 

We have demonstrated a technique to produce tunneling-electron-induced 
photon emission from Cu-TBP porphyrin molecules on a Cu(lOO) surface. 
The emission intensity and optical spectral data show not only the enhance- 
ment effects of the molecule, but also new features revealing the molecular 
origin of light emission. The envelope of the spectral peak associated with 
molecular fluorescence appears to exhibit a fine structure that is consistent 
with the vibrational transition of the molecule. The quantum efficiency is typ- 
ically 10”^ photons per electron for Cu-TBP porphyrins. Since Cu-centered 
porphyrin molecules do not fiuoresce under photoluminescence, the observa- 
tion of unexpected light emission from the molecule may suggest new physics 
in the nanoscale STM excitation. 

Although the data from the present work and others are still not sufficient 
to reach a definite conclusion about the nature of the observed phenomena, 
the general physical process of STM-induced molecular light emission can be 
pictured as follows. When an STM is positioned above a molecule and a bias 
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voltage is applied across the tunneling gap, a net tunneling current flows from 
the tip to the substrate via the absorbed molecule. A tunneling electron or 
an absorbed molecule can excite a surface plasmon in the metal substrate, 
which causes an enhancement of the local electromagnetic held and in turn 
helps to excite the molecules. The electromagnetic excitation of the molecule 
occurs via the resonant coupling of molecular electronic states with the local- 
ized electromagnetic modes. The excited molecules relax vibrationally and 
then decay to the ground states with emission of light. The decoupling of 
luminescent porphyrin cores from the substrate, a means to suppress fast 
nonradiative energy dissipation, is found to be important for efficient light 
emission. 

On the other hand, the combination of STM with optical techniques 
shows promise for carrying out electronic and vibrational spectroscopy on 
the single-molecule level. This hybrid technique will facilitate the identifica- 
tion of individual molecules (‘molecular recognition’) and luminescent defects 
on surfaces, permitting chemical mapping on the molecular scale. The in- 
tense photon emission from single molecules may be a potential point source 
of light that could be used for high resolution imaging and even for flat- 
display technology. The idea of having a variety of molecular chromophores 
in a component to yield specific responses upon excitation may be a way to 
address the input and output issue of electronic signals via localized electro- 
magnetic modes. The use of plasmons from nanoclusters to amplify the gain 
in light transistors has recently been demonstrated for optical circuits [76] 
and similar concepts of plasmon enhancement may be useful for the design 
of molecular-scale amplifiers and transistors. In addition, an understanding 
of electron transport in the tip-molecule-substrate geometry may provide 
insights into some biological processes since electron transfer from electron- 
donor to electron-acceptor molecules via a molecular bridge is an important 
feature of many biological and chemical systems [77]. 
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3.1 Introduction 

Combined electrical and optical measurement methods applied to a sample 
are able to supply us with comprehensive information concerning condensed 
matter. However, up to now, such a combination has only been carried out for 
semiconductors. Optical methods are ineffective in metals, whereas methods 
based on conductivity are ineffective in organic systems, which are gener- 
ally insulators. The invention of tunneling microscopy [1,2] opened up new 
possibilities for combining electrical and optical methods because scanning 
tunneling microscopy (STM) enables us to use electric currents to excite sin- 
gle organic molecules. 

Up to now the main sources of information concerning physical charac- 
teristics of small junctions in STM are DC current-voltage curves for these 
junctions [3]. It should be stressed that STM provides the best space res- 
olution attainable by physical methods today. In fact, we can see a single 
atom using STM [4] . As far as the time resolution of scanning tunneling mi- 
croscopy is concerned, it is not particularly high and can only distinguish on 
the millisecond scale. STM does not ‘see’ any faster dynamics. Consequently, 
STM displays an averaged picture when faster relaxation processes exist in 
a junction. 

On the other hand, impressive progress has been made in single molecule 
spectroscopy (SMS) over the last decade [5-9]. SMS uses an excitation flu- 
orescence method that was widely employed before in the traditional spec- 
troscopy of molecular ensembles to measure absorption spectra. However, a 
laser source of light can provide space resolution comparable with the wave- 
length. Therefore, we have to use very low concentrations of guest molecules 
in SMS. When the distance between individual guest molecules becomes large 
we are able to excite a single molecule. By using near-held scanning optical 
microscopy (NSOM), we can raise the space resolution of the setup, although 
it remains inferior to that reached by STM [10]. 

However, SMS has some advantages over STM. Firstly, the time resolu- 
tion of SMS can reach the nanosecond level [11]. Secondly, SMS can probe 
guest molecules other than those lying on the surface. The third advantage 
stems from the possibility of using monochromatic laser excitation because 
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the frequency of laser excitation is an additional external parameter that can 
easily be changed in an experiment. 

The majority of the reviews in this book discuss problems of electric cur- 
rent in STM. We therefore decided to concentrate our discussion on problems 
concerning photons emitted in STM or registered in SMS with an emphasis 
on SMS methods. We intend to discuss in detail theory, experimental data 
and possible applications of SMS because this information may be useful in 
applying photon counting methods in STM and interpreting experimental 
data when STM probes organic molecules on surfaces. 

3.2 Light Emission from STM 

Let US start with a short review of a problem which concerns light emission 
from small tunnel junctions. A schematic view of STM is shown in Fig. 3.1. 

If a bias of a few volts is applied between the metallic tip and surface, 
an electric current from the tip to the surface emerges when the distance d 
between the tip and the surface falls to nanometric scales. The current is 
due to a single-electron transition through the potential barrier between the 
tip and the surface. The tunneling nature of the current is confirmed by its 
strong exponential dependence on the distance d. The tunneling current I is 
increased if the bias V is raised. 

The I{y) dependence shown in Fig. 3.2 is the main physical characteristic 
of the junction in STM. This is a nonlinear function. It reveals a breakdown 
of Ohm’s law because the junction conductivity depends on the bias. 

The bias-dependence of conductivity sometimes exhibits non-monotonic 
behavior, as can be seen in Fig. 3.2. The conductivity is proportional to the 
rate of single-electron tunneling through the potential barrier in the junction. 
As has been demonstrated by direct numerical calculation in a double-well 
potential model [13,14], the tunneling transition probability is increased by 
one order of magnitude if there is a resonance between levels in the left- and 
right-hand wells separated by a potential barrier. The oscillations observed 




Fig. 3.1. Tip-surface tunneling junction in STM [12] 
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Fig. 3.2. Conductance of the tunneling junction between W tip and Cu(lll) surface 
vs. bias [12] 



in this rate and shown in Fig. 3.2 can be ascribed to this resonant effect. 
This means that single-electron tunneling is accompanied by the creation of 
electromagnetic field quanta in accordance with the following energy balance 

eV = nTiQ , (3.1) 

where the left-hand side describes the change in the tunneling electron energy 
and n is the number of quanta created by the tunneling electron. 

Due to its high space resolution, STM is a very promising tool for probing 
single molecules on a surface. Various guest molecules have been investigated 
using tunneling current excitation, including Ceo [15-17], porphyrin deriva- 
tives [18,19] and molecular complexes [20]. Fluorescence excited by tunneling 
electrons was the main source of information concerning these molecules. 

However, there is a general problem in this type of spectroscopic exper- 
iment. The tip-surface tunneling junction can emit light even if there are 
no guest molecules in the junction. Transformation of electronic energy into 
light energy is detected in junctions which do not include any guest molecules 
[21,22]. How can we separate molecular fluorescence from light emitted by the 
junction itself? What is the mechanism for transforming electronic energy 
into light? 

There are various points of view concerning mechanisms for this energy 
transformation. From the point of view of classical electrodynamics, a con- 
stant bias creates a constant electric current. In this case only radiation of 
the Cherenkov type is possible [23]. This idea is discussed in [24]. However, it 
is unlikely that tunneling currents can be considered in this classical fashion. 

Johansson et al. [25] used formulae for quantum transitions in order to 
take the tunneling current into account. The spectral shape of light emission 
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was calculated in [22] with the help of the theory in [25] and good agreement 
with experimental data was found. Unfortunately, this theory neglects retar- 
dation effects in the electromagnetic field. Smolyaninov [24] believes that it 
is impossible to neglect retardation and that the agreement between experi- 
mental data and calculations declared in [22] seems to be fortuitous. Indeed, 
if we neglect the time which determines the difference in retardation of 
the electromagnetic wave emitted by different charges or currents, all pos- 
sible radiation disappears. For instance, the dipolar radiation only emerges 
if, expanding electromagnetic potentials in power series in the parameter to, 
we go to first order ujto in the approximation, where lo is the characteristic 
frequency of the classical current. 

At the present time the most plausible point of view is that the energy of 
a single tunneling electron is transformed into plasmon energy and that it is 
plasmons which then create photons [12]. However, the theory based on this 
idea in [12,22,25] is not only questioned by Smolyaniniv [20]. Xiao [27] also 
criticizes the approach developed in [12,22,25], arguing against any important 
role for plasmons in the photon emission process. Xiao tries to explain light 
emission from the tip-surface junction using a simplified Hamiltonian 

H = Ho eV cos out , (3.2) 

where Hq is the Hamiltonian of the electron moving through the potential 
barrier between the tip and surface and V is the applied bias. It is claimed 
that a system with this simplified Hamiltonian is able to describe emission 
from the tip-surface junction and that plasmons are not needed to explain 
emission. The arguments of both groups can be found in letters [28,29]. Even 
this brief review of approaches to the problem of light emission from the tip- 
surface junction clearly shows that the theory of such emission is far from 
settled. 

Since tunneling transitions are of a quantum nature, they occur randomly 
in time. Hence emitted photons accompanying this tunneling will likewise be 
distributed randomly in time. Put another way, time intervals between suc- 
cessively emitted photons will fiuctuate. At a current of InA, the average 
time between successively emitted photons is of the order of 10“^^ s. A PMT 
will reliably detect a train of photons even if the quantum efficiency of de- 
tection is about 10“^. The distribution of time intervals between emitted 
photons contains information about the dynamics of the tunneling junction. 
To our knowledge, this information concerning STM junctions has not yet 
been taken into consideration. In SMS, on the contrary, the temporal be- 
havior of fiuorescence measured in the photon counting regime is the main 
source of information about single-molecule dynamics. One- and two-photon 
counting methods in SMS serve as reliable sources of information concerning 
the dynamics of single molecules interacting with their environment [30,31]. 
The main goal of this review is to discuss achievements using these methods 
in the context of possible applications to scanning tunneling microscopy. 
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3.3 Role of Photon Counting 
in Time-Resolved Spectroscopy 

Over the last two decades the methods of photon counting have spread from 
the narrow confines of nuclear physics applications to the much broader area 
of common spectroscopy. The new generation of photomultipliers (PMTs) 
with microchannel plates and avalanche diodes has given additional impetus 
for applications of photon counting methods to far- and near-field optics. Pho- 
ton counting has important advantages over traditional measurement meth- 
ods for weak optical signals, such as modulated spectroscopy (e.g., the lock-in 
amplifier method and other analog photocurrent measurements). 

Upon photo-excitation an organic molecule is promoted to a non-equili- 
brium state from which it relaxes back to the ground state. Time-resolved 
photoluminescence spectroscopy is the method of choice for obtaining infor- 
mation on relaxation pathway (s). If the primary excitation event creates a 
vibrationally hot singlet (^i) state or a higher- lying Sn state, then relaxation 
can in principle involve vibrational cooling, internal conversion, intersystem 
crossing, readjustment of molecular selection, and finally, radiative or non- 
radiative decay to the ground state. In dense media, solvation effects as well 
as energy transfer and dissociation into a pair of charge carriers can also oc- 
cur. Which of these processes can be monitored by time-resolved spectroscopy 
depends on the time resolution of the experiment. 

Conventional methods with nanosecond time resolution yield the lifetime 
of an emitting state, whilst the preceding history of that state remains con- 
cealed, however. The advent of picosecond (ps) and femtosecond (fs) laser 
systems as excitation sources and of streak camera and up-conversion tech- 
niques on the detection side allow extension of the time resolution of fluores- 
cence spectroscopy into the sub-picosecond regime. This has opened the way 
to studying ultra-fast photophysical relaxation processes such as the solva- 
tion of an excited chromophore in solution, dephasing of optical excitations, 
and energy transfer in condensed systems. 

Photon counting makes it possible to apply statistical methods for pro- 
cessing measurement results. It means that we can measure kinetic processes 
in optics using pulse trains rather than single shots. Well known achievements 
of picosecond and femtosecond laser techniques give impressive examples of 
the high repetition of short laser pulses (less than 10 fs). The modern tech- 
nique of time-correlated single photon counting gives physical resolution less 
than 20 ps, but this value can be improved by an order of magnitude by 
employing special deconvolution methods. 
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3.4 Ensembles 

and Single-Molecule Spectroscopy Studies 

The huge progress in semiconductor nanophysics has led to studies of tiny 
molecular structures right down to single molecules. Electrical and optical 
properties of molecular nanostructures are under investigation around the 
world. Of course, the study of fast processes for electron system relaxation in 
restricted geometries is a key point. Effective combination of methods with 
high spatial and time resolution can achieve the desired result. 

Two main approaches are used in this field - confocal and near- field scan- 
ning optical microscopy. Both methods can combine with time-resolved spec- 
troscopy based on single photon counting. The main advantage of confocal 
microscopy is large value of optical signal (number of emitting photons) in 
contrast with NSOM approach. In case of NSOM the best spatial NSOM 
resolution (better that 20 nm) can be reached and the optical signal weak- 
ness can be overcome by increasing the time of experiments for good photon 
statistics assurance. 

In practice, it is often necessary to identify different analytes in a given 
sample. The authors of [32] reported measurements of spectroscopic proper- 
ties in a matrix at the single-molecule level. In these experiments, differences 
in the spectral emission properties of individual molecules of Rhodamine 6G 
and Texas Red in combination with two separate detectors were used to dis- 
tinguish the two dyes. 

To distinguish different fluorescent dyes or to identify tagged analyte 
molecules, it is not necessary to collect nearly as many photons as are needed 
for an exact lifetime measurement with the time-correlated single-photon 
counting (TCSPC) technique. The problem then consists in comparing raw 
data with certain well-defined fluorescence decay patterns measured with 
high precision. Hence, from a statistical point of view, one is confronted with 
a classification problem. An optimal classification method has already been 
developed in the framework of information theory and is described in detail 
elsewhere [33-36]. 

An important question arises for the time-resolved identification of single 
dye molecules in a matrix: how many photons are necessary to distinguish 
two or more dye labels with different lifetimes? A first approximation for 
the probability of misclassiflcation can be obtained by Cramer’s inequality 
[36]. Hence, in combination with an efficient maximum likelihood estimator 
(MLE) algorithm (3.4) [37], the identification of different dye molecules on 
the single-molecule level should be possible due to their characteristic fluo- 
rescence lifetimes. 

Zander et al. [38] first demonstrated the identification of single rhodamine 
dye molecules in water by their characteristic fluorescence lifetimes of 1.79 ± 
0.33 ns (Rhodamine B) and 3.79 ± 0.38 ns (Rhodamine 6G) excited by a 
frequency-doubled titanium/sapphire laser emitting at 514 nm. This was 
achieved using confocal microscopy in combination with the technique of 
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time-correlated single-photon counting. Using an efficient modified MLE al- 
gorithm, the classification probability was > 90% with 100 collected photons 
per burst. 

The detection volume is usually defined by the intersection of the focused 
excitation laser and the image of a spatial filter. Reduction of the detection 
volume improves the signal-to-background ratio by many orders of magni- 
tude without measurable photodestruction of the dye molecules under study. 
Poisson statistics predicts, for a concentration of less than 10“^^ M, that the 
number of molecules fluctuates predominantly between 0 and 1 in an ap- 
plied detection volume in the femtoliter region. Hence, the probability of two 
or more molecules being present in the detection volume simultaneously is 
negligible at concentrations below 10“ 

On the other hand, a small detection volume reduces the transition time, 
i.e., the measurement time. During the transition through the detection vol- 
ume, the dye molecule is excited by the laser light from the ground state 

into high- lying vibrational levels of the first excited state ^i, which un- 
dergoes rapid non-radiative internal conversion to low-lying S\ levels. The 
optical saturation limit is the maximum rate that a dye molecule can be cy- 
cled between and S\ and is dependent on the fluorescence lifetime Tf of 
the dye. 

Besides irreversible photodestruction, several depopulation pathways, 
such as intersystem crossing into the triplet state, compete with fluorescence 
emission thereby reducing the number of emitted photons. Since no photon 
is obtainable during the triplet lifetime of the excited dye, the triplet quan- 
tum yield and lifetime is a very important photophysical constant for single- 
molecule experiments. In other words, only those fluorescent dyes which 
exhibit very small triplet yields and short triplet lifetimes are suitable for 
single-molecule experiments. Consequently, most measurements on the single- 
molecule level in a matrix have been performed with rhodamine dyes. Ad- 
ditionally, rhodamine dyes emit up to 10^ photons before photodestruction, 
even in aqueous solution [39]. 

To sum up, in order to identify dye molecules in a matrix, both spectral 
and time-resolved detection methods can be used. A detailed description of 
modern optical studies of single molecules can be found in other sections of 
this chapter. 

3.5 Importance of Ensembles for Applications. 

Novel Ideas with Nanosize Ensembles 

There is considerable scientific and technological interest in design, synthesis 
and characterization of organic materials that possess unusual solid-state 
properties. Of practical interest are nanosize molecular ensembles such as 
linear molecular aggregates, regular and non-regular dendrimers and other 
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organic materials with restricted geometry in which cooperative effects have 
been observed. 

The interest to exciton dynamics in dendritic polymers (also known as 
Cayley tree or Bethe lattices) and molecular agregates (J- and H-, or mixed 
type) arises from their peculiar dimensionality [41] and can be satisfied by 
complex steady state and time resolved optical experiments. The dynamics 
of electronic excitations in the fractal dimensional systems (restricted or un- 
limited), the direct exciton transport in hierarchically built macromolecules 
[42], the interaction of these large molecules to each other are some of the 
fundamental problems to be modeled. Highly symmetric and perfectly hy- 
perbranched macromolecules with well controlled structures are good candi- 
dates for real exciton funnels. The approach [43] considers the high-branched 
wedges of dendritic molecules as light harvesting antennae, which pass the 
absorbed energy to initial core of the molecule playing a role of emitter. It is 
reasonable to suggest that such materials may be applied as a emitting layer 
of light emitting diodes (LEDs). 

Most of the polymers reported for electroluminescence (EL) applications 
have linear structures and little attention has been paid to using three- 
dimensional polymers as the active emitters in LEDs [43] . Compared to linear 
polymers with similar molecular weights, three-dimensional hyperbranched 
polymers have lower intrinsic viscosity [44,45] and better solubility. However, 
the EL properties of hyperbranched polymers have not been reported. In con- 
sidering the high functionality and desirable physical properties, combined 
with the robustness and ease of preparation of hyperbranched polymers, EL 
devices made from three-dimensional polymers are likely to have improved 
stability and efficiency compared with those made from general linear poly- 
mers. We try to extend two novel findings in this field: the light-harvesting 
effect and segregation of light-emitting Ln cations by dendritic shells. Both 
processes provide a way of enhancing luminescence by vectorial exciton mi- 
gration to the core and isolation of luminescent core molecules from other 
such molecules. Applications of nanosize molecular ensembles (dendrimers or 
aggregates) can thus provide some advantages for practical applications. 

By this approach, new features have recently suggested a method for 
electroluminescence study of thin films by applying a voltage between metal- 
covered NSOM tip and molecules/ensembles on a conductive substrate. Once 
again, photon counting can provide a wealth of information about the optical 
behavior of molecules/ensembles excited by tunneling currents. 

3.6 Nanoensembles 

3.6.1 Regular Dendrimers 

As a new class of materials, dendrimers or ‘starburst’ polymers have gener- 
ated great interest throughout the scientific community. Despite synthesizing 
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Fig. 3.3. Construction of typical regular dendrimer 



difficulties, a wide range of these new materials have been produced and 
characterized [46]. 

However, the study of their properties is still in its infancy (compared 
with classical polymers) and there remains a large gap in our knowledge 
which must be filled if the potential uses of these new materials are to be 
realized. Dendrimers are highly branched structures, produced either by di- 
vergent or convergent synthetic methodologies which allow the formation of 
monodispersed nanoparticles of controlled dimensions and specific surface 
fictionalization. 

The rationale for our studies on starburst dendrimers was based on the 
view that restrictions on the scale of single molecules of dimensions of the 
order 100 nm may result in novel electronic or optical phenomena not seen at 
other dimensions. For example, some correlated electronic effects on the den- 
drimer ‘surface’ at limiting generation would be observed such as fluorescence 
quenching or excimer formation. 



3.6.2 Hyperbranched Polymers 

Two distinct approaches to the synthesis of these hyperbranched polymers 
have evolved from the seminal work of Vogtle [47] and Denkewaller [48]. 
Tomalia’s ‘dendrimers’ [46] and Newkome’s ‘arborols’ [49] were independently 
developed through a process termed divergent synthesis. This approach is 
characterized by polymer growth emanating from a central core via an it- 
erative protection/deprotection reaction scheme. Polymer growth typically 
begins from a ‘core’ molecule which undergoes exhaustive reactions with com- 
plementary monomers having two or more protected branch sites. Removal of 
the protecting groups and subsequent reaction of the liberated reactive sites 
leads to first-generation polymers. Repetition of this reaction process leads 
to polymers of the desired molecular weight, molecular size and topology. 
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The convergent synthetic approach independently developed by Hawker 
and Frechet [50,51] begins at what will eventually become the outer surface 
of the dendrimer. Polymer ‘wedges’ are synthesized via sequential reactions 
and contain single reactive functionalities at their loci. Wedges are then at- 
tached to a polyfunctional central core to complete dendrimer formation. A 
rapidly increasing number of cascade molecules has been synthesized and 
investigated. 

3.6.3 J-aggregates 

J-aggregates of cyanine dye molecules are nanostructures with sizes inter- 
mediate between molecular crystals and isolated molecules. Their optical 
properties are distinctly different from those of their monomeric constituents 
[52,53]. Interactions between the molecules in J-aggregates give rise to delo- 
calized electronic excitations (Frenkel excitons), which are associated with a 
collective effect in the optical response [54-58] . 

Since they are delocalized over several monomer units, the excitations 
in J-aggregates are thought in turn to migrate over large distances (thou- 
sands of monomers) within the J-aggregates by an incoherent transfer mech- 
anism [59]. Although these exciton transport properties in J-aggregates are 
rather important for applications, they have been understood to a much 
lesser extent than those of molecular crystals [60]. The same is also true for 
the geometry of the molecular arrangement within the J-aggregate. Regu- 
lar brickwork [61], twisted brickwork [62] and several kinds of herringbone 
arrangements [63,64] have been proposed for J-aggregates. The unit cells of 
J-aggregate nanocrystals were assumed to contain one [65—67], two [68,69] or 
four molecules [70]. Moreover, different samples of J-aggregates consisting of 
the same dye molecules can contain a different number of molecules per unit 
cell depending on the preparation procedure [71,72]. 



3.7 Time-Correlated Single Photon Counting (TCSPC) 
Method 

The marriage of fluorescence microscopy with time-correlated single photon 
counting has contributed to our understanding of nanosize ensembles. 

Currently, most fluorescence microscopic imaging is performed as a mea- 
surement of emission or excitation intensity. These types of measurements 
have several limitations: 

• intensity-based fluorescence imaging can be difficult to quantify, 

• fast dynamic events (of the order of nanosecond to picosecond) in nanosize 
molecular systems cannot be studied, 

• the emitted fluorescence from such samples is often complex and repre- 
sents contributions from a number of molecular species which cannot be 
individually analyzed. 
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• autofluorescence and background fluorescence limit detection sensitivity. 

In contrast, the measurement of fluorescence lifetimes does not suffer from 
these limitations. Fluorescence lifetimes possess the added benefits of being 
independent of local intensity, concentration and photobleaching of the fluo- 
rophore. Fluorophores with similar spectra may exhibit significant differences 
in their lifetimes, and the same fluorophore may display distinct lifetimes 
in different environments. Because fluorescence lifetimes are not affected by 
scattering, measurements of fluorescence lifetimes provide more sensitive and 
quantitative information about complex structures. 

In most cases, periodic excitation pulses are used with a duty factor such 
that the luminescence between pulses manages to decay almost completely. 
Luminescence decay can be observed using an oscilloscope, either conven- 
tional or stroboscopic. Flash lamps or special discharges may be used as an 
excitation source. 

An important variant of pulsed methods is the statistical method of pho- 
ton counting. In this method, excitation is carried out with short pulses (1 ns 
or shorter) repeating at a frequency of 10 kHz or higher. Luminescence is 
detected by a photomultiplier operating in the photon counting mode. The 
luminescence intensity incident on the photocathode of the photomultiplier 
is weakened to such an extent that at most one photoelectron is knocked free 
from the photocathode after every excitation pulse. The time period between 
the excitation pulse and the moment of emission of this photoelectron is gov- 
erned by the probability of light emission. The number of excitation pulses, 
i.e., those in which the first photoelectron they release has been recorded, is 
considerably smaller than the total number of pulses, and the mean number 
of these first photoelectrons as a function of time reflects the luminescence 
decay. 

In the setup used for this method, the excitation pulse triggers the tim- 
ing circuit (start signal) and the first photoelectron stops this circuit (stop 
signal). The circuit records the period between the start and stop signals and 
transforms it to the amplitude of the output signal (time-amplitude con- 
verter) The amplitude of signals is analyzed using a multichannel amplitude 
analyzer. As a result, kinetic curves are obtained whose ordinate corresponds 
to the number of pulses in the given channel of the analyzer, i.e., the lumines- 
cence intensity, and the abscissa corresponds to the channel number which is 
proportional to time. 

As in other pulse methods, the resultant kinetic curves are convolutions 
of the true decay curves with the excitation pulse. Various methods have 
been proposed for deconvolution of these curves using computers and pro- 
ducing an elementary decay law upon ^-excitation. Owing to unavoidable 
noise, the deconvolution operation can generally only be carried out approx- 
imately. The problem is simplified if it is known in advance that the decay is 
exponential or can be approximated by a small number of exponents or by 
some predetermined function. 
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Fluorescence decays are determined in the Lebedev Physical Institute 
(Moscow) by time-correlated single photon counting (TCSPC), using a pi- 
cosecond dye laser for excitation at 575 nm as described elsewhere [44]. The 
instrumental response function possesses a FWHM of 0.7 ns while the typ- 
ical fluorescence lifetime for THIATS J-aggregates is 6 ns. All experiments 
are performed at very low excitation intensities (less then lOphotons/cm^ 
per pulse) to avoid exciton-exciton annihilation and photodegradation of the 
samples. 

The fluorescence lifetime can be obtained by analyzing decay curves within 
the framework of a multi-exponential model. The model tends to give values 
for in the range of 1.05-1.50 for different temperatures. When the data 
are analyzed as the sum of two exponentials, is usually below 1.30. 

3.8 Application of Photon Counting 
to Kinetic Measurements on Nanoensembles 

3.8.1 Dendrimer Rotation and Viscosity 

Time-resolved fluorescence anisotropy measurements can provide important 
information about the rotation of nanosize ensembles (e.g., dendrimers) and 
assess the role of the environment. Dendrimers with probes are shown in 
Fig. 3.4. In fluorescence anisotropy decay experiments, the time-correlation 
function of the emission transition is determined. For details concerning the 
experiment and data analysis, see [72,73]. 

The excitation transition probability of a molecule which absorbs inci- 
dent light is a function of the angle ol between the directions of its absorp- 
tion transition moment and the electric vector of the incident light, that is, 
E[a) — £"0 cos^ OL. Therefore, fluorescence emitted from molecules with an 
opto- anisotropic fleld ought to be observed as remarkably polarized light. 
However, we cannot observe this polarized light in a general optical system. 

The origin of the fluorescence depolarization phenomena can be classifled 
into four categories: the gap between absorption and emission transition mo- 
ments, intermolecular energy transfer. Brownian motion, and non-uniformity 
of fluorescent molecular orientation. Consequently, it is by the fluorescence 
depolarization method that we measure the degree of fluorescence depolar- 
ization and analyze particular physical factors based on the origin of the 
fluorescence depolarization phenomena in this experimental system. If there 
is no intermolecular energy transfer in an experimental system, we can extract 
only Brownian motion effects from the experimental emission anisotropy ratio 
data. In recent years, this technique has received attention as a way of mea- 
suring the Brownian motion of molecules and ensembles (or macromolecules) 
as this method has some useful properties, such as high sensitivity, molecular 
labeling, and real time measurement. 

Assuming that the dendrimer rotates as a spherical particle, the overall 
rotational correlation time obtained from fluorescence anisotropy decay would 
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Fig. 3.4. Examples of a few generations of dendrimers: Si}! (Me), Si®® (Me), and 
Sii(All) 
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lead to the coefficient of viscosity of the environment via the Stokes-Einstein 
relationship. 

The rotational reorientation of a fluorescent probe molecule that is bound 
to the dendrimer was investigated at the Lebedev Physical Institute with a 
view to evaluating the rotational mobility of the dendrimer. Polarized flu- 
orescence decay was measured with an Edinburgh Instruments spectrofluo- 
rimeter Model 199 by the time-correlated single photon counting technique. 
Fluorescence steady-state spectra were recorded using a multichannel optical 
analyzer OMA-2 (EG&G) with a A = 337 nm nitrogen laser. 

We measured the fluorescence anisotropy decay of pyrenyl attached to 
SiJKMe) and Si 2 g(Me) dissolved in cyclohexanol [74,75]. The kinetics of po- 
larized fluorescence and fluorescence anisotropy of pyrenyl in Si 29 (Me) is de- 
picted in Fig. 3.5. The fluorescence anisotropy decays with a lifetime of about 
2 ns reflecting the rotational reorientation of the fluorescent probe molecule 
attached to the dendrimer. It should be noted that in contrast to pyrenyl 
bound to Si 2 g(Me), the fluorescence for pure pyrenyl in cyclohexanol is com- 
pletely depolarized and does not show any kinetics in the time range investi- 
gated. The fluorescence anisotropy decay for pyrenyl bound to Si}! (Me) was 
found to deviate from that for Si 2 g(Me) only at early times, coinciding at the 
tail of the decay curve. Additional experiments with better time resolution 
are required to elucidate the difference in reorientation mobilities of Si^KMe) 
and Si 2 g(Me). 
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3.8.2 J-aggregate Radiative Lifetime and Structure 

Experimental data which can give unbiased information about the dimension 
of J-aggregates are those concerning the temperature dependence of the ra- 
diative lifetime of the exciton [76,77]. To date, this relationship has been 
determined for only two systems: J-aggregates of l,l'-diethyl-2,2'-cyanine 
(PIC) [78,79] and l,l'-diethyl-3,3'-bis(sulfopropyl)-5,5,6,6'-tetrachlorobenzi- 
midacarbocyanine (BIC) [80]. The temperature dependence of the fluores- 
cence lifetime has also been measured for J-aggregates of some other dyes 
[81,82]. However, an accurate analysis of the data within the framework of 
the problem of exciton dimensionality has not been undertaken. 

The effect of the temperature-induced population of exciton states has 
been calculated explicitly for a 1-D lattice by including the exciton-phonon 
interaction [83]. The transition dipole of a J-aggregate will also be limited by 
exciton scattering caused by static disorder. However, when these effects were 
considered in a theoretical analysis of the experimental data, it was impossi- 
ble, within the framework of 1-D excitons, to explain the greater than 8-fold 
increase in Trad in PIC aggregates when the temperature is increased from 1.5 
to 150 K [84]. An acceptable fit to the experimental data was obtained only 
within the framework of a 2-D exciton model, which allowed the authors to 
conclude that J-aggregates of PIC in solution possess a 2-D exciton coherent 
volume. 

Results of photon counting investigations into the temperature depen- 
dence of the radiative lifetime of molecular aggregates are shown in Fig. 3.6. 
The main goal of the current photon counting application is to present and 
analyze new data for the temperature dependence of the radiative decay 
time obtained for J-aggregates of THIATS, which differ strongly from the 
well-known PIC aggregates. 

To compare the temperature-dependent exciton dynamics of aggregates 
of THIATS to that of J-aggregates of PIC, it is convenient to consider the 
experimental data (radiative decay time) expressed in units of the inverse 
coherent length N~^. The temperature dependence of the inverse coherent 
length observed for J-aggregates for THIATS and PIC is shown in Fig. 3.7. 
The main difference between these two systems is the rate of increase of 
N~^ with increasing temperature up to 130 K. For aggregates of THIATS, 
no strong temperature dependence of was observed. 

To explain this slow increase of N~'^ with temperature, it is not necessary 
to introduce assumptions about a possible 2-D character of the excitons; 
a 1-D crystal possesses enough ‘dark states’ which can become populated 
upon increasing the temperature. This example demonstrates the relation 
between dimensionality and experimental data obtained by the correlated 
single photon method. 
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Fig. 3.6. Top: temperature dependence of the fluorescence decay time r {trian- 
gles with size corresponding to experimental error) obtained upon excitation at 
575 nm, and of the radiative exciton lifetime Trad {circles). Bottom: temperature 
dependence of the fluorescence quantum yield ^ for J-aggregates of THIATS in 
3/2-water/ethylene glycol, with excitation at 532 nm 



3.9 Photon Counting in SMS 

By pulse laser excitation we can study fast electronic relaxation of excited 
nanostructures even if we deal with their ensembles. However, slow relaxation 
in the ground electronic state of nanostructures cannot be studied by pulse 
excitation of the ensembles. Slow relaxation on time scales from microseconds 
to days is typical for impurity centers in polymers and glasses [85,86]. Re- 
laxation of the molecular ensemble in the ground electronic state will display 
some kind of averaged picture although experiments with single molecules 
reveal that the relaxation of each impurity center in an amorphous matrix 
displays individual features. 

In order to carry out spectroscopy on single molecules, it is clear that we 
need to decrease the concentration of the guest molecule. When this concen- 
tration decreases, two circumstances work in favor of SMS. Firstly, distances 
between guest molecules increase. If the distance exceeds the size of the laser 
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Fig. 3.7. Temperature dependence of inverse coherent length ^ for J-aggregates 
of THIATS and PIC. Solid lines are theoretical curves for 1-D excitons 
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Fig. 3.8. Appearance of statistical fine structure in the inhomogeneously broadened 
optical band when the guest molecule concentration is reduced [87] 



spot, the laser light will only excite a single molecule. Different types of 
near-field microscopy can reduce the size of the laser spot by a factor of 
ten compared with the diffraction limit. Secondly, when the guest molecule 
concentration is reduced, the inhomogeneously broadened optical band ac- 
quires a fine structure because in some frequency domains optical lines will 
be absent, as shown in Fig. 3.8. 

It is obvious that large inhomogeneous broadening helps to reveal well 
resolved optical lines belonging to single molecules. Polymers and glasses 
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Fig. 3.9. Train of photons emitted by a single molecule under resonant and off- 
resonant laser excitation 



serve as good solid matrices in SMS because of the large inhomogeneous 
broadening, reaching several hundred cm“^ in these matrices. 

Since a single molecule can absorb only one photon, an absorption mea- 
surement is not feasible because it is difficult to detect the disappearance 
of one photon in excited laser light. It is possible only if special modulation 
methods developed for measuring weak absorption are used [5]. 

However, if the molecule studied has fluorescence, the simple method of 
excitation fluorescence can be used [6]. At the present time, this is the main 
method in SMS. In this case the molecule is embedded in a cooled polymer 
matrix and irradiated by CW monochromatic laser light. The molecule con- 
sequently absorbs and then emits one photon. A train of photons emitted by 
a single molecule is shown in Fig. 3.9. 

In one-photon counting methods we count every photon from this train. 
The total number N{Aq) of photons counted by PMT for a predefined time 
interval, say one second, is proportional to the molecular absorption cross- 
section. It is a function of the detuning Aq — uj — cjq, where coq is the Bohr 
frequency of the molecular transition. 

At short counting times the number N{Aq) is small and will fluctuate 
because absorption and emission events occur at random times. The ratio of 
the fluctuation 6N{Aq) to the number N{Ao) of counted photons decreases 
as 5N{Ao)/N{Ao) = when N{Ao) increases. This means that the 

counting time is determined by the rate of photon emission and the desired 
signal-to-noise ratio. 

By scanning the laser frequency, we can change the detuning Aq and 
thereby measure the absorption line of a single molecule. An example of such 
an experiment is shown in Fig. 3.10. 

This exemplifies the case when one-photon counting allows us to measure 
the probability of absorption of a photon with energy Tiuo. This measurement 
takes a long time if we hope to reach a good signal-to-noise ratio. Conse- 
quently, the impurity center must be stable. However, in polymers we deal 
with impurity centers of another type. These can change their resonant fre- 
quency as time goes by. Figure 3.11 shows just such an example. The line 
position was measured in a 2.5 s laser scan. Here we witness rare jumps be- 
tween four spectral positions. 

A random function uo{t) like the one shown in Fig. 3.11 is called a spectral 
trajectory. Spectral trajectories of various types have been observed, with fast 
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Fig. 3.10. Excitation fluorescence line of a single pentacene molecule in a para- 
terphenyl crystal at 1.5 K [88] 
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Fig. 3.11. Jumps in the absorption line of a single terrylene molecule doped 
polyethylene at 1.5 K [89] 



and slow jumps between different spectral positions [8,89-93]. It is obvious 
that the spectral trajectory is a random function of time which contains 
information about some kind of relaxation processes. These processes are 
characterized by probabilities. Below we discuss the way these probabilities 
can be determined from the random trajectory. 



3.10 Intermittency in Single-Molecule Fluorescence 

The ‘spectral trajectory’ concept can be generalized by introducing the more 
general concept of a ‘quantum trajectory’. Indeed any transition in a quan- 
tum system like a single molecule happens randomly in time. Therefore, any 
physical variable connected with the transition will fluctuate on the time 
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Fig. 3.12. Fluorescence from single FAD [92] 




scale. The frequency fluctuations shown in Fig. 3.11 are only one of many 
possible quantum trajectory variations. 

Figure 3.12 shows an example of quantum trajectory of another type when 
the trajectory describes fluorescence intensity fluctuations. These fluctuations 
resemble a random telegraph signal. This trajectory enables one to monitor 
the chemical reaction on a single-molecule level. This has been done by Lu, 
Xun and Xie in [92]. 

These authors studied enzymatic turnovers of oxidized single cholesterol 
molecules in real time by monitoring luminescence from the enzyme’s flu- 
orescence site, flavin adenine dinucleotide (FAD). FAD can emit light with 
a wavelength of 520 nm when in oxidized form. However, the fluorescence 
disappears after its reduction by cholesterol. Single FAD was irradiated by 
He-Cd laser light at a wavelength of 442 nm. Fluorescence was collected with 
an inverted fluorescence microscope by raster-scanning the sample. At large 
cholesterol concentrations, reduction and oxidation processes were in dynamic 
equilibrium. In the molecule ensemble, they observed equilibrium concentra- 
tions of both forms. However, on the single-molecule level, oxidation and 
reduction happen randomly in time, as can be detected via the appearance 
and disappearance of fluorescence. This random process is described by the 
quantum trajectory in Fig. 3.12. The lengths of ‘off’ and ‘on’ intervals take 
random values. 

Lu et al. processed the measured quantum trajectories statistically, includ- 
ing up to 500 random off- and on-intervals. The histogram of on-intervals is 
shown in Fig. 3.13. The distribution of on-intervals can be described by the 
following function: 

pit) — [ exp(-at) - exp(-6t)l , (3.3) 

0 — a '■ 

where a = 2.9 s“^ and b = 17 s~^. This is a probability density for finding 
the molecule in oxidized form. It is characterized by two rate constants a 
and b which describe oxidation and reduction rates. This example shows how 
we can find the rate constants of the process from a quantum trajectory of 
random character. 
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3.11 Spectroscopy of Single Molecule 
with Triplet Level. Photon Bunching 

Figure 3.14 shows the energy diagram for the lowest electronic levels of a 
typical organic molecule. As a rule, these molecules have a triplet level be- 
tween the ground state singlet level and the first excited singlet level. The 
probability 7 st of transition between the triplet and the ground state singlet 
level is small so the triplet lifetime ranges from microseconds to seconds. The 
long-lived triplet level dramatically alters the train of emitted photons. 

When the molecule irradiated by CW laser light lies in the singlet state, it 
can jump between the ground and excited singlet levels absorbing and emit- 
ting photons at random times. The probability Tts of intersystem crossing in 
the excited electronic state is quite large and in fact is comparable with the 
probability 1/Ti of spontaneous photon emission. Therefore, the intersystem 
quantum yield r] — TiIts/(1 + FiTts) is of the order of unity. When the 
excited molecule jumps to the triplet state 2, emission stops and radiation is 
absent while the molecule remains in the triplet state. After a random time 
of the order of l/ysT, the molecule jumps to the ground electronic state and 
then emission appears again, as shown in Fig. 3.15. 
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Fig. 3.14. Typical tree-level diagram for an organic molecule and rate constants 
which determine the dynamics of the molecule 
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Fig. 3.15. Photon bunching in photon emission due to the influence of the triplet 
electronic level 
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Fig. 3.16. Intermittency of single-molecule fluorescence 



In fact, the train of photons consists of bunches and this phenomenon 
is called photon bunching. In accordance with Fig. 3.15, on-intervals when 
photons are emitted alternate with off-intervals toff when there is no emission. 

Let us now consider how to calculate probabilities for finding definite 
values of on- and off-intervals. Let the PMT count photons for a time interval 
t. If t is longer than the average interval between two photons, the PMT signal 
will be constant when fluorescence is emitted. When the molecule jumps 
to the triplet state, fluorescence disappears. We will measure a quantum 
trajectory, as shown in Fig. 3.16. 

This quantum trajectory, which looks like a random telegraph signal, is 
related to the fluorescence intensity. From the trajectory, we can find average 
Ton and Toff intervals by substituting the measured ton and toff values for ti 
in the formula 



r 



1 

N 




(3.4) 



A quantum trajectory can also characterize other physical variables like fre- 
quency, polarization, and so on. In all cases the quantum trajectory is a 
random function of time. Any trajectory is an irreproducible characteristic 
of the molecule and cannot only be calculated theoretically with the help of 
the Hamiltonian of the system. Indeed, quantum mechanics allows us to cal- 
culate probabilities but not trajectories. Let us consider how we can calculate 
the distributions of on- and off-intervals. 

It is obvious when a molecule emits light that it is in singlet electronic 
states. When there is no emission the molecule occupies a triplet electronic 
level. The dynamics of singlet states is described by the following rate equa- 
tions for probabilities of finding the molecule in the ground and excited singlet 
states: 



h = - 



Ti 



T jTts “k k 



pi -k kpo , 



Po 




Pi - kpo . 



(3.5) 
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Figure 3.14 explains the physical meaning of all the rate constants. It is 
obvious that 

Pon = Pi + Po (3-6) 

is the probability of finding the molecule in a singlet electronic state, i.e., in 
an on-state. Summing the two equations, we find 

Pon = -A^SPl* (3.7) 

At times exceeding l/Trs, we may put pi = 0 in (3.5). We then find from 
the first equation the relation 

Pi — T , /T V r — 

1 + kTi -h i tsm 

Making use of (3.7) and (3.8), we find 

~ 1 1 /^o7 I r \T^ ^ p ^ ' ( 3 - 9 ) 

1 + ( 2 fc + FtsiTi 1 + FtsTi 

With the help of this equation, we can transform (3.7) to the form 



where 



1 _ FtsTi 
^on 1 + 7 ts7i 



The ratio here describes the intersystem crossing quantum yield. A solution 
of (3.14) is the function 



This is the probability density for finding a random value of ton equal to 
t. It is a distribution function of on-intervals which can be measured after 
carrying out statistical processing of the quantum trajectory. It is obvious 
that Ton determined by (3.11) is an average value of on-intervals. It depends 
linearly on the laser light absorption intensity k. 

When the molecule goes into the triplet state, emission stops. This means 
that the molecule is in the off-state. Hence we may write p 2 = Poff- The 
dynamics of the off-state is described by the following equation: 

P2 = ~7stP2 , (3.13) 



Solving (3.13), we find 
^ofF = 7 st exp(-7sT^) • 



This function describes a distribution of off-intervals. The value Toaf deter- 
mined by (3.14) is an average value of off-intervals and is independent of the 
intensity. 
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3.12 Two-Photon Correlators 



By counting photons, we measure a quantum trajectory and then by process- 
ing the measured trajectory, we can find probabilities. A question arises: can 
we avoid the measurement of a random function like a quantum trajectory 
and measure transition probabilities directly by counting photons in real- 
time experiments? In fact, it is possible if we use the so-called two-photon 
counting method. As the name suggests, in these methods we count pairs of 
photons. The photon pair count rate is called the two-photon correlator. The 
theory for two-photon correlators has been developed recently in [93-95]. 

Figure 3.17 shows photons emitted by a single molecule irradiated by 
CW laser light. Let us consider pairs of photons with a definite time delay 
between them. The figure shows that there are two types of photon pair. For 
instance, we can count in the so-called start-stop regime by taking into ac- 
count only those pairs consisting of consecutively emitted photons like pairs 
(2,3) and (12,13). The count rate s{t) of such pairs is called the start-stop 
correlator. However, there is another possibility, namely, to count all pairs of 
photons with a delay, including additional pairs like the pair (4,6) with one 
intermediate photon 5. The count rate for all pairs with an arbitrary number 
of intermediate photons is called the full two-photon correlator. The count 
rate for photon pairs with a definite delay is proportional to the probability 
of finding this type of pair. It can be calculated quantum mechanically. The 
theory of two-photon correlators for two- and three-level molecules is devel- 
oped in [93]. The theory has been generalized to molecules interacting with 
huge numbers of tunneling systems as happens in a polymer or glass matrix 
[94,95]. 
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Fig. 3.17. Two types of photon pairs 



Start-Stop Correlator. Let us begin by considering the simplest case re- 
lating to a two-level molecule interacting with a transverse electromagnetic 
field. Let Wq, Wi and be the probabilities of finding the molecule in the 
ground electronic state, in the excited electronic state and again in the ground 
electronic state after spontaneous emission of a photon with wave vector fc. 
If we take t = 0 a,s the time of recording the photon 2 shown in Fig. 3.17, we 
find that VFo(O) = 1 and IFi(O) = Wk{0) = 0. Due to laser excitation, these 
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probabilities will change. In time r we find an emitted photon. It is obvious 
that 

= '^Wk{T)dT (3.16) 

k k 

is the probability of photon 3 being emitted in the time interval from t = 0 
to t = dr. The count rate for pairs like the pair (2,3) is proportional to the 
probability Wk{r), which refiects the molecular dynamics. It is shown in 
[93] that probabilities are governed by the conservation law 



Wo + Wi + J2^k = l- (3.17) 

k 

This means that one of the probabilities need not be considered. It was also 
shown that 

^^ = Y^Wk{t) = s{t), (3.18) 

k 

and that the probability W\ can be found from the following equations 

li^io = -i(^ - ir)w,o - x(iyo - i^i) , 

1^01 = i(^ + '^r)Wm - x(Wo - Wr) , 

W, (3.19) 

W, = -x{Wio + Wox)-^. 

-Ll 

Wo = x{Wio + Woi)^ 

where A is the difference between the laser and resonant frequencies, F = 
l/2Ti, and y = jE • d//i is the Rabi frequency. Equation (3.18) can easily 
be derived by summing the third and fourth equations in (3.19) provided we 
take into account (3.17). The function s{t) is called the start-stop correlator. 



Full Two-Photon Correlator. If we count all pairs with a time delay 
and take into account the pair (4,6) with one intermediate photon (Fig. 3.17) 
and also other pairs with various numbers of intermediate photons as shown 
in Fig. 3.18, we arrive at the full two-photon correlator. The first term in 
Fig. 3.18 relates to the start-stop regime where there are no intermediate 
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Fig. 3.18. Summation of photon pairs with various numbers of intermediate pho- 
tons 
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photons. The second term relates to the pair with one intermediate photon 
fc, and so on. The counting rate in this case equals [93] 

k kk' kk'k" ^ V k 

= p{t) = ^, (3.20) 

where pn satisfies the optical Bloch equations 

pio = -i{A - ir)pio - xipoo - Pn ) , 

Poi = i(^ + iT)pio - xipoo — Pn ) , 

Pn = -xiPio T Poi) - ^ 5 (3.21) 

Poo = x(pi 0 + Poi) + ^ • 

Equation (3.20) determines the full two-photon correlator. 

Using (3.19) and (3.21), we can a find simple relation between the start- 
stop and full two-photon correlators [93] 

t 

p{t) = s{t) + y* 5(t - x)p{x) dx . (3.22) 

0 

Solving this integral equation by an iteration procedure, we find the solution 
in the form of a sum with an infinite number of terms. These terms relate to 
the terms of the sum in (3.20). For instance, the first term 

= (3.23) 

k 

describes the counting rate in the start-stop regime, whilst the second term 

t 

Wkk' — f s{t-x)s{x)dx (3.24) 

kk' Q 

describes the contribution to the full two-photon correlator from photon pairs 
like the pair (4,6) with one intermediate photon 5, as shown in Fig. 3.17. It is 
obvious that the nth term of the sum describes the contribution from pairs 
with n intermediate photons. 

Figure 3.19 shows start-stop and full correlators calculated in [93] with the 
help of (3.19) and (3.21) for various values of the Rabi frequency. The start- 
stop correlator aproaches zero on time scales similar to the fluorescence time 
Ti . It is impossible to research the slow dynamics typical of chromophores in 
polymers and glasses by means of the start-stop correlator. Slow dynamics 
in polymers and glasses must therefore be probed by the full two-photon 
correlator. 
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Fig. 3.19. Dependence of correlators s{t) and p{t) on time delay between two pho- 
tons at Z\o = 0 and = 0-1 (I): 0-2 (2) and 1 (3) 



Because we cannot find photon pairs with very short delays, both types 
of correlators approach zero at these short delays. This phenomenon is called 
anti-bunching. The full two-photon correlator is a constant at large times. 
This means the absence of any correlation between two photons with large 
delay. The second photon ‘forgets’ about the first photon if the time delay 
is much longer than the fluorescence time Ti. However, if some kind of slow 
process exists in the system and this process influences absorption, the prob- 
ability of flnding the second photon will feel this slow process. It will be 
reflected in the long-term behavior of the full correlator. In the next section, 
we consider just such an example. 

3.13 Full Two-Photon Correlator 
for Single Molecule with Triplet Level 

We have already seen in Sect. 3.11 that the triplet level strongly influences 
the photon train emitted by a single molecule irradiated by CW laser light. 
Let us calculate the full two-photon correlator for such a single molecule. In 
order to describe this case, we need to add a flfth equation to the four we 
already have in (3.21). This equation should describe a transition from the 
excited singlet level to the triplet level and a transition from the triplet level 
to the ground electronic state in accordance with the diagram in Fig. 3.14. 
The rate of light-induced transitions is 

where Aq is the difference between the laser frequency and the Bohr frequency 
of the singlet-singlet resonant transition. The expression for the full two- 
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photon correlator for such a molecule was found in [93]: 



p{t) = 



7st 



k 



Ti U-Rl 



+ 1 



where 



7st 



3-(7o-i?o)^ 



7o — i?o / 2i?o 

7ST ^ 0~(7o+-Rq)^ 



- 1 - 



7o + Rq 



2Ro 



r + 2fc + 7sT D 
7o = , ^0 = 



r + 2k — 7sT \ 









(3.26) 



(3.27) 



The temporal dependence of this correlator is shown on a logarithmic time 
scale in Fig. 3.20. 

Specific features of this correlator have a definite physical meaning. The 
value of the ordinate at small delays, say 10~^ s, corresponds to the count 
rate of photon pairs in a bunch. The value of the ordinate at long delays, 
say 1 s, corresponds to a situation where we count photon pairs consisting of 
photons from various branches. In other words, it is the count rate for photon 
branches. 

Exponential decay looks like a smooth step over one order of magnitude 
of time on the logarithmic time scale. The two steps on the curve in Fig. 3.20 
can be related to the exponential growth of emission with time delay 10“^s 
and exponential decrease in emission with time delay 10“^s. These two times 
correspond to the rate of population of excited electronic level 1 and to its 
rate of depopulation due to the transition to triplet level 2 in Fig. 3.14. This 
example clearly shows that the full two-photon correlator enables us to in- 
vestigate slow relaxation in the system. This is a very important feature 
of the full two-photon correlator because it can be used to investigate slow 
relaxations in polymers and glasses. 




Fig. 3.20. Full two-photon correlator for a single molecule with triplet level cal- 
culated by means of (3.26) with T 2 /T 1 = 10“^, ysxTs = 10“®, FtsTl = 9 and 
xT2 = 2 X 10"^ 
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3.14 Tunneling Dynamics of Polymers Probed by SMS 

It is well known that low temperature properties of polymers are due to the 
existence of so-called two- level systems (TLS) [96,97]. They are described by 
a double- well potential as shown in Fig. 3.21. 

The change in the energy splitting of the double-well adiabatic poten- 
tial resulting from electronic excitation of the chromophore is the param- 
eter of a chromophore-TLS interaction of the Franck-Condon type. This 
chromophore-TLS interaction determines anomalous temperature broaden- 
ing of spectral holes burned in the chromophore optical bands. Such broad- 
ening depends linearly on temperature in the mK temperature range [98] in 
contrast to the temperature broadening in crystals which follows either 
or Arrhenius-type temperature laws [99]. 

In accordance with Fig. 3.21, the optical band of the chromophore in- 
teracting with a TLS consists of two optical lines with separation A in the 
frequency scale. 

These doublets could serve as direct evidence for the existence of TLSs 
in polymers. Nevertheless, they have never been observed in hole-burning 
spectroscopy of polymers because there is a distribution over frequency sepa- 
rations in molecular ensembles. The first direct observation of such doublets 
was made in SMS [100]. 

There are seas of TLSs in every polymer. Their relaxation rates R are 
distributed over a very wide time scale, covering more than ten orders of 
magnitude from microseconds to days. It is clear that the slow relaxation 
of TLSs must manifest itself in the full two-photon correlator. The theory 
for two-photon correlators [93] was therefore generalized to take into account 
the interaction of a single chromopore with a sea of TLSs. This was done in 
[94,95]. 

The theory was built up by taking into account the following hierarchy of 
rates 

1/Ti > A: > A, a, B, 6 , (3.28) 




Fig. 3.21. Tunneling system of a polymer in the ground and excited electronic 
states of a chromophore 
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as shown in Fig. 3.21. Here T\ is the chrornophore lifetime, k is the rate of 
light-induced transitions and A, a, and b are rates of tunneling in the 
ground and excited electronic states of a chrornophore. As shown in [94,95] 
tunneling rates in the excited electronic state of the chrornophore are pro- 
portional to kBTi and kbTi. Hence if A, a, B and b are of the same order, 
we can neglect light-induced tunneling in the excited electronic state. This 
conclusion is correct for all the TLSs inherent in a polymer. In this case the 
formula for the full two-photon correlator takes the form 

p{t) = k{Ao, t, T) (l - , (3.29) 

where the function fc(Z\o,t,T) of detuning Aq = u — ljq and temperature 
T describes the light absorption coefficient of the chrornophore interacting 
with TLSs, some of which had not reached thermal equilibrium by the time 
t when the second photon of the photon pair was emitted. For the case when 
a chrornophore interacts with a single TLS, this function takes the form 

k{Ao, t, T) = [1 - p{t, T)] k(Ao) + p{t, T)k{Ao - A ) , (3.30) 

where k{Ao) is a Lorentzian described by (3.25), and 

Pit, T) = /(T) + [/)(0) - /(T)] expi-Rt) (3.31) 

determines the probability of finding the TLS in the excited state with the 
energy e. Here R = A + a. It is obvious that the probability in thermal 
equilibrium is p(co,T) = f{T) = [exp(e:/fcT) + 1]“^. 

How will TLS relaxation manifests itself in the two-photon correlator? 
The full two-photon correlator calculated using (3.29) and (3.30) is shown 
in Fig. 3.22. For excitation at frequency uq we must take p(0) = 0 in (3.31) 
because at this moment of time the system occurs in the left-hand well of the 
double- well potential shown in Fig. 3.21. For excitation at frequency cjq -j~ A ^ 
we must take p(0) = 1 in (3.31). 




logCt,G) 



Fig. 3.22. Temporal dependence of full two-photon correlators and excited elec- 
tronic state of a guest molecule for excitation at frequencies coo (1) and cjq + A (2). 
/ = 0.4, A — S/T 2 , Ti = 10“^ s and R = lO^^s"^ 
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Fig. 3.23. Dependence of two-photon correlator on detuning for the same values of 
parameters as in Fig. 3.22. Left: at short delays with tR 1. Right: at long delays 
with tR^ 1 



The correlators shown in Fig. 3.22 increase their values at times of order 
Ti and decrease their values at times of order 1/i? which characterizes TLS 
relaxation. The ratio of the correlators 1 and 2 at long delays is approximately 
equal to the ratio (1 — /)//. Comparison of Figs. 3.22 and 3.20 reveals that 
relaxation due to singletMriplet transitions and relaxation due to tunneling 
transitions manifest themselves in the full two-photon correlator in a similar 
fashion. Let us consider now the correlator described by (3.29) and (3.30) as 
a function of the laser frequency. 

The result of our calculation is shown in Fig. 3.23. These pictures of the 
two-photon correlator with short and long delays between photons are conve- 
nient when discussing spectral pictures obtained in one-photon experiments 
with short and long laser scans. 

For fast laser scans, when t^cR <C 1, we observe a single optical line jump- 
ing randomly between two spectral positions from one laser scan to another. 
The spectral picture is similar to that shown in Fig. 3.23 (left). Which line 
will appear is determined by a random factor, depending on whether the TLS 
is found in the left- or right-hand well. 

For slow laser scans, when tscR ^ 1, we observe just two optical lines 
in every laser scan as shown in Fig. 3.23 (right). The line intensities are pro- 
portional to the probabilities 1 — / and / of finding the TLS in the left- or 
right-hand well. For slow scans, we cannot detect relaxation if it is faster than 
the duration of the scan. The duration of the laser scan determines the time 
resolution of the experimental setup. In STM measurements of the I~V char- 
acteristic the time of the current measurement plays the role of in optical 
experiments. However, it is worth noting that in two-photon experiments the 
shortest delay t between two photons will characterize the time resolution of 
the setup. This delay can be several orders of magnitude shorter than the 
best time resolution in experiments with electric current. 

The physical picture discussed in connection with Fig. 3.23 will manifest 
itself in the train of emitted photons. This train of photons is shown in 
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Fig. 3.24. Train of photons emitted by a single molecule excited with various laser 
frequencies 



Fig. 3.24. When a single molecule is excited at frequency cjq? a pause appears 
in the emission at the random time when the TLS jumps from the left- to 
the right-hand well. The chromophore in the right-hand well has a different 
resonant frequency, ujq-\- A. Therefore, if the laser frequency is not tuned to 
the new value of the resonant frequency, the emission disappears and the lower 
train of photons in Fig. 3.24 will be absent. In this special case Fig. 3.24 will 
be similar to Fig. 3.15 which shows the train of photons emitted by a single 
molecule with a triplet level. If we tune the laser frequency to uoq + A, we 
will observe bunches from the lower train in Fig. 3.24. It is important to note 
that pauses in one train of photons exactly coincide with bunches in the other 
train, and vice versa. This is a telling feature of the chromophore interacting 
with a TLS. 

3.15 Relation Between One- and Two-Photon 
Counting Methods in a More Complex Case 

In the preceding section we considered the relationship between one- and two- 
photon methods for the simplest case when a chromophore interacts with a 
single TLS. The relevant spectral trajectory will display random jumping 
between two spectral positions. It is obvious that this physical model cannot 
describe the spectral trajectory shown in Fig. 3.11 where the optical line 
jumps between four spectral positions. This trajectory can be understood 
with the help of a model which describes a chromophore interacting with two 
TLSs. 

Two TLSs have four levels. Therefore, instead of (3.30), we should take a 
coefficient k with the form 

fc(Z\o, T) = [(1 - />)(! - p')fc(^o) + P{1 - p’)k{Ao - A) 

+(1 — p)p'k{A{) — A') + pp'k{A() — A — A^)\ , (3.32) 

where p and p' describe the probabilities of finding the first and second TLS to 
be excited. They are described by equations similar to (3.30). The Lorentzian 
fc(Z\o) is described by (3.25). The frequency of the laser light determines 
initial conditions for the probabilities p and p' . For example, if we excite a 
single molecule with the detuning Z\q = cj — cx;q — 0, only the first Lorentzian 
in (3.32) works. Consequently, the system consisting of a chromophore with 
two TLSs occurs after emission of a photon in the quantum state with p(0) = 
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p'(0) = 0. If we excite the system via the second Lorentzian in (3.32) we must 
take p(0) = 1, p^(0) = 0, and so on. 

Figure 3.25 shows the results of calculating the full two-photon correlator 
by means of (3.32) upon excitation via the first, second, third and fourth 
Lorentzian, i.e., for the initial conditions p(0) = p'(0) = 0, p(0) = 1 and 
p'(0) = 0, p(0) = 0 and p'(0) == 1, and p(0) = p'(0) = 1. The curves in 
Fig. 3.25 were calculated for three different time delays between photons. 
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Fig. 3.25. Full two-photon correlator upon excitation via various optical lines and 
at various time delays between photons t = 10~^ s (top), t = Is {middle), t = 10^ s 
(bottom). TLS parameters are as follows: £ = 0.5 K, R = 6 x 10^s“^, T 2 A = 4, 
e' = IK, R' = T 2 A' = -6, L7K 
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For t < 1/i? < 1/i?' we obtain Fig. 3,25 (top left). The relaxation in 
both TLSs is very slow and relaxation takes longer than the delay between 
photons in the pair. In this case the correlator as a function of the frequency is 
described by a single Lorentzian. Which of the four Lorentzians will manifest 
itself depends on the value of the laser frequency. As already shown in the 
preceding section, the spectral picture for the two-photon correlator can be 
related to a spectral trajectory measured in one-photon experiments with 
laser scans. If the laser scan time satisfies the inequality tsc < l/R < l/R' 
in one-photon experiments, we obtain a spectral trajectory with jumping 
between four spectral positions like that shown in Fig. 3.11. 

For l/i?<t<l/i?'we obtain Fig. 3.25 (top right). In this case relaxation 
in one TLS occurs faster than the laser scan. The second photon does not 
feel this fast relaxation. We therefore see the doublet related to the TLS 
with fast relaxation upon excitation via 1,2 lines or 3,4 lines. In one-photon 
experiments, if 1/i? < tgc < l/R' we will see a doublet of optical lines jumping 
randomly on the frequency scale from one laser scan to another. Figure 3.25 
(top right) shows this jumping. 

For l/i?<l/i?'<twe obtain Fig. 3.25 (bottom). In this case, we will see 
all four optical lines in each laser scan just as we see the same stable spectral 
picture in one-photon counting experiments with slow laser scans. 

3.16 Local Dynamics 
of Individual Molecules in Polymers 

The spectral trajectory shown in Fig. 3.11 is just one example from the wide 
variety of spectral trajectories measured in experiments with single terrylene 
molecules in polyethylene crystals. In fact every guest molecule has its own 
individual spectral trajectory even when used to dope crystals [88,89]. It 
reveals the wide range of local environments even for guest molecules existing 
in crystals. 

It is worth noting that some of the differences in spectral trajectories can 
result from the random character of each trajectory. Differences of this type 
do not reflect any differences in the local environment of guest molecules. 
Only differences in probabilities measured in experiments can prove that 
there really is a variety of local environments. 



Line Width Distribution. The half- width 2 /T 2 of a zero-phonon line de- 
pends on temperature and exceeds the fluorescence rate 1/Ti. This half- width 
results from electron-phonon and electron-TLS coupling. One might not ex- 
pect such variety in this type of coupling for guest molecules in crystals. 
However, SMS revealed that a distribution over half-widths exists even in 
the anthracene single crystal [101]. Such distributions over line widths of 
single molecules are larger in polymers. An example is shown in Fig. 3.26. 




3 Photon Counting Methods in STM and SMS 



83 




Homogeneous linewidth, /’(MHz) 



Fig. 3.26. Distribution of the 176 line widths of terylene molecules in a polyethylene 
matrix at T = 1.8 K [90] 



For terylene molecules doping the Shpol’ski matrix, poly (vinyl butyral), 
poly (methyl ethacylate) and polystyrene, such distributions can be found 
in [101]. The existence of distributions over line widths means that optical 
dephasing of single molecules is determined by the interaction of the guest 
molecule with localized phonons and TLSs in the immediate neighbourhood 
of the molecule. 



Temporal Behavior of Full Two-Photon Correlators. The temporal 
behavior of the two-photon correlator characterizes probabilities of relaxation 
processes. Therefore differences in two-photon correlators can be related to 
differences in the local environment of individual guest molecules. In ex- 
periments with single molecules, an autocorrelation function (AF) is 

measured. The AF is related to the fluorescence intensity I{t) and can be 
expressed via the full two-photon correlator as follows: 



p{to) 

p{oo) 






i-oc {I{t)I(t + Oo)) 



lim 

t—^oo 



{I{t)I{t + to)) 

mf 



(3.33) 



The infinite time is determined in practice as the time when all relaxation is 
over. The AF thus approaches unity at large time delays to- 

In practice, it is more convenient to use the AF than the two-photon 
correlators, because all factors relating to experimental conditions beyond 
the control of the researcher are excluded by taking a ratio. However, if the 
count rate of photon pairs with large delay is very small, the value of p(cx)) 
is very small as well. Therefore, the value of the AF will be found with large 
error. This type of situation emerges when we count photons emitted by a 
single molecule with triplet level. If the probability 7 st of the triplet-singlet 
transition is small, the value of p{oc) will be small and it will be measured 
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log( t, ns) 

Fig. 3.27. Autocorrelation functions for three terylene molecules in polyethylene at 
T = l.SK [90] and theoretical curves calculated using (3.29), (3.30) and (3.32) [31] 



with large inaccuracy. This will dramatically influence the ratio in (3.33). 
The situation can become worse still if a molecule has undergone spectral- 
hole burning. In this case emission disappears at large times. It is then unclear 
what value of p{oo) we should use in (3.32). 

Figure 3.27 shows the AF of three individual terylene molecules in a 
polyethylene matrix. They are quite different. Since the formula for the AF 
involves probabilities, the difference can be related to a difference in the 
immediate neighbourhoods of the three guest molecules. The curves show 
the difference in local dynamics for three individual molecules. 

On a logarithmic time scale, a smooth step by one order of magnitude 
means an exponential decay. Therefore, we may say that in the immediate 
vicinity of molecule 1, there are two TLSs with relaxation rates 10^ s~^ and 
ls~^, whilst there is one TLS with relaxation rate 0.1 s~^ and 100 s~^ in the 
vicinity of molecules 2 and 3, respectively. 

3.17 Conclusions 

3.17.1 Fast Photon Counting Measurements 
vs. CW Measurements 

The time dependence of light emission from photoexcited molecules is stud- 
ied extensively in physics and chemistry to understand the interactions and 
reactions of those molecules. However, until now, nearly all such experiments 
have studied large numbers of molecules (ensembles) simultaneously, thus 
providing data that represents the average of the behavior of the individual 
molecules. This severely limits the study of these phenomena. 

The new approach involves a matrix containing fluorescent molecules and 
exposed to focused laser light. The molecules are spatially separated in such 
a way that only one is exposed to the light at any time. When the sample is 
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scanned with the laser, a fluorescence image can be obtained in which each 
bright spot corresponds to a single molecule. 

By continuously sweeping the laser polarization from 0 to 90 degrees while 
examining a single molecule, the orientation of the dipole of the molecule was 
determined with high precision and its changes in time recorded. 

In some cases, the light emission of a single molecule is observed to switch 
on and off as the molecule is presumed to convert between light-emitting and 
non-light-emitting or Mark’ states. By measuring the dipole orientation before 
and after the ‘dark state jump’, researchers have been able to show that the 
molecule remains stationary as it switches state. In other cases researchers 
observe changes in the polarized emission signal indicating that the molecules 
experienced ‘rotational jumps’ in which they reoriented themselves. 

The fundamental diflference between optical studies of molecule ensembles 
(molecular aggregates or macromolecular dendrimers) and optical studies of 
single molecules should be emphasized. In the first case it possible to apply 
the technique of time-correlated single photon counting to exciton kinetic 
measurements under short laser pulse excitation. This is to be contrasted 
with the second case where spectral measurements and molecular jumps can 
be studied under CW laser excitation. In both cases it is possible to apply 
microscopy methods. 

3.17.2 Combination of TCSPC and NSOM: New Advantages 

At the Lebedev Physical Institute in Moscow, a unique setup has been in- 
stalled. The new setup combines the facilities of two techniques: NSOM and 
TCSPC. The TCSPC system is based on a femtosecond Ti-Sapphire laser 
pumped by an Ar ion laser with T-format online detection. 

In this way, doped or pure thin film luminescence kinetics can be studied 
with maximal for optics limit. It gives information about exciton dynamics 
between nanosize species or between macromolecule fragments. This method 
is extremely important for studying the above-mentioned antenna and shell 
effects associated with dendrimers. 



3.17.3 Kinetics of Energy Transfer Between the NSOM Tip 
and Single Molecules/Ensembles on Surfaces 

A new microscopic technique [43] is demonstrated that combines attributes 
from both near-held scanning optical microscopy (NSOM) and fluorescence 
resonance energy transfer (FRET). The method relies on attaching the ac- 
ceptor dye of a FRET pair to the end of a near-held flber optic probe (tip). 
Light exiting the NSOM probe, which is non-resonant with the acceptor dye, 
excites the donor dye introduced into a sample. As the tip approaches the 
sample containing the donor dye, energy transfer from the excited donor to 
the tip-bound acceptor produces a redshifted fluorescence. 
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The main information about FRET can be taken from the fluorescence 
kinetics of the donor or acceptor. Hence in this case a combination of NSOM 
and correlated single photon counting methods is needed. 

This method is shown to provide enhanced depth sensitivity in fluores- 
cence measurements, which may be particularly informative in studies on 
thick specimens such as cells. The technique also provides a mechanism for 
obtaining high spatial resolution without the need for a metal coating around 
the NSOM probe and should work equally well with non- waveguide probes 
such as atomic force microscopy tips. This may lead to dramatically improved 
spatial resolution in fluorescence imaging. 
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4 Scanning Tunneling Spectroscopy 
and Electronic Properties 
of Single Fullerene Molecules 
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The discovery of an effective preparation procedure for obtaining macro- 
scopic quantities of Ceo, a third form of carbon besides graphite and dia- 
mond, has triggered a significant research effort to understand the physi- 
cal and chemical properties of Ceo and Ceo-related fullerene molecules [1]. 
Fullerene molecules are closed cage molecules containing only hexagonal and 
pentagonal faces. Chemically stable, cage-structured Ceo molecules form an 
interesting new family of adsorbates on surfaces. They differ significantly 
from the elemental or simple molecular adsorbates because of their three- 
dimensional character on the atomic scale. A unique fundamental property 
of this type of adsorbate is molecular orientation with respect to the host 
substrate. When an isolated molecule ceases its rotational motion on a sur- 
face, it may, in general, adopt a number of binding configurations and hence 
a range of different orientations. Interaction between molecules could yield 
still other orientational arrangements. Understanding site- or orientation- 
dependent Ceo'Substrate interactions therefore lies at the heart of the design 
of new catalysts with functionalized cage molecules or the fabrication of thin 
films with desirable orientational orders. However, most surface analysis tech- 
niques (e.g., HREELS, PES) can only provide statistically averaged results 
rather than spatially resolved results. 

On the other hand, device sizes will ultimately decrease down to the 
nanometer scale and the building blocks and functional units of nano-devices 
will be individual molecules and even individual atoms in the near future. 
Recent breakthroughs have been achieved through scientists’ skill in manip- 
ulating individual atoms and molecules with STM, and a variety of molecular 
devices and concepts like atom relays, single- atom switches and molecular am- 
plifiers have been proposed and realized [2-6] . Cage-structured Ceo molecules 
possess many interesting physical and chemical properties and have shown 
their amazing potential in applications for molecular devices. For example, 
the first single-molecule electromechanical amplifier with a gain of five has 
been produced by a metal tip pressing a Coo molecule [6]. This is based on 
the modulation of virtual resonance tunneling through the Coo molecule by 
electromechanical deformation of the Coo cage. The interplay between the 
single-electron tunneling effect and splitting of degenerate molecular levels 
has also been observed for Coo molecules deposited on an insulating layer at 
4.2 K [7]. 
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In this chapter, the method of scanning tunneling spectroscopy (STS) 
is introduced first and then recent STM/STS studies on Ceo molecules are 
reviewed. We will focus on the following two aspects: 

• local electronic structures of individual Ceo molecules on Si substrates, 

• electronic transport properties of individual Ceo molecules. 

4.1 Scanning Tunneling Spectroscopy 

Scanning tunneling microscopy (STM) was invented by G. Binnig and 
H. Rohrer [8-10]. The basic principle of the instrument is quantum mechani- 
cal tunneling of electrons between a metallic tip and the sample which can be 
a flat surface of any conducting or semiconducting material. The magnitude 
of the tunneling current is exponentially proportional to the inverse distance 
between the tip and the surface, so it is possible to obtain atomic resolution 
images by scanning the tip over the surface. If the tip is held over a given 
point on the surface of the sample, and the voltage is swept to make a series 
of current measurements as a function of voltage, an I-V curve can then be 
measured corresponding to that location. According to the theory of Tesoff 
and Hamann [11], the tunneling current in the STM can be expressed as 



I{V) oc 


rEp -f-e V 

/ p{r,E)dE, 


(4.1) 


p{r, E) = 


= Yl\Mr)\^SiE-Ei), 


(4.2) 



i 



where p{r,E), '0i(^), and are the local density of states (LDOS) of the 
sample, the sample wave function with energy Ei^ and the Fermi energy, 
respectively. The normalized derivative {dI/dV)/{I/V) of this curve, where 
/ is the tunneling current and V the sample bias voltage, reflects the LDOS 
of the samples in a first order approximation [12]. 

Consequently, scanning tunneling spectroscopy (STS) was expected to 
offer a novel view of the electronic structure of surfaces because it could 
probe the energy dependence of the sample electronic states with unprece- 
dented spatial resolution. Although early studies fulfilled this promise in the 
dangling-bond semiconductor [13], the progress of STS in metals has been 
slower than expected. The difficulties encountered arise on both the exper- 
imental and the theoretical sides. The results for the tunneling current as 
well as the spectroscopy depend on the geometric and electronic structures 
of the tip. A blunt tip with a featureless LDOS is required for reproducible 
spectroscopy reflecting the LDOS of the sample in the first approximation 
mentioned above [12]. But a blunt tip is easily ‘sharpened’ by accidentally 
attaching a foreign adatom or small cluster of atoms during the scanning. 
Although this is good for producing a high resolution image, such sharpening 
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Fig. 4.1. Calculated local density of states (s states) at the apex atom for bcc (line 
composed of circles) and fee (solid line) (lll)-oriented W pyramids. The arrow 
indicates the position of the Fermi level. The inset shows the two tip geometries 
considered [14] 



leads to a tip structure exhibiting resonant electronic states that strongly 
modify the spectra. On the other hand, even the clean tip may have resonant 
states around the Fermi energy if the tip ends with a specific structure. For 
example, Vazquez et al. [14] have shown both experimentally and theoret- 
ically that while the (111) bcc W-tip yields a smooth, bulklike density of 
states (DOS), the (111) fee tip exhibits two narrow resonances within the 
upper energy range of the bulk spectrum on both sides of (see Fig. 4.1). 
Despite the lack of complete experimental control over the tip structure and 
consequent difficulties of reproducibility in the practical context, STS is still 
a truly unique way to probe local electronic properties. 

Recently, applications of STM to molecules or superamolecular assem- 
blies have attracted considerable interest. The STS of individual adsorbed 
molecules is expected to provide valuable information about the electronic 
structure and quantum transport properties. However, this information de- 
pends on the molecule-substrate interaction and coupling of the molecule to 
its environment. This is therefore a challenging task and several key issues 
have to be addressed in addition to the problems for semiconductor and metal 
surfaces. Firstly, STS measurements for molecules are usually conducted on 
the molecule-substrate system, so suitable atomically flat and conducting 
substrates have to be found. Secondly, interpretation of the tunneling spec- 
trum of molecular species requires a deep understanding of their electronic 
structure and transport properties. If the interaction between molecule and 
substrate is strong, it will be more difficult to understand the peaks in the 
spectrum since the interaction gives rise to strong chemisorption bonds and 
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a corresponding alteration of the energy levels of the adsorbed molecules. 
Thirdly, it is important to eliminate the surface mobility of the molecules to 
allow stable STM imaging and spectrum measurement. This may be achieved 
in different ways, for instance, by physisorption at low temperatures or self- 
organization of molecular structures. 

In the following sections, we will use the Ceo molecule as the model system 
in order to show how to obtain information on molecular electronic struc- 
ture, site- and orientation-dependent molecule-substrate interactions and 
electronic transport properties from elaborate STS experiments combined 
with theoretical calculations. We will also discuss single- molecule devices us- 
ing Ceo as functional unit. 

4.2 Electronic Structures of Fullerene Molecules 

The structure model of the Ceo molecule was proposed by Kroto et al. in 1985 
[15]. It consists of twelve pentagonal rings and twenty hexagonal rings, and 
no pentagonal rings make direct contact each other (see Fig. 4.2). According 
to NMR experiments [16], all C atoms are identical in a Ceo molecule. The 
bond lengths are 1.46 A and 1.40 A, respectively, for the C-C single bonds 
between a pentagon and a hexagon and C=C double bonds between two 
hexagons. The bonding character of Ceo is predominantly with a small 
admixture of character due to the non-zero curvature. For this reason, 
the electronic states can be decomposed into approximate tt and a states. 
The bonding a states reside well below the highest occupied level, which is 
composed of orbitals having primarily tt character [17,18]. The nearly spheri- 
cal structure of the Ceo suggests a labelling of these electronic states in terms 
of spherical harmonics, with the a and tt electrons corresponding to different 
radial quantum numbers. Theoretical calculations [19] for the isolated Ceo 
molecule show the relative positioning in energy of these molecular orbitals 
(Fig. 4.3). The tt electrons pair up to fill 30 energy states. The degeneracy of 
the angular momentum L — h level, which contains 11 states, is removed by 
the icosahedral symmetry of the molecule, resulting in three separate levels: 




Fig. 4.2. Ball-and-stick (a) and bond model (b) of a Ceo molecule. Each bond on 
the pentagonal rings is a single C-C bond, while those shared by two neighboring 
hexagonal rings are double C=C bonds 
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Fig. 4.3. Density of states for the isolated Ceo molecule calculated using the discrete 
variational LDA method. Vertical lines denote energy levels of the molecule. The 
DOS is obtained by broadening the discrete levels with a Lorentz function for better 
comparison with experimental data 



• the ftu, or highest occupied molecular orbital (HOMO), containing 10 
electrons, 

• the triply degenerate tiu, or lowest unoccupied molecular orbital 
(LUMO), which can accommodate 6 electrons, 

• a similar triply degenerate t 2 n orbital (the LUMO +1) at higher energy. 

In order to study the electronic structure of an isolated Ceo molecule us- 
ing STS, it is necessary to put the molecules onto a substrate that has little 
or no interaction with Ceo- A Ceo monolayer supported on a self- assembled 
monolayer (SAM), as illustrated in Fig. 4.4a, is a good structure because the 
SAM is atomically flat and interacts very weakly with the Ceo molecules. The 
preparation of the SAM substrate takes two steps. A 160 nm thick Au(lll) 
film is first grown on a freshly cleaved and well-outgassed mica sheet by vac- 
uum deposition at 300° C. The mica-supported gold film is then quenched in 
ethanol and, while still covered with an ethanol droplet, immersed in a 2mM 
ethanol solution of alkylthiol (CH 3 (CH 2 )sSH) for about 48 hours at room 
temperature. This procedure routinely yields high quality SAM as verified 
by STM. A piece of the SAM substrate is cut and transferred into the ultra- 
high vacuum and low-temperature STM chamber where a submonolayer of 
Ceo is then deposited onto the substrate. The STM investigation shows that 
the Ceo molecules form close-packed hexagonal arrays on the SAM substrate 
(see Fig. 4.4b). 

At room temperature, the Ceo arrays are unstable. Molecules at the edge 
of an array can detach readily and diffuse to another part of the same array 
or other nearby arrays. This observation is truly indicative of a very weak 
interaction between the Ceo and the alkylthiol. The STS measurements were 
performed at 78 K using an UHV LT-STM with a Pt-Ir tip which was sub- 
jected to a careful cleaning treatment. Figure 4.5 displays a typical DOS that 
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(a) 





Fig. 4.4. (a) Schematic of a Ceo overlayer supported on the thiol SAM. (b) STM 
image of Ceo close-packed hexagonal arrays on the SAM substrate (40 x 40 nm^) 



was derived from I~V spectra measured by positioning the tip above an in- 
ner molecule in a selected island. The low temperature and the fact that Ceo 
molecules are assembled together both ensure that the molecules do not move 
around during tip scanning and spectrum measurement. Note that adjacent 
Ceo molecules do not modify each other’s electronic structures due to van 
der Waals intermolecular interactions, so STS data is expected to contain 
information relevant to an isolated molecule. Indeed, we find from Fig. 4.5 
that the main features of the experimental results are in good agreement 
with calculated theoretical results. 




Bias Voltage (V) 



Fig. 4.5. Density of states (A) for a Ceo molecule in a close-packed hexagonal 
island on alkylthiol (derived from the experimental I-V curve). Curve B is the 
theoretically calculated DOS for the isolated Ceo molecule 
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4.3 Scanning Tunneling Spectroscopy 

of Single Ceo Molecules Adsorbed on a Si Surface 

It is very important to study fullerene-surface interactions if we are to un- 
derstand the fundamental physical and chemical properties of fullerene and 
develop potential applications. The adsorption characteristics and molecule- 
substrate interaction of Ceo semiconductor surfaces as well as on metal 
surfaces have been investigated extensively by STM [20] and other techniques 
such as high-resolution electron-energy-loss spectroscopy (HREELS) [21,22] 
and photoelectron spectroscopy (PES) [23-25]. In this section we will discuss 
the scanning tunneling spectroscopy of single Ceo molecules on an Si surface, 
and focus on local electronic properties that depend on the adsorption site 
and molecular orientation. 

4.3.1 Ceo Molecules on an Si(lll)-(7x7) Surface 

The structure of the Si(lll)-(7 x 7) surface has been extensively studied 
and is now well established. According to the accredited model. Si atoms 
take several different bonding configurations within the 7x7 unit cell, and 
each corner hole atom, adatom and rest atom has one dangling bond for the 
Si(lll)-(7 X 7) surface. A strong interaction between the silicon surface and 
fullerenes could therefore be expected due to the high density of dangling 
bonds. Early STM studies of Ceo on Si(lll) indicated that there are four 
different adsorption sites (see Fig. 4.6) when Ceo coverage is much less than 
monolayer, and STM images showed similar strip-like intramolecular struc- 
tures with positive bias voltage [26]. Moreover, Li et al. [27] found that the 
height of Ceo is strongly dependent on the bias voltage. In principle, the pen- 
tagonal and hexagonal faces or edge atoms of the Ceo molecule could bond 
diflFerently with dangling bonds on surface Si atoms, and this results in vari- 
ous bias-dependent characteristics. Apart from showing that the interaction 




Fig. 4.6. Schematic of four possible Ceo adsorption sites on Si(lll)-(7x7): A 
faulted half, A' unfaulted half, B corner holes, C dimer lines. About 80% of the Ceo 
molecules are on sites A and A', 13% on site B, and 7% on site C 
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Fig. 4.7. A 25 X 25 nm^ STM topography image of 0.05 monolayer Ceo deposited 
on an Si(lll)-(7x7) surface 



between adsorbed Ceo molecules and substrate is strong, and that differences 
in bonding configuration (or orientation) exist for individual Ceo molecules 
on the different sites, these experiments have left the challenge of identifying 
the specific molecular orientation and its local electronic structure largely 
unanswered. 

Recently, we used low temperature STM/STS to identify the molecu- 
lar orientation of individual Ceo molecules and probed their local electronic 
structure on Si(lll)-(7 x 7) surfaces [28,29]. 

The experiments were carried out using LT-STM with an electrochemi- 
cally etched tungsten tip. All STM topographies and STS were obtained at 
78 K in order to reduce thermal noise. The tip-sample separations for STS 
measurements were controlled by setting the sample bias voltage at 2.0 V 
and setting the tunneling current in the range 0.1-0.2nA. Figure 4.7 shows a 
25 X 25 nm^ STM topography image of 0.05 monolayer Ceo deposited on an 
Si(lll)-(7 X 7) surface. 

Curves A, B and C in Fig. 4.8 are typical STS results corresponding to Ceo 
adsorbed at sites A, B, C, respectively, on an Si(lll)-(7 x 7) surface. Curve 
D is the average STS result on the bare Si surface. In order to discuss the 
energy band structure of the adsorbed Ceo, these STS results are represented 
in a graph of (d//dV)/(//V) vs. V. It is already known that the tip status 
is of great importance in STS investigations. In order to ensure that our STS 
results are reliable and reproducible, we use the Si(lll)-(7 x 7) surface as 
a reference. Before and after the STS measurements on a Coo molecule, we 
performed STS measurements on the Si surface. The results indicate that the 
STS results on the Si surface are reproducible and have the same features as 
reported previously [13]. 

We found that the band gaps of Coo adsorbed at sites A, B and C are 1.4, 
0.8 and 1.3 eV respectively. These values are much smaller than the value for 
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Fig. 4.8. STS results: curves A, B and C are the averaged STS results corresponding 
to Ceo adsorbed at sites A, B, and C, respectively, on an Si(lll)-(7 x 7) surface. 
Curve D is the averaged STS result on the Si(lll)-(7 x 7) surface 



free Ceo (1.9 eV, see Fig. 4.2 and [30,31]). Previous studies have shown that 
Ceo molecules adsorbed on an Si(100)-(2 x 1) surface have a band gap of 
about 1.8 eV [33,34]. The reduced band gap of Ceo adsorbed on the Si(lll)- 
(7x7) surface can be explained as follows. Because of the strong interaction 
between Ceo and dangling bonds of the Si(lll)-(7 x 7) surface, Ceo molecules 
are distorted by bonding with Si atoms and their symmetry is reduced. The 
splitting orbitals derived from the symmetry reduction will mix with the 
surface states of the Si(lll)-(7 x 7) surface [35,24]. In this way, the LDOS 
of the Si surface is incorporated into the band gap of the Ceo molecules. 
Therefore, strictly speaking, the STS results in Fig. 4.8 no longer originate 
from Ceo molecules, but from a 'compound’ made from a Ceo molecule and 
Si atoms. For convenience, we will still denote these curves as the STS of Ceo 
molecules in the following unless otherwise specified. 

Moreover, since the surface structures of Si(lll)-(7 x 7) differ from one 
another at different sites, so the mixing process of the Ceo orbitals with the 
LDOS of Si(lll)-(7 x 7) is also site-dependent. It is observed that the STS 
characteristics of sites A, B, and C are more similar at unoccupied states than 
at occupied states. This suggests that occupied orbitals of the Ceo molecule 
mix more strongly than its unoccupied orbitals with Si surface states. 

The HOMO and LUMO are the orbitals that dominate the physical and 
chemical properties of adsorbed Ceo molecules, so it is important to identify 
them in our STS results. Even though as discussed above the interaction 
between the Ceo molecule and the Si surface is strong, it is still reasonable 
to expect the original orbitals of free Ceo molecules to remain traceable. The 
reason for this is that Ceo molecules are rather large and stable. The cage 
structure of the Ceo molecule cannot be changed, so distortion of the Ceo 
molecules should be small. Moreover, the interaction between Ceo and Si 
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surface is localized at the bonding interface. Taking into account the fact 
that the original band gap between HOMO and LUMO for free Ceo is about 
1.6-1.9eV [30,31] and comparing STS curves with free Ceo energy levels, the 
HOMO and LUMO of Ceo molecules are assigned to the peaks of the STS as 
labeled in Fig. 4.8. 

To obtain a deeper understanding of the STS results and verify the accu- 
racy of the above discussion concerning electronic properties of adsorbed Ceo 
molecules, further theoretical calculations are necessary. However, we must 
first identify molecular orientations from the high resolution STM images 
before carrying out any calculation. It was found that the internal pattern of 
the Ceo depends strongly on the bias voltage and the tip-sample distance. 
Figure 4.9 shows the bias dependent results of two molecules adsorbed on an 
A-site (Figs. 4.9al, a2) and a B-site (Figs.4.9bl, b2), respectively. When the 
sample is positively biased, i.e., electrons tunnel from the tip to unoccupied 
molecular orbitals and empty states of the substrate, an intramolecular pat- 
tern with one bright pentagon ring plus two curved strokes is observed at the 
A-site. Similar reproducible results were also obtained for Ceo molecules on 
A'-sites. This is not surprising since the binding configurations are similar 
for these two sites (ignoring the faulted/unfaulted sublayers). At an A-site 
(or A'-site), the molecule appears as four bright stripes (Fig. 4.9al) when the 
sample is biased at —1.8V with respect to the tip, and the electrons tunnel 
from the occupied molecular orbitals and the filled states of the substrate to 
the tip. On site B, Fig. 4.9b2 shows that the internal pattern of the Ceo at 
positive sample bias is a bright pentagon ring plus three curved strokes to its 
right, which is significantly different from the pattern on site A (Fig.4.9a2). 
By contrast. Fig. 4.9bl shows that the negative bias STM image on site B 
is quite similar to that on site A (Fig. 4.9al), both exhibiting four slightly 
curved bright stripes. 

Reproducible high resolution STM images, such as Figs. 4.9al, 4.9a2, 
4.9bl and 4.9b2, showing clear intramolecular features of the Ceo molecules, 
raise the possibility of identifying the molecular orientations on Si(lll)-(7 x 7) 
surfaces. In general, however, the observed internal structures are not directly 
related to the internal atomic configurations, but rather to the electronic 
structure of the adsorbed fullerenes and surfaces. Mapping from the latter to 
the former clearly entails theoretical modeling. 

We adopted Tersoff and Hamann’s formula [11] and its extension in order 
to simulate STM images. Equations (4.1) and (4.2) assume a constant den- 
sity of states of the tip, which allows us to obtain STM images only from the 
LDOS of the sample surface. This assumption is tenable when the separation 
between tip and sample is big enough. The experimental conditions for ob- 
taining Figs. 4.9al, 4.9a2, 4.9bl and 4.9b2 meet this requirement. It is known 
that the LDOS near Ep makes a dominant contribution to STM images, so 
we can obtain the LDOS by adding up several orbitals around the highest 
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Fig. 4.9. High resolution STM images of individual Ceo molecules recorded with 
different bias voltages on an A-site of the Si(lll)-(7 x 7) surface: (al) Vs = —1.8 V, 
(a2) Vs = 2.5 V. On a B-site: (bl) Vs = —1.8 V, (b2) Vs = 2.3 V. The tunneling 
currents are 0.15 nA for positive bias images and 0.1 nA for negative bias images 



occupied molecular orbital and the lowest unoccupied molecular orbital with 
the same weight. 

We used cluster models to mimic an individual Ceo molecule adsorbed 
on Si(lll)-(7 X 7) with a number of binding configurations. In these cluster 
models, the Si dangling bonds irrelevant to Ceo adsorption are saturated by 
hydrogen atoms. We only considered two adsorption sites A and B. For Ceo 
adsorbed on site A, our cluster model is CeoSi 5 rH 42 , while it is CeoSierH 45 
for the B site. Figure 4.10 shows schematic views of the cluster models. 

The electronic structures of the clusters were calculated using the dis- 
crete variational LDA method, which has been described in detail elsewhere 
[36,37]. In formulating the Kohn-Sham equations, we used an exchange and 
correlation potential of the von Barth-Hedin form [38] with parameters taken 
from Moruzzi et al. [39]. The atomic basis functions representing the valence 
electron orbitals were 2s-~2p for C, 3s-3p for Si, and Is for H, and the rest of 
the core orbitals were treated as frozen. Using 600 sample points per atom 
in the numerical integration, we achieved sufficient convergence for both the 
electronic spectrum and the binding energy. A self-consistent charge model 
density was used to fit the electron density. We only optimized the height of 
Ceo above the Si surface to reduce computational efforts. Finally, the LDOS 
distribution on the upper hemisphere was calculated with a 16 a. u. radius 
(Rt-s) from the center of the Ceo- 

Our simulated STM images show that the positive bias images depend 
strongly on the orientation of the Ceo ori the substrate, and less on the 
adsorption sites. As examples. Figs. 4.11a and b show the calculated images 
with the sample biased at +2.5 V (Vg = 2.5 V) and the Ceo molecule with one 
of its 5-6 bonds adsorbed on A and B sites, respectively. Figures 4. lie and f 
are the calculated images with the same sample bias (Vg = 2.5 V) but the Ceo 
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Fig. 4.10. (a) Schematic view of the curved surface on which the tip moves around. 

(b) Cluster model of Si 5 ?H 42 for Ceo adsorbed on site A of an Si(lll)-(7x7) surface. 

(c) Cluster model of Sie 7 H 45 for Ceo adsorbed on site B 
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Fig. 4.11. Simulated STM images of Ceo on an Si(lll)-(7 x 7) surface, (a) Vs = 
2.5 V, 5-6 bond adsorption on A-site. (b) Vs = 2.5 V, 5-6 bond adsorption on B- 
site. (c) Vs = —1.8 V, 5-6 bond adsorption on A-site. (d) Vs = —1.8 V, edge-atom 
adsorption on B-site. (e) Vs = 2.5 V, edge-atom adsorption on A-site. (f) Vs = 2.5 V, 
edge-atom adsorption on B-site. Rts = 16a.u. for all cases 



molecules have the edge atom adsorbed on A and B sites, respectively. On the 
other hand, the simulated negative bias images are only weakly dependent 
on either the adsorption site or the orientation, and all display four bright 
stripes. Two calculated images with the sample biased at —1.8V on sites A 
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(5-6 bond adsorption) and B (edge-atom adsorption) are given in Figs. 4.11c 
and d, respectively. 

Comparing our simulations with experiment, we find that the experimetal 
image of Ceo on the A-site (Fig. 4.9a2) best matches with Fig. 4.11a (a bright 
pentagon ring plus two bright curved strokes to its left), while the image 
of Ceo on the B-site (Fig. 4.9b2) matches well with Fig. 4. Ilf (one bright 
pentagon ring plus three bright curved strokes to its right). We thus conclude 
that Ceo adsorbs on site A (or A') with one of its 5-6 bonds facing towards 
the surface, and on site B with one of its edge atoms facing towards the Si 
substrate. 

The ability to identify the orientation of Ceo molecules on specific adsorp- 
tion sites allows us to calculate the local electronic structures and theoretical 
STS results of oriented Ceo at those sites, and thus leads to a deeper under- 
standing of the Ceo~Si interaction. Figure 4.12a shows the energy states of a 
Ceo molecule adsorbed at an A-site with one of its 5-6 bonds facing down- 
wards. We noted that the highest occupied energies of a CeoSiesHse cluster 
and a free Ceo molecule are — 7.67eV and — 9.64eV, respectively. This means 
that electrons transfer from the Si substrate to Ceo- From the energy dia- 




Fig. 4.12. (a) Energy levels of a Ceo molecule adsorbed at an A-site with one of its 
5-6 bonds facing downwards, (b) Energy levels for Ceo molecule providing a major 
Mulliken population, (c) Energy levels for Si atoms providing a major Mulliken 
population, (d) Energy levels of a free Ceo molecule. For all cases, solid and broken 
lines denote occupied and unoccupied molecular orbitals, respectively, and units 
are eV 
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Fig. 4.13. Experimental STS and calculated LDOS of Ceo adsorbed at site B on an 
Si(lll)-(7 X 7) surface. The vertical lines at the bottom are the calculated energy 
levels of a free Ceo molecule 



gram, we also find that the valence and conduction bands are wide. This 
may be caused by the hybridization of C and Si atoms. 

The energy diagram for Ceo adsorbed on the B-site is as complicated as 
that for the A-site shown in Fig. 4.12a. We have calculated the LDOS of Ceo 
adsorbed on the B-site and the result is shown in Fig. 4.13 together with 
the experimental STS. We find that the theoretical result has reproduced 
the main features of the measured STS. We also plotted the energy levels 
of a free Ceo molecule calculated by the discrete-variational local density 
approximation (DV-LDA) method in Fig. 4.13. For every energy level near 
the Fermi energy of the free Ceo 5 we find a corresponding peak from the 
calculated LDOS, as well as from the experimental STS. This suggests that 
the basic characteristics of free Ceo molecules have not been changed after 
they adsorb on the Si(lll)-(7 x 7) surface. It also proves that it is feasible 
to determine the HOMO' and LUMO' of Ceo adsorbed on the Si surface. In 
the calculated LDOS, the mixing states are also found and are similar to the 
experimental STS result. The calculation gives the ionicity of the entire Ceo 
as about —0.5, indicating charge transfers from the Si surface to the adsorbed 
Ceo molecule, as confirmed by our STS results. Since it is believed that the 
bonding between Ceo ^md Si is of covalent type [24,40], so the charge transfers 
from the Si surface to the Ceo molecule should be implemented through a 
polar bonding mechanism. 

It is interesting to note that although site B is the one with the strongest 
interaction between Ceo and the Si surface, it is not the most commonly 
chosen adsorption site. Based on the above discussion, we suggest here a 
possible explanation for this. At site B, the two bonding sp^ hybridized C 
atoms are separated by a larger distance on the Ceo molecule, so the formation 
of these two bonds will need some extra energy to cause the originally nearest 
two sp^ hybridization orbitals, initiated by breaking a double bond, to relax 
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onto two separated C atoms in the Ceo molecule. This means that, even 
though the interaction between Ceo and the Si surface is strong at the B-site, 
it is difficult for a Ceo molecule to become chemically adsorbed there. 

4.3.2 Ceo Molecules on an Si(100)-(2x 1) Surface 

Hashizume et al. [32] first studied Ceo adsorption on the Si(100)-(2 x 1) 
surface. It was observed that almost all the Ceo molecules are absorbed in 
the trough between two dimer rows without any preference to the step edges 
or nucleation into islands at defect sites. Along the troughs, Ceo molecules are 
distributed randomly without forming islands. Indeed the interaction between 
Ceo and Si substrate is quite strong since Ceo molecules stay almost at the 
place they first land with little mobility. 

The STS of individual Ceo molecules on two different adsorption sites of 
the Si(lOO) surface was investigated by Yao et al. [33,34]. In Fig. 4.14 spec- 
trum A is the STS data measured on the molecule located at the top of 
a single dimer (site A in Fig. 4.15). The occupied state (HOMO) peak and 
the unoccupied state (LUMO) peak are located at — 0.7eV and — 2.3eV, re- 
spectively. The HOMO-LUMO energy gap is thus about 1.8 eV, which is 
about the same as that for free Ceo molecules (1.6-1.9eV). The wide peak 
at — 2.3eV (labeled I) corresponds to the two closely spaced hg and gg lev- 
els, which are separated from the HOMO by about 1.6 eV. The tig peak is 
found to be about 0.9 eV (labeled IV) above the LUMO level. STS data for 
Coo molecules on site A shows that they have roughly the same energy lev- 
els as free Coo molecules. This result implies that the breaking of one C-C 
bond does not have a major influence on the electronic structure of a type-A 
molecule. 

Spectrum B in Fig. 4.14 is the STS data measured for another type of 
Coo molecule which is usually located at a missing dimer defect (denoted as 




Fig. 4.14. STS measurements on Ceo molecules located at the top of a single dimer 
and a missing dimer defect of an Si(100)-(2 x 1) surface, respectively [33] 
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Fig. 4.15. Geometrical model for Geo adsorbed on an Si(100)-(2 x 1) surface. Open 
circles show possible adsorption sites 



a B molecule). We find a strongly occupied state located at — 1.5eV (labeled 
I) below and a wide band structure above (labeled II). The gap is 
about 2.2 eV, which is close to the bandgap of SiC [41,42]. Local atomic re- 
arrangement of the Si(lOO) surface, which includes a missing dimer defect 
accompanied by two C-shaped defects in neighboring dimer rows, can be 
found in the vicinity of every type B molecule. STS results indicate that the 
interaction between a type B molecule and the surface is much stronger than 
that for an A molecule. The result also suggests that B molecules form cova- 
lent bonds not only with second layer Si atoms, but also with the surrounding 
Si atoms in the first layer. 

Hashizume et al. [32] observed three or four bright strips for each Ceo in 
their high resolution negative bias STM images. The STM images reported 
by Yao et al. [33,34] provided a few non-strip patterns. Their results indicated 
that the intramolecular structure of type A molecules is very different from 
that of type B molecules and supported the STS conclusion that the local 
electronic structure of the A-site differs from that of the B-site. However, Yao 
et al. did not discuss Ceo molecular orientations at different sites. 

Several investigations have been done to try to understand the origin of 
these strip-like patterns. Kawazoe et al. [43] calculated the electronic struc- 
ture of the Ceo molecule on the Si(100)-(2 x 1) surface by means of a first- 
principles calculation. They took one layer of Ceo molecules, and the Si(lOO) 
substrate was treated as a positive charge background. According to their 
calculation, the electron charge density of distribution for each level can be 
mapped out and the charge density distribution of the 180th level (HOMO) 
and the 181st level (LUMO) exhibit strips similar to those in the STM im- 
ages. Yamaguchi [44] calculated the electronic structure of a model cluster of 
Ceo molecules on an Si(100)-(2 x 1) surface by the DV-X^ method. Yajima 
et al. [45] pointed out that it is the local density of states which should be 
compared with STM images, and that it is not sufficient to calculate the elec- 
tron density contour map of some particular energy states. Although their 
results reproduced the strip-like patterns in a way similar to the experimen- 
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34 Bohr 

Fig. 4.16. (a) Simulated STM image of the Ceo molecule on an Si(100)-(2 x 1) 
surface (sample bias -2.0 eV, tip height 12bohr) [45] and (b) the corresponding Ceo 
orientation 



tal results (see Fig. 4.16), it is still difficult to identify molecular orientations 
since the negative bias images are not dependent on molecular orientation 
and adsorption site, as shown above. Further experiments are needed to ob- 
tain high resolution positive bias images of the C6o/Si(100) system if we are 
to understand the local Ceo-Si interaction mechanism. 

4.4 Electronic Transport Properties 
of a Single Fullerene Molecule 

Until recently, the electronic properties of individual molecules had not been 
measured. The main difficulty is the question of how to connect individ- 
ual molecules to electrodes. The scanning probe microscope has dramati- 
cally changed this situation by introducing a molecule into tunnel junctions 
in which the tip and supporting surface are two electrodes. The STM tun- 
nel junction with a molecule at the interface has the advantage that single 
molecules can be observed on an atomically flat electrode (surface) before or 
after the electrical measurement. 

Molecules are typically characterized by their HOMO and LUMO states 
separated by an energy gap. The discrete levels of the HOMO and LUMO 
states are broadened and shifted through their interaction with the surface 
when the molecule is adsorbed on a substrate. When a voltage is applied 
across a molecule by two electrodes, electrons tunnel directly through the 
LUMO states and are therefore in resonance with these states if the bias 
voltage is larger than the energy gap. If the applied bias voltage is not high 
enough to allow electrons to tunnel directly through the molecular levels, 
electrons can still be transported by means of the tails of the broadened 
molecular levels. This transport mechanism, known as virtual resonance tun- 
neling (VRT), can be used to produce an electronic amplification phenomenon 
in single molecules. 

Recently, electronic properties of individual Ceo molecules have been in- 
vestigated using STM [7,46,47]. The resistance of a Ceo molecule was de- 
termined by defining the Tip contact’ with the molecule and carrying out 
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tunneling spectroscopy on an isolated Ceo in the presence of charging ef- 
fects. Single-molecule amplification by electromechanical modulation of vir- 
tual resonance quantum mechanical tunneling was also investigated, and it 
was shown that the STM tip can make, visualize and operate such a device 
rather than just function as an electrode. 

4.4.1 Electrical Resistance of a Single Ceo Molecule 
Determined by STM 

When the STM tip is very close to a conducting surface, the electrical con- 
ductance (resistance) of the metal-vacuum-metal junction is given by the 
slope of the I~V characteristic at low applied voltages where a linear rela- 
tionship is observed. Typically, the voltage (Vt) is less than 100 mV [48]. If a 
molecule is introduced into such a junction, the I-V characteristic will cer- 
tainly be changed, and in principle the electrical resistance of the molecule 
can be determined accordingly. However, some basic concepts regarding the 
electrical transport of a single molecule between two metallic electrodes have 
to be addressed. 

Buttiker and Landauer [49] have developed a concept of ‘electronic trans- 
parency’ to characterize the ability of a molecule to let electrons fiow through 
it. High transparency corresponds to a low resistance of the metal-molecule- 
metal junction. In contrast to macroscopic conductance determined by 
electron-phonon scattering, the transport of electrons through a molecule 
is determined by the quantum mechanical probability of one electron moving 
from one electrode to the other. As discussed above, if a molecule’s HOMO 
and LUMO states do not lie in the range of the Fermi level, tunneling occurs 
by a process known as VRT when the bias voltage is low enough. In this 
case, transparency is determined by the efficiency of coupling between the 
molecule and its metallic substrate or the non-zero weight of the tails of the 
broadened molecular levels at the electrode’s Fermi level. 

Consequently, the issue of the electronic conductance or transparency of 
a single molecule requires us to understand the molecule-substrate interac- 
tion and its inffuence on the broadening and shifting of molecular levels. In 
particular, when the STM tip approaches the molecule and the tip-molecule 
junction becomes smaller and smaller, molecular orbitals can be modified 
by mechanical deformation of the molecular conformation or shape, thereby 
also affecting the electronic coupling of molecule and electrodes. It is expected 
that such deformation and modification will significantly change the slope of 
the I-V characteristics, so that the meaning of electronic and mechanical 
contacts with an individual molecule should be defined in such a way as to 
get a reasonable value of single-molecule conductance (resistance) from the 
experimental results. 

Joachim et al. [46] did the first study of electrical contact with an individ- 
ual Ceo molecule. The experiment was performed using an STM at 300 K in 
ultrahigh vacuum. Submonolayer Ceo molecules were deposited by thermal 
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Fig. 4.17. Variation of the tunneling current intensity I{s) through a Ceo molecule 
and corresponding resistance as a function of the distance 5 from tip apex to surface. 
Open circles indicate experimental values. The solid line indicates calculated values 
with full deformation of the Ceo upon tip approach. In both cases, a bias voltage 
of 50 mV was chosen [46] 



evaporation onto an Au(lOO) surface. The Ceo HUMO-LUMO gap is not too 
wide in its adsorbed state on An (about leV) [50], so a reasonable current 
intensity can be measured at low bias voltage within the VRT regime. The 
tip was positioned above a Ceo and the electrical current I was studied as a 
function of tip displacement s towards the molecule [I{s) characteristic]. This 
was achieved by superimposing a series of offsets A^ and recording individ- 
ual I-V characteristics over a range of resistances varying from 5 MQ to 
R ^ 5 000 M^2. It was confirmed that, to within the bias voltages of ±200 mV, 
I-V characteristics recorded at fixed s are linear for a series of tip altitudes. 
This linearity allows us to define the tunneling resistance by I{s) measure- 
ments at constant Vt = 50 mV. Figure 4.17 shows a representative I{s) curve 
averaged over some 20 measurements at a bias voltage of 50 mV. All the I{s) 
curves display a linear log I variation as the tip is ramped forward, followed 
by a large and reversible positive deviation at smaller gap values. 

Joachim et al. [46] also calculated the theoretical I{s) curve using the 
elastic scattering quantum chemistry technique configured for STM (STM- 
ESQC) [51]. The full atomic structure and valence orbitals of the tip, the Ceo 
molecule, and the Au(llO) substrate were considered in these calculations. 
The Ceo deformation was included in STM-ESQC using the MM2 routine 
with a standard sp‘^ carbon parametrization [52], and the Ceo conformation 
was optimized for each s. The current I through the molecule was then cal- 
culated for each conformation. 
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Fig. 4.18. Optimized structure of the Ceo molecule in the tunneling junction for a 
tip apex to surface distance of (a) s = IdA and (b) s = 7.35 A. The W tip apex 
was considered to be rigid during approach [46] 



The deformation of Ceo by the W tip apex is presented in Fig. 4.18, as 
optimized during STM-ESQC calculations. The calculated I{s) curve (solid 
line in Fig. 4.17) agrees very well with the experimental results when junc- 
tion distances are between 16 and 14 A. As s decreases down to 13 A, a 
small deformation of the Ceo cage occurs due to van der Waals tip attrac- 
tion and the deformation leads to a deviation from the linearity of log/(s). 
In the range 16-13 A, the junction is in a tunneling regime. According to 
the optimized energy of Ceo iii the junction, mechanical contact is estab- 
lished at s ~ 13.2 A, corresponding to zero force on the Ceo in the junction 
[d(£'(Ceo)/d5 = 0]. At this point the Ceo is slightly attracted by the tip. At 
s = 12.3 A, the van der Waals expansion is compensated by tip compression 
and the Ceo recovers approximately its original shape. If this point is defined 
as the point of electrical contact, the resistance of the total junction with 
Ceo is 54.807 ± 13 MQ as determined from experimental measurement (see 
Fig. 4.17). When s < 12.0 A, the cage undergoes compression and R{s) de- 
creases rapidly towards the quantum resistance limit of Rq = hj2e^. At this 
limit the molecule is almost totally transparent to the tunneling electrons. 

Is the resistance of 54.807 MQ the intrinsic electrical resistance of an indi- 
vidual Ceo molecule? At low bias voltage, all the I-V characteristics remain 
linear within the HOMO-LUMO gap and the tunneling current through Ceo 
is dependent on the symmetry and spatial extension of its molecular orbits 
[51]. In this sense, the resistance may be considered to contain an intrinsic 
characteristic of the molecule. Nevertheless, the current flows because the Ceo 
molecular orbitals are broadened and shifted by their electronic coupling with 
the tip and surface. Joachim’s observation of a pressure dependence of the 
electronic transparence of an individual Ceo is a step toward a realistic char- 
acterization of the electrical properties of a single molecule in a nanoscopic 
environment. Moreover, this study indicates key differences in the physics of 
Ceo confined in such a situation compared to bulk solid Ceo [53,54]. 
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4.4.2 Resonant Tunneling of Single Ceo Molecules 
in a Double Barrier Tunneling Junction (DBTJ) 

Metal clusters with a size of a few nanometers have been intensively studied 
[55-66]. When the metal cluster is placed as the center electrode of a DBTJ, 
the current-voltage (I~V) curves of the structure show Coulomb blockade 
and Coulomb staircase behavior due to the discreteness of the charge and 
the significant charging energy Ec of the cluster as compared to the thermal 
energy fcpT [55-63,66]. When the size of the metal cluster is small enough, 
the quantum size effect further quantizes the energy levels which can be ob- 
served in the I~V characteristic of a DBTJ as well [64-66]. The interplay 
between the single-electron tunneling (SET) effect and the quantum size ef- 
fect can be experimentally observed most clearly when the charging energy 
Ec of the cluster by a single electron is comparable to the electronic energy 
level separation JEp, and both energy scales are larger than fepT. Both ef- 
fects are relevant to the development of nanoscale electronic devices, such 
as single-electron (and single-molecule) transistors, and have thus attracted 
considerable experimental and theoretical attention [67-69]. 

An isolated Ceo molecule provides an ideal model system for this re- 
search. It is unique due to its size (~ 8 A in diameter), much smaller than the 
nanoparticles investigated so far, its spherical shape, and its rich electronic 
spectrum, which is affected by charging. Moreover, this system enables a si- 
multaneous study of the interplay between two regimes, one where the level 
spacing AEp > Eq (due to molecular orbital spacing), and the other where 
AEl < Ec (due to level splitting). 

Porath et al. studied the tunneling spectroscopy of isolated Ceo molecules 
in the presence of charging effects using a cryogenic STM [7]. They realized 
a DBTJ configuration in which a Ceo molecule is coupled via two tunnel 
junctions to the gold substrate and the tip of the STM (see Fig. 4.19a), and a 
thin insulating polymethyl methacrylate (PMMA) layer (less than 40 A thick) 
was used as the barrier between Ceo and gold substrate. The Ceo molecules 
were studied spectroscopically by taking tunneling I-V characteristics. Fig- 
ure 4.19b plots four I~V characteristics for a single molecule at 4.2 K with 
different STM settings (shifted for clarity). The main features in the I-V 
curves are as follows. Firstly, there is a nonvanishing gap in the curves around 
zero bias that oscillates periodically between minimum and maximum widths 
as we monotonically change A for fixed (and vice versa). This oscillation 
is due to the periodic variation of ‘fractional charge’ (Qo), consistent with 
the ‘orthodox’ theory [70,71]. The fact that the minimum gap is nonzero 
indicates an additional contribution related to the molecular level spacing 
besides the ‘classical’ Coulomb blockade. Moreover, in contrast to pure SET 
characteristics, the curves become highly asymmetric and exhibit steps with 
variable widths and heights. 

Figure 4.20a plots three experimental I-V characteristics (thick curves) 
taken with different junction parameters such that Ci < C 2 , Ci ^ C 2 and 
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Fig. 4.19. (a) Schematic and equivalent circuit of the experimental setup: a Ceo 
molecule [quantum dot (QD)] coupled via two tunnel junctions to the STM tip and 
gold substrate, (b) Tunneling I-V characteristics at 4.2 K taken for a Ceo molecule 
in a DBTJ configuration, with different STM settings. Inset: an I-V trace taken 
near a Ceo molecule showing Ohmic behavior [7] 
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Fig. 4.20. Tunneling I-V characteristics and their d//dF traces at 4.2 K [thick 
curves in (a) and (b)] and corresponding theoretical fitting [thin curves in (a) and 
(b)] for Ceo in a DBTJ. Curves group I: Ci < C 2 . Curves group II: C\ ^ C2. Curves 
group III: Cl > C2 [7] 

Cl > C 2 5 respectively, where Ci is the capacitance of the molecule-electrode 
junction and C 2 the capacitance of the tip-molecule junction. According to 
the ‘orthodox’ theory [70,71], the step width from the SET effect is propor- 
tional to the reciprocal of the larger of C\ and C 2 . From Fig. 4.20, we find 
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that these curves contain additional steps beyond those related to the SET 
effect, which reflect the discrete electronic structure of the Ceo molecule. If 
the data is plotted in the dl/dV traces, shown in Fig. 4.20b, the structure 
is exhibited in a more pronounced way in which each peak corresponds to 
resonant tunneling through a discrete level. 

Porath et al. also did theoretical calculations to explain the experimental 
spectrum. They used Mrthodox’ theory [70] but modified the formalism to 
include the electronic spectrum of Ceo- Information on both the junction 
parameters and the electronic molecular spectrum was obtained from the 
calculation. For an unperturbed, neutral Ceo molecule, the spherical L = 5 
levels are split into three groups; fivefold degenerate hiu levels (HOMO), 
threefold degenerate tiu (LUMO), and threefold degenerate t 2 u (LUMOd-2). 
The LUMO+1 (tig, threefold degenerate) originates from the L — 6 levels. 
These degeneracies can be lifted by two different means that break icosahedral 
symmetry: the electric field between the tip and the substrate, and the Jahn- 
Teller (JT) effect [72,73], which is prominent in the ionized and excited states 
of the molecule. Accordingly, they fitted all the curves using three groups 
of levels related to the HOMO, LUMO, and the LUMO+1 orbitals. The 
results are also represented in Fig. 4.20. One can observe a full splitting of 
the five HOMO levels and the three LUMO levels as well as a partially lifted 
degeneracy in these levels. The fittings were performed by taking the HOMO 
levels in the range —0.8 to —0.6 eV, the LUMO in the range 0-0.15 eV, and the 
LUMO+1 onset at ~ 0.6 eV. It should be noted that the observed HOMO- 
LUMO level spacing is ~ 0.7 eV, much smaller than the level spacing in the 
free neutral molecule. This result is reasonable because they are actually 
measuring the gap between the LUMO (C^q) and the HOMO (C^q), and 
because the threshold for tunneling in the DBTJ is determined by the energy 
of the final, ionized state of the molecule. The typical level splitting they 
obtained, ~ 0.05 eV, is somewhat smaller than predictions for the JT effect 
[72,73]. 

The effect of the junction capacitances on the shapes of the I-V curves 
is also an interesting question. The calculated results indicate that both the 
charging energy and the level separation exceed the thermal energy ksT even 
at room temperature, so the SET effect and the interplay with the molecular 
levels can be clearly observed [74]. When Ci < C 2 , tunneling through the 
DBTJ commences at junction 1, for the level configuration at hand, where 
electrons tunnel into the LUMO for positive bias polarity, and off the HOMO 
for negative bias. This is a major source for the pronounced asymmetry ob- 
served in the curves. The above interpretation is supported by the fact that 
in the opposite case, when C\ > C 2 , the I~V curves obtained are nearly an 
inversion (with respect to U = 0) of those in Fig. 4.20. 

In summary, the discrete tunneling spectra of isolated Ceo fullerenes in 
the DBTJ configuration manifest the interplay between charging and quan- 
tum size effects simultaneously in two regimes: AEi^ > Eq and AjFl < Ec- 
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The degenerate HOMO, LUMO, and possibly the LUMO+1 levels of the 
unperturbed molecule were split and fully resolved spectroscopically. The de- 
gree of splitting is somewhat smaller than predictions for the JT effect. Local 
fields may contribute to this degeneracy-lifting as well. Theoretical curves, 
calculated using the ‘orthodox’ model for SET, modified to account for the 
discrete levels spectrum of Ceo, agree well with the experimental data. The 
work of Porath et al. indicates that, in spite of the small size of the molecule, 
one can associate with it an effective capacitance similar to that of metallic 
particles of comparable size. A quantum-mechanical approach is probably 
more suitable for treating the capacitance in this molecular size regime. 

4.4.3 Molecular Device 

Using a Single Ceo Molecule as the Functional Unit 

Design and fabrication of molecular devices based on their amazing functional 
quantum properties using STM has been an inspiring field recently [75,76]. 
Cage-structured Ceo molecules have many interesting properties and have 
shown their amazing potential in applications to molecular devices. 

Joachim and Gimzeweski have proposed and demonstrated a single- 
molecule amplifier [77]. It is based on an extension of their measurements 
discussed in Sect. 4.4.1 as regards the electromechanical modification of VRT 
for a single Ceo molecule contacted by two electrodes. The principle of the 
device can be described in more detail as follows. When the tip apex ap- 
proaches the top of the Ceo cage, electrical contact between the molecule and 
the electrode is defined as the point when the tunneling current intensity is 
maximum through Ceo while maintaining the condition of minimum defor- 
mation of the Ceo cage. Such a condition is realized by selecting Vg (the gate 
voltage of the amplifier) in such a way that the HOMO-LUMO gap is almost 
that of an unperturbed Ceo molecule. When the tip is further lowered by 
0.1 nm toward the molecule, cage deformation leads to an enhancement of 
current intensity by two orders of magnitude. Within this 0.1 nm range, the 
I-V characteristic of the device remains linear if Up (the voltage through the 
polarization resistance) is kept below 200 mV. By optimizing the load and po- 
larization resistances (J?l and i?p) according to the Ceo resistance, under and 
out of compression, a gain of five has been obtained experimentally in such 
a molecular amplifier. This work demonstrates a new mechanical means for 
controlling quantum mechanical processes on the subnanometer scale. More- 
over, such a device has the amazingly tiny operation energy of one atto joule 
(10~^^ J) and this value is four orders of magnitude smaller than the best 
known solid state switch [78]. 

The Esaki tunneling diode [79], which has a wide variety of applications in 
the semiconductor industry, such as fast switches, oscillators, and frequency- 
locking circuits, has also been exploited at nanometer scales. The negative 
differential resistance (NDR) characterized by the phenomenon of decreasing 
current with increasing voltage in the I~V curves is the essential property of 
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the Esaki diode. It has been shown experimentally [80-85] and theoretically 
[85,86] that NDR could arise in the tunneling structures with an STM tip and 
a localized surface site when the LDOS of both the tip and the surface site 
have sufficiently narrow features. The sharp LDOS features of the surface site 
have been realized by, e.g., doping B atoms on the Si(lll) surface, [80-82], 
depositing Ag atoms on a semiconductor surface [83] , using rebonded B-type 
steps on the Si(001)-2 x 1 surface [84], or synthesizing a self-assembled mono- 
layer of small organic molecules on a gold substrate [85]. The narrow LDOS 
features of the tip are usually realized by ‘sharpening’ the tip apex, which 
means attaching a foreign adatom or small cluster. We note that the structure 
of such a ‘sharpened’ tip and the associated LDOS was uncontrollable and 
unstable to a certain extent in previous experiments [80,81,83]. Thus NDR is 
not always observable over a particular type of site, but depends on the tip 
status which may vary even within a single experiment [80,81]. 

We propose here an NDR molecular device involving a Pt-Ir tip with a 
Ceo molecule on its apex and another Ceo molecule physically adsorbed on 
a hexanethiol SAM. Figure 4.21a shows the schematic diagram of the NDR 
device. Previous experiments showed that a Ceo molecule could easily be 
attached to the tip when a Pt-Ir tip was brought near the sample surface 
and rastered through Ceo islands [87]. Because of the enhanced metal-Ceo 
interactions [87,88], the attached Ceo molecule can be stable at the apex 
of the tip. It is well known that a free Ceo molecule has /h symmetry and 
narrow- featured electronic structure. The interaction between the adsorbed 
Ceo molecule and the Pt-Ir tip is not strong enough to destroy the molecular 
symmetry [88,89], so the narrow features of the electronic states of Ceo are 
preserved. As for the Ceo molecules adsorbed on the hexanethiols, they are 
weakly bound to the substrate by van der Waals forces. Consequently, the 
sharp features of the LDOS of the Ceo on the hexanethiol surface can also be 
preserved. We therefore suggest that the NDR phenomena can be observed 
due to the narrow structure in the LDOS of the Ceo molecules on both 
electrodes. 

In order to confirm that NDR could be realized from the tunneling be- 
tween two Ceo molecules, we performed a theoretical simulation for the tun- 
neling structure shown in Fig. 4.21a. As a first order approximation, the tun- 
neling current can be expressed as [90] 

f‘E-p-\-eV 

I{V)(x / pt{E - V)ps{E)exp [ - 2s^/2{W - E) + V]dE , (4.3) 

where Pt{E) and ps{E) are the LDOS associated with the tip and sample 
respectively, W is the height of the tunneling barrier, and s is the tip-sample 
distance. 

The tunneling structure in Fig. 4.21a consists of two junctions: junction 
1 is between Ceo on the tip and Ceo on the hexanethiols, and junction 2 
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Fig. 4.21. (a) Schematic diagram of the NDR device involving one Ceo molecule 
physically adsorbed on the SAM and another adsorbed on an STM tip. (b) The 
simulated LDOS {upper part) of Ceo on the tip and Ceo on the SAM, and corre- 
sponding calculated I-V curve {lower part) 



between Ceo on the hexanethiols and the Au(lll) substrate. The voltage 
drop across junction 1 is 






VC2 
Cl + C2 



(4.4) 



where V is the applied bias voltage, and C\ and C 2 are the electric capaci- 
tances of junction 1 and junction 2, respectively. Considering that the sepa- 
ration between the centers of Ceo on the tip and Ceo on the hexanethiols is 
assumed to be 10.5 A, the thickness of the hexanethiol SAM is about 10.2 A 
[91], the dielectric constant for SAM on Au(lll) is taken as 2.6 [92], and 
Vi/V is estimated to be 0.72. In the simulation, the energy levels of a free 
Ceo molecule were calculated by the DV-LDA method. Then the LDOS of the 
Ceo molecules on both electrodes were derived from the Lorentz broadening of 
these discrete energy levels [93]. The Lorentzian width depends on coupling 
with the environment, as well as on the temperature. Figure 4.21b shows 
the LDOS of the Ceo molecules obtained by assuming a constant Lorentzian 
width (0.18 eV) for all the energy levels. The shift of the LDOS of the Ceo 
on the tip is due to charge transfer from the tip to the Ceo molecule [87,88]. 
In the present work, we assume a 0.15eV shift of the Fermi level toward the 
unoccupied states. The simulated I-V curve clearly shows the NDR effect at 
both positive and negative bias. 
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From the simulation, we can conclude that an NDR molecular device can 
be realized by the tunneling structure involving a Ceo adsorbed on a metal tip 
and another Ceo weakly bound to the hexanethiol SAM. This is attributed 
to the narrow-featured LDOS near the Fermi level of both the Ceo molecules. 
Compared with the previous atomic scale NDR device, the advantage of the 
proposed Ceo molecular NDR device is that it has a controllable and relatively 
stable tip structure and therefore stable and reproducible NDR is expected 
in such a device. 

4.5 Conclusion 

We have reviewed progress in the study of single fullerene molecules using 
STS. From the STS data we can determine site- and orientation-dependent 
electronic structures of single molecules adsorbed on a surface. In princi- 
ple, the method described in this section can also be applied to other kinds 
of molecule-substrate system. However, we should interpret STS data with 
caution since the electronic structure of the STM tip is to some extent un- 
controllable and so may strongly modify the measured spectrum. Elaborate 
spectroscopic experiments should therefore be combined with high resolution 
STM images and theoretical modeling in order to eliminate the tip effect and 
thus extract intrinsic information about the adsorbed single molecule. More- 
over, experiments on individual molecules using STM have permitted deeper 
insight into the quantum electronics of molecular systems such as the inter- 
play between the single-electron tunneling effect and splitting of degenerate 
molecular levels. At the level of an individual molecule, unique information 
has been obtained about conformational and mechanical properties and their 
relations to quantum-electronical behavior. All this has extended our knowl- 
edge of single molecules and may lead to the design and construction of new 
artificial molecular devices and machines. 
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5 Interaction 

Between Nonrelativistic Electrons 
and Optical Evanescent Waves 



J. Bae, R. Ishikawa, and K. Mizuno 

Scanning near-field optical microscopy (SNOM) has attracted much atten- 
tion as a new optical technique to circumvent diffraction limits of the wave- 
length in conventional optics [1,2]. Near-field microscopes [3,4] have demon- 
strated subwavelength resolution of several tens of nanometers for imaging 
various kinds of objects, including semiconductors, superconductors, and bio- 
logical tissue. In SNOM, near-field images of objects are produced by scatter- 
ing the optical near-field on objects with a small probe tip such as a tapered 
optical fiber and a sharpened metal whisker. The following question always 
arises from the SNOM operating principle: how do SNOM images reflect op- 
tical properties of an object? The probe tip must be small compared to the 
object, not only to achieve higher resolution, but also to avoid strong distur- 
bance of the near-field distribution. From this point of view, free electrons in 
vacuum are most suitable as a near-field probe, because electrons have the 
smallest size and therefore the least influence on the field. 

Nonrelativistic electron beams in vacuum have been widely utilized for 
amplifying and oscillating coherent electromagnetic waves and for operat- 
ing electron microscopes with subatomic resolution. Operating frequencies in 
these beam devices lie mostly in the microwave or millimeter wave region [5] . 
Much effort has been made to increase the operating speed of devices using 
the (inverse) Smith-Purcell effect [6,7], but the frequencies still remain in 
the submillimeter wave region [8,9]. Consequently, the interaction behavior 
between electrons and light, including quantum effects that occur at optical 
frequencies, is still under investigation. 

In this chapter, we discuss the energy modulation of nonrelativistic elec- 
trons with evanescent waves contained in optical near-fields. In the next sec- 
tion, we present microgap interaction circuits that can be used at optical 
frequencies. In Sect. 5.2, a metal microslit among the microgap interaction 
circuits is adopted for interaction with an optical near-field and is discussed 
theoretically. Section 5.3 describes experiments performed at infrared wave- 
lengths. The results demonstrate that the infrared laser can modulate the 
electron energy for a low energy electron beam using the microslit circuit. 
The extension of infrared investigations to the visible light region is discussed 
in Sect. 5.4. Section 5.5 reviews results presented in this chapter. 
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5.1 Microgap Interaction Circuits 



5.1.1 Energy and Momentum Conservation 



In the interaction between electrons and light, energy and momentum con- 
servation must be satisfied [?]. Consider the case where an electron absorbs 
a photon. The dispersion relations for the electron and the photon are 



We — \Jnn?c^ c^pe^ 5 (5.1) 

= h,uj = cpp , (5.2) 



where W and p indicate the energy and the momentum, respectively, m is the 
electron rest mass, c is the speed of light, h = /i/ 27 t, h is Planck’s constant, 
and uj is the angular frequency. When We increases by Wp, the electron 
momentum Pe increases by 

A 

Ap=— 5.3 

V 

where v is the electron velocity. Since Pp = hu/c is always smaller than 
Zip, momentum conservation is not satisfied and consequently the interac- 
tion between an electron and a photon never occurs in free space. For the 
interaction, some circuits that supply the shortfall in the electron momentum 
are indispensable. 

Figure 5.1 shows possible interaction circuits having a smaller gap d than 
the wavelength: 



(a) a thin dielectric film, 

(b) a metal film gap, 

(c) a metal microslit. 



The dielectric film is the interaction circuit used by Schwarz and Hora in 
1969 [11]. The metal film gap is very similar to conventional circuits used in 
klystrons [5] . The metal microslit was proposed for the study of electron-light 
interactions in an optical near-field [12]. 

In Fig. 5.1, the effective interaction length in the three circuits is limited 
so that an uncertainty exists in the momentum. This additional momentum 
Pc from the circuits is given approximately by 



Pc ~ 



h 

d 



(5.4) 



If Pc > Ap^ momentum conservation can be satisfied for the interaction. In 
the case of the metal film gap, the condition for the gap d is derived from 
(5.3) and (5.4) as 



d< f3X (5.5) 

where (3 — v/c and A is the wavelength of light. For the experimental pa- 
rameters, A =488 nm and /3 =0.5 which corresponds to an electron energy of 
about 80 keV, the gap must be narrower than 244 nm. 
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(c) Metal Microslit 

Fig. 5.1. Schematic diagram of three different gap circuits. The gap width d is 
smaller than the wavelength of the incident laser 




The same considerations can be applied in a different way to the metal 
microslit shown in Fig. 5.1c, using an evanescent wave theory [12] rather than 
the uncertainty principle. A number of wave components with different wave 
numbers key are induced in the slit by the incident light. On the basis of 
the evanescent wave theory, the wave numbers key in the ^-direction could 
extend from zero to infinity. From momentum conservation, key must be equal 
to Ap/h. Using (5.3), the following relationship between key and v is found: 

/Cev = - • (5-6) 

V 

Since key is larger than the wave number ko = lo/c of light in free space, this 
wave is evanescent. This equation also represents phase matching between 
the electron and the optical evanescent wave, that is, Vp = n, where Vp is the 
phase velocity of the evanescent wave given by Vp = uu/key. The validity of 
these theoretical considerations will be discussed in more detail later. 

This type of metal microslit is more suitable than the other gap circuits 
for precise measurement of the electron-light interaction and for investigating 
quantum effects. This is because there are no disturbances caused by electron 
scattering in the metal film or the dielectric medium. 
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5.1.2 Transition Rates of Electrons 



The transition rates of electrons in a metal film gap and a thin dielectric 
film have been estimated theoretically [12]. The calculational model is shown 
in Fig. 5.2. The electron beam passes through the gap and travels in the 
positive x-direction, whilst the laser light wave is propagated in the positive 
^-direction. To simplify the calculation, the following assumptions were made: 

• the initial velocities of the electrons are same, 

• propagation modes of the light wave in the metal film gap and the thin 
dielectric film are the TEMqo and TMqi modes, respectively, 

• the incident light wave is polarized in the x-direction, 

• the gap materials have no rf loss, 

• the electrons interact only with the light. 



The second assumption is valid when the gap width is smaller than A/2 and 
A/(2\/n — 1) for the metal and dielectric film gaps, respectively, where n is 
the refractive index of the dielectric film. 

The transition rate for the electron energy in the metal film gap was 
calculated in accordance with analyses by Marcuse [14]. The calculated rate 
rcm is expressed as 






2qcf3^ 



iPi sin^ 



ujd 



(5.7) 



where q is the electron charge, i the electron current density, Eq the dielectric 
constant of free space, and P[ the power density of the incident light. The 
value of Wm represents the rate at which an electron absorbs one photon of 
energy huu. The rate for emission of the photon is almost the same. 

Similarly, it is possible to calculate a transition rate tCd for the dielectric 
film by quantizing the laser field of the fundamental TMqi propagation mode. 
The expression giving is 




Fig. 5.2. Gap circuit configuration used for theoretical analysis 
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(5.8) 



where 



f^ix + K = {nikof , kl^ + kl^ {riekof , 



\kex\=(- 

\rii 



k'lx tan k‘ixd/2 , 



and ri[ and rie are the refractive indexes inside and outside the dielectric film, 
respectively, ky is the wave number in the ^-direction in free space, and k[x 
and kex are the wave numbers in the x-direction inside and outside the film. 

From (5.7) and (5.8), the transition rates can be estimated as a function 
of the gap width. Figure 5.3 shows the calculated results used to find the 
optimum gap width. The parameters used in the calculation are (3 = 0.5, 
A = 780 nm, n[ = 1.45 (Si02), and ne = 1. As can be seen from (5.6), Wm 
changes sinusoidally and peaks at the gap width satisfying the equation 




Fig. 5.3. Calculated transition rates for an electron with velocity /3 = 0.5 as a 
function of the gap width of the metal film gap and the Si02 dielectric film. The 
gap width is normalized to the laser wavelength 
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d = PX 0 ^ (5-9) 

where m is an integer. Substituting /3 = 0.5 and m = 0 for the first peak, an 
optimum gap width of A/4 is obtained. The variation of in Fig. 5.3 differs 
from the predictions of the klystron theory in which the maximum value of 
Wm occurs at d ~ 0 [5]. The difference comes from the different treatments 
of the photon density in the gap. The klystron theory assumes that the total 
number of photons stored in the gap is constant, whereas in our treatment 
it is assumed that the photon density in the gap is constant and determined 
by the incident laser power. 

From Fig. 5.3, it can be seen that for the Si02 film, the first peak value 
of Wd is 0.18 times that of w^. The optimum film thickness of 0.38 is also 
longer than the optimum width of the metal film gap. In the dielectric film, 
the laser field is distributed outside the film as an evanescent wave so that 
the number of photons inside the gap is smaller than that for the metal gap. 
The longer gap width increases the number of photons inside the dielectric 
film. 

Figure 5.4 shows the calculated transition rates as a function of the light 
intensity for the metal film gap and the Si02 film. These transition rates 
represent the probability per unit time of one electron absorbing a photon. In 
the calculation, the optimum gap widths of 0.25 A and 0.38 A were used for the 
metal film gap and the Si02 film, respectively. Other calculation parameters 
take the same values as described previously. From Fig. 5.4, it can be seen 
that the transition rates are 1.1 x 10“^ and 2xl0“^/sfora power intensity of 
10® W/cm^ in the metal film gap and the Si02 film, respectively. The power 




Fig. 5.4. Calculated transition rates as a function of the power intensity of the 
incident light for the metal film gap and the Si 02 dielectric film at a wavelength of 
780 nm 
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intensity of the laser corresponds to an output power of 790 mW focused onto 
an area of diameter 10 fim. 

The transition rate just in front of the metal slit is intended to be the 
same as that of the metal film gap because the field distribution at the surface 
of the metal slit is similar to that in the film gap. A detailed analysis for the 
metal slit is described in a classical manner in the next section. 



5.2 Metal Microslit 

The metal microslit is used for generating optical near-fields in which an 
electron beam is modulated with light at optical frequencies. Optical near- 
held distributions on the slit were determined using the method of moments. 
The theory was checked by comparing with the measured held distributions 
at a microwave frequency. The energy changes of electrons in the near-held 
were evaluated numerically through computer simulations using the Lorentz 
force equations. From the results, the relationship between wavelength, slit 
width, and electron velocity was found. 



5.2.1 Near-Field Distributions 

The calculational model for a metal slit is shown in Fig. 5.5. In accordance 
with Chou and Adams’ analysis [15], the electric held components and 
Ey and a magnetic held component Hz in scattered waves from the slit were 
determined for an incident wave with held vectors E[ and H[. In the calcu- 
lation, the following assumptions were made: 

• a metal slit consists of two semi-inhnite plane screens with perfect con- 
ductance and zero thickness, 

• the incident wave is a plane wave with x-axis polarization and normal 
incidence upon the slit. 




Fig. 5.5. Calculational model for the metal slit circuit 
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Fig. 5.6. Experimental setup for measuring the near-field distribution on the metal 
slit at 9 GHz 



In order to check the theory, a near-field distribution on a metal slit was 
measured using a macrowave-scale model at 9 GHz (wavelength A = 33 mm). 
The experimental setup is shown in Fig. 5.6. The metal slit consists of two 
aluminum plates with a height of 400 mm, a width of 190 mm, and a thickness 
of 1mm. The antenna probe with a length of 1.6 mm was placed at the end 
of a thin coaxial cable with a diameter of 0.8 mm connected to a spectrum 
analyzer. This antenna detects an electric field intensity Ex in the x-direction, 
the dominant field in the interaction with electrons. 

Figures 5.7a and b compare the calculated and measured field intensity 
distributions of Ex on the slit. The slit width d was 0.5A. The field intensities 
were normalized to E[x which is the measured field Sit x = y — 0 without the 
slit. In Fig. 5.7, the theoretical and experimental field distributions are quite 
similar except for small ripples appearing at x/A ~ -0.5. The ripples could 
be caused by interference of waves scattered from the coaxial probe and the 
slit. 

The theoretical and measured variations oi Ex dit x — 0 are plotted in 
Fig. 5.8. The theory for d — 0.75 A agrees well with the measurement. When 
d decreases to 0.1 2 A, the deviation from the measurement increases due to 
the probe having a finite size of about 0.05A. However, the theory accurately 
predicts the measured variations of the field intensities. These results indicate 
that the theory is valid within the bounds of experimental error. 

From Fig. 5.8, it can be seen that the smaller the slit width, the steeper 
the field decay. The near-fields are localized within a distance of y ^ d from 
the surface. These near-field distributions at the slit are similar to those for 
small apertures used as optical near-field probes in SNOM [16]. 
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0.5X 




(b) Experiment 



Fig. 5.7. (a) Calculated and (b) measured field distributions of the electric field 
component Ex in the x-direction on a slit of width d = A/4. Ex is normalized to the 
incident field Eix and the positions x and y are also normalized to the wavelength 



5.2.2 Optimum Slit Width 

Using the theoretical near-field distributions, energy changes of electrons 
passing close to the slit surface were estimated through computer simulation. 
The calculational model is the same as the previous one shown in Fig. 5.5. 
Electrons with velocity v\ move in the x-direction at distance y\ from the 
slit surface. All field components in the near- field, i.e., electric fields E^^ Ey 
and magnetic field were taken into account in the calculation. The total 
energy changes of the electrons were determined by integrating small energy 
changes with the Lorentz force over a small distance along the electron tra- 
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Fig. 5.8. Comparison between calculated and measured field intensities as a func- 
tion of the distance y away from the slit surface for different slit widths at x = 0. 
The field intensities Ex are normalized to Eix without the slit 




Fig. 5.9. Calculated maximum energy spread of electrons as a function of their 
transit angles 0 for different electron velocities (3 = Ui/c from 0.3 to 0.6 at y\ = 0.01 A. 
In the figure, the calculated energy spreads are normalized to the maximum value 
in the curve for /3 — 0.6 at ^ ~ 2.57t 



jectory [17]. The integral length of ten times the slit width was chosen to 
cover the near-field region of the slit completely. 

Maximum energy changes (spread) AW for electrons with different veloc- 
ities /? at = O.OIA were calculated as a function of slit width d. Results are 
shown in Fig. 5.9. The abscissa is the transit angle 9 = kod/f3 corresponding 
to d. The ordinate is AW normalized to the peak in the curve for j3 = 0.6. 
In Fig. 5.9, the AW curves have two peaks, at 6 ^ 0.57t and 2.57t. Note that 
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these optimum transit angles do not depend on (3. From the transition angles, 
it was found that optimum slit widths d exist for (3 and are given by 

d = , (5.10) 

where m is an integer. Comparing (5.10) with (5.9), it follows that the opti- 
mum width in the slit is narrower than that in the metal film gap by /3A/4. 
This could be due to the non-uniform distribution of the near-held over the 
slit. In Fig. 5.9, the second peaks at wider widths are always larger than the 
hrst peaks at narrower widths. When the slit width becomes wider, the held 
intensity on the slit decreases as shown in Fig. 5.8 and the effective interaction 
length increases. The results thus indicate that the increase in the interaction 
length is a dominant effect for AW, at least up to 9 = 2.57t. 



5.2.3 Phase Matching Condition 

In Sect. 5.1.1, we said that electrons only interact with an evanescent wave 
contained in an optical near-held on a microslit. From energy and momentum 
conservation, the interaction condition, i.e., the phase matching condition 
between electron and evanescent wave, was derived as shown in (5.6). Let us 
check the validity of this condition using a computer simulation. 

Figure 5.10 shows the calculated wave number {k) spectra of the near- 
helds Ex on the slits with d/X — 0.38, 0.5, an 0.62 at y{ — 0.01 A. The abscissa 
is the wave number kx normalized to fco and the corresponding electron veloc- 
ity calculated from (5.6). Amplitudes of wave components are normalized to 
the maximum value in the curve for d/\ — 0.38 at kx = 0. The wave compo- 
nents with kx > ko in Fig. 5.10 are evanescent waves that cannot propagate 
in free space because these waves have imaginary kz. 

In Fig. 5.10, the thick solid curve represents the calculated AW for the 
slit with d/X = 0.62. In the calculation, a CO 2 laser with power density 
10^ W/cm^ at A = 10.6 |im was taken as the incident wave. This power den- 
sity corresponds to a 10 kW output power focused onto an area of diameter 
100 pm. As can be seen from Fig. 5.10, AW is proportional to the ampli- 
tude of the evanescent wave. Calculation results indicate that electrons with 
initial velocity V[ interact only with the evanescent wave of phase velocity 
Vp{= uj/kx) = so that the interaction condition derived from conservation 
theory is valid. 

5.2.4 Interaction Space 

Figure 5.11 shows the calculated AW for different electron velocities as a 
function of the position y from the slit surface. Slits having the optimum slit 
widths oi d/X — 0.38, 0.5, 0.62 for (3 = 0.3, 0.4, 0.5, respectively, were used 
for calculation. Other parameters used are the same as before. 
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Electnon velocity, ^=v//c 



Fig. 5.10. Calculated wave number spectra for the near-field distributions of Ex on 
the slits with different widths d/X = 0.38, 0.5, and 0.62. In the figure, kx is the wave 
number in the a:-direction and j3 is the electron velocity calculated by substituting 
kx for kev in (5.5). The thick curve represents the maximum electron energy spread 
for different initial velocities from (3 = 0.25 to 0.8 which corresponds to an initial 
electron energy of 17-340 keV 




yf X 



Fig. 5.11. Calculated maximum electron energy spreads as a function oi y/\ for 
different electron velocities j3 = 0.5, 0.4, and 0.3 



In Fig. 5.11, when y increases from zero to 0.5A, AW falls off exponentially 
to near zero. Since AW is proportional to the field intensity at the slit, these 
curves represent effective field distributions over the slit. The shapes of the 
field distributions in Fig. 5.11 are quite different from those of the near-fields 
shown in Fig. 5.8. The field distributions are those of evanescent waves inter- 
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acting with the electrons. The field decay constant a of the evanescent wave 
in the ^-direction is given by 



a = k{) 




(5.11) 



This equation is easily derived from the relationship, = k^^ + ky^. Using 
(5.6) and (5.11), values of a for /3 = 0.3, 0.4, 0.5 are calculated to be 3.2fco, 
2.3fco, 1.7/co, respectively. The decay constants agree with those estimated 
from the curves shown in Fig. 5.11. 

The calculation results shown in Figs. 5.10 and 5.11 suggest that an elec- 
tron beam could be used to measure a A:-spectrum of an optical near-field 
appearing on a nano-object. This is a distinct feature of near-field microscopy 
with an electron beam as compared with conventional SNOMs. 

The interaction space of the slit can be defined diS = 1/ a because the 
field intensity of the evanescent wave falls off by e~^ . From (5.6) and (5.11), 
ye is expressed as 



_ A p 
“ 27T yi-/?2 



(5.12) 



As can be seen from this equation, the interaction space in the slit circuit 
is strongly limited, particularly for a low energy electron beam. For /3 = 0.5 
(~ 80keV) and A = 10.6 pm, ye is about 1pm. The electrons must therefore 
pass very close to the slit surface to obtain significant energy exchanges with 
the light. 

In practical measurements, the interaction space at the slit also depends 
on the field intensity of the evanescent wave. The energy changes of electrons 
are usually measured by an electron-energy analyzer with finite resolution. 
If the resolution is 1 eV, only electrons that gain or lose an energy of more 
than 1 eV can be detected. Therefore, the effective interaction space where 
detectable electrons are produced widens as the field intensity increases. For 
instance, the effective space in the curve for /? = 0.5 in Fig. 5.11 is about 
3.2 pm. 



5.3 Modulation with Infrared Evanescent Waves 

5.3.1 Experimental Setup 

The metal microslit circuit was demonstrated experimentally at the infrared 
wavelength of 10.6 pm [18]. The experimental setup is shown in Fig. 5.12. 
An electromechanical Q-switched (EMQ) CO 2 laser [19] oscillates in TEMqo 
mode at A = 10.6 pm. This laser provides regular pulses with an output peak 
power of 11 kW and a width of 140 ns at a high repetition rate of 1 kpps. The 
laser beam focused through a ZnSe lens has a diameter of about 200 pm at 
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Fig. 5.12. Experimental setup for measuring the electron-light interaction in the 
infrared wave region 



the slit surface. The metal slit consists of two polished copper blocks, and 
the width is 8.4 |im. The electron beam has a diameter of 240 |im. The slot 
aperture placed in front of the slit confines the beam area on the slit to 10 fim 
in height and 100 |Lim in width. The initial energy W[ of the electron beam was 
adjusted between 40 and 90keV. The electron energy was measured using a 
retarding field analyzer [20]. 

The pulsed laser output modulates the energy of the electron beam so that 
the number of electrons through the analyzer changes during the pulse. The 
electrons passed through the energy analyzer were detected by a secondary 
electron multiplier (collector) connected to a gated counter which is triggered 
by the laser pulse. 

5.3.2 Electron Energy Spectrum 

Figures 5.13a and b show typical temporal changes in the laser pulse and the 
corresponding response from the gated counter for an electron beam with 
W[ = 80keV. The collector current was i = InA and the peak power of the 
laser was P[ = 9kW. The output response was measured by the counter in 
box-car averager mode with a time resolution of 20 ns and integration time 
2 s (accumulation number = 2 000). In this experiment, Vf was set to —5.5 V 
so that the counts in Fig. 5.13 represent the number of electrons gaining more 
energy than \qVi\ from the laser. 

As seen from Fig. 5.13, the shape of the output response differs consider- 
ably from that of the laser pulse. The response time of 380 ns in Fig. 5.13b 
is much longer than the laser pulse width of 140 ns. This output response 
results from the fact that the number of electrons interacting with the laser 
light is proportional to the laser field, not to the power. Note that a high 
signal-to-noise ratio (S/N) of more than 1 000 is achieved for the measure- 
ments. 
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Fig. 5.13. Measured temporal changes of (a) the input laser pulse and (b) the 
output response from the gated counter for electrons with Wi = SOkeV at A = 
10.6 pm 



Figure 5.14a shows the measured energy spectra of electrons with (A) and 
without {B) laser illumination, while Fig. 5.14b shows the difference between 
the two spectra A — B.ln the experiment, W[ = 80 keV and P[ = 10 kW were 
used. The electron current was held below 1 pA to maintain a linear response 
from the electron multiplier even for a larger current input at Ff > 0. To 
increase S/N in the measurement, a longer gate width of 1.5 ps was chosen 
and an integration time of 10 s was set in the counter. 

From the measured spectrum B in Fig. 5.14a, it can be seen that the en- 
ergy analyzer has resolution better than 0.8 eV for an 80keV electron beam. 
The output count decreases gradually as Vf increases from +1 eV, due to the 
dispersion of the energy analyzer [20]. When the laser beam irradiates the 
electrons, the spectrum B changes to the spectrum A with a wider energy 
spread. Spectrum A still contains a number of electrons that have not inter- 
acted with the light. In order to remove these electrons and the dispersion 
effect from the measured spectrum A, the output counts in B were subtracted 
from those in A. Figure 5.14b thus indicates the energy spectrum only for 
the electrons that interacted with the light. From Fig. 5.14b, it can be seen 
that the 10 kW laser beam can give an energy spread of more than ±10 eV 
to the electrons. 

Since the energy analyzer passes all the higher energy electrons, it is 
expected that for large Vf the counts of electrons with laser illumination 
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Fig. 5.14. Measured energy spectra for electrons with and without laser illumi- 
nation (a) and the difference between the two spectra (b) for a laser beam with 
peak power 10 kW. The zero filter bias corresponds to the initial electron energy of 
SOkeV 



should be the same as those for electrons without laser illumination. However, 
as shown in Fig. 5.14b, the counts with laser illumination are slightly smaller 
than those without laser illumination even at Vf > +10 V. This is due to 
deflection of the electron beam by the laser illumination. Consequently, a 
part of the electron beam has been clipped by the aperture before reaching 
the collector. 

In Figs. 5.14a and b, about 42 000 electrons passed the slit, and about 
9 000 electrons among them interacted with the laser beam. Since the height 
of the electron beam at the slit is 10 |Lim, this ratio of signal electrons to total 
electrons implies that the interaction space of the slit is about 2 |im, which 
agrees with the theoretical prediction described in Sect. 5.2.4. 

5.3.3 Laser Field Dependence 

In order to conhrm that the electron energy spread is caused by the laser held 
on the slit, we measured the dependence of electron energy spreads on the 
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Fig. 5.15. Measured electron energy spreads versus polarization angle 0 of the 
incident laser beam. The solid curve is the theoretically predicted change with 0. 
In the inset^ E[ is the electric vector of the laser beam 




Laser peak power (kW> 

Fig. 5.16. Measured energy spreads of the electrons with an energy of SOkeV as 
a function of the peak power of the incident laser beam. The solid curve is the 
theoretical fit to the measured data 



polarization of the laser beam and the power. The results for the polarization 
angles 6 are shown in Fig. 5.15. The ordinate is the electron energy spread AW 
normalized to the maximum value at 0 = 0. The experimental parameters 
used are the same as the previous ones. The solid curve is the theoretical 
change in AW calculated assuming the relation AW ex E[ cos 0, where E[ is 
the field intensity of the incident laser beam. Measurements agree well with 
theory. 

Figure 5.16 shows the measured dependence of AW on the laser power. 
The solid curve is the theoretical fit to the measurements calculated by as- 
suming that ALE is proportional to the field intensity of the incident wave, 
i.e., the square root of the laser power. 
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In Fig. 5.16, the measured AW is 13 eV at Pi = 10 kW, to be compared 
with the theoretical value of 22 eV estimated through the computer simu- 
lation. The reduction in the electron energy spread arises from differences 
between the actual slit and the ideal one used in the theoretical analysis. 
Since the actual slit consisted of two thick copper blocks with finite conduc- 
tance, the amplitude of the evanescent wave may be small compared to the 
theoretical one. The results shown in Figs. 5.15 and 5.16 show that electrons 
are accelerated or decelerated by the infrared laser field at the microslit. 

5.3.4 Electron Velocity Dependence 

Figure 5.17 shows the measured and calculated energy spreads AW of the 
electrons as a function of the initial electron energy W\. The initial electron 
velocities j3 corresponding to W\ are also indicated in the figure. The the- 
oretical energy changes were calculated for a slit width of 7.2 |im to obtain 
the best fit to the measurements. In the figure, AW was normalized to the 
maximum values of 15 eV for the measurement and 34 eV for the theory, with 
W{ = 90 keV. Due to the limitations of the electron gun, the measurements 
were made with W{ above 40 keV. 

As described in Sect. 5.2.3, the measured variation of the electron energy 
spread shown in Fig. 5.17 represents the wave number spectrum of the optical 
near-held at the slit. The measured wave number spectrum implies that the 
effective width of the slit is not 8.4 pm as measured by the optical microscope, 
but 7.2 pm for the laser wave at A = 10.6 pm. 

As can be seen from Fig. 5.17, electrons with lower energy, less than 
50keV, were modulated by the laser beam in the microslit. The electron 




Initial electron energy (keV) 

Fig. 5.17. Measured electron energy spreads as a function of the initial electron 
energy. In the hgure, the electron velocities (3 corresponding to the electron energy 
are also indicated. The solid curve is the theoretically predicted electron energy 
spread for the slit with width 7.2 pm 
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energy spread is about 4eV at W[ — 40keV which is still enough to be able 
to measure the interaction. These results indicate that nonrelativistic elec- 
trons can be modulated using the metal microslit at optical frequencies, as 
predicted by theory. 

5.4 Interaction with Visible Light 

Figure 5.18 shows a conceptual drawing of the experimental system for the 
interaction of nonrelativistic electrons with light at shorter wavelengths than 
about 1.2 pm. The metal microslit is fabricated at the end of an optical fiber so 
that the laser beam is guided to the slit without any precise adjustment. Since 
light with shorter wavelengths has a photon energy higher than 1 eV, and the 
electron energy analyzer can resolve the energy changes of the electrons due 
to the photons, a quantum step in the electron energy is detectable. 

On the basis of the theoretical and experimental results, the number of 
signal electrons in the interaction was estimated in both classical and quan- 
tum treatments, assuming a laser with A = 780 nm and output power of 
30 mW, and an electron beam with velocity P = 0.5 and density ImA/cm^. 
It was also assumed that the diameter of the laser beam in the fiber was 6 pm. 

According to the classical theory described in Sect. 5.2, the energy change 
of the electrons is estimated to be less than 0.1 eV. This energy change is 
too small to allow observation of the interaction in practical experiments. In 
order to give more than 1 eV to the electrons, a higher output power than 1 W 
is required for the laser. However, this laser power is available in experiments. 

In contrast to the conclusion from classical considerations, a different 
result is derived from a quantum mechanical treatment. Assuming the tran- 
sition rates in the metal film gap and the microslit to be the same, the tran- 
sition rate of about 1 x 10~^/s, which is the rate for an electron to absorb a 
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Fig. 5.18. Conceptual drawing of the experimental setup for measuring the 
electron-light interaction in the near-infrared and visible light regions 
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Fig. 5.19. Metal microslit fabricated at the center of the end of a single-mode 
optical fiber, (a) Ridge structure with a taper angle of 80°, (b) metal slit with a 
width of 270 nm fabricated on the ridge structure 



laser photon, was obtained using (5.7). From (5.10) and (5.12), the optimum 
slit width and the interaction space are 490 and 72 nm, respectively. Then 
2.7 X 10^ electrons per second can pass through the space of 72 nm x 6|Lim 
at the slit. The number of signal electrons is thus more than 20000 particles 
per second. This number is enough to observe the electron-light interaction 
experimentally. Since the interaction behavior between electron and light in 
a transition region from the classical to quantum regimes is still vague, the 
final conclusion for the above discussions should be drawn from experimental 
evidence. 

In the experiment at visible light wavelengths, a key device is the metal 
microslit with submicron width. The microslit is fabricated at the top of the 
ridge, as shown in Fig. 5.18. This ridge structure is needed to avoid a collision 
of the electrons with the slit due to the image force acting on the electrons in 
the proximity of the slit surface. Such a microslit can be fabricated using a 
chemical etching technique in a similar way to conventional near-field probes 
with nanometric apertures [21,22]. 

Figure 5.19 shows a prototype of the microslit fabricated at the center of 
the core of an optical fiber with a 125 fim diameter. Figure 5.19a shows the 
ridge structure and Fig. 5.19b the same ridge with aluminum coating and a 
270 nm gap at the top. The ridge structure has a length of 5.5 |im for a core 
diameter of 8.8 jam, a height of 7.2 |im, and a taper angle of 80°. Both the 
flatness along the length of the ridge and the radius of curvature at the top 
of the ridge are less than 30 nm. Another probe with a 180 nm wide slit was 
also fabricated using gold instead of aluminum as a coating metal. 
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5.5 Conclusion 

Energy modulation of an electron beam with light waves was analyzed theo- 
retically using classical and quantum treatments. Three gap circuits, a metal 
film gap, a dielectric film, and a metal microslit were discussed as possible 
interaction circuits. Theory showed that the transition rates of electrons in 
the gap circuits were high enough to demonstrate the electron-light inter- 
action experimentally. To achieve the interaction with an optical near-field, 
the metal microslit structure was adopted as the interaction circuit. Classical 
theory revealed that electrons passing through the slit can interact with an 
evanescent wave contained in the near-field induced on the microslit by light 
illumination. The relationship between the slit width, electron velocity, and 
light wavelength for the strongest interaction was determined through com- 
puter simulation. The microslit interaction circuit was successfully demon- 
strated experimentally at the infrared wavelength of 10.6 pm. Experiments 
showed that the energy of electrons can be changed by more than 10 eV with 
10 kW laser illumination for a nonrelativistic electron beam with low energy, 
less than 50keV. These results provide experimental backing for theoretical 
predictions. From the theoretical and experimental results, it was concluded 
that an optical near-field on a small object can be detected using an electron 
beam, so that a new type of near-field microscopy with free electrons may be 
feasible. 
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6.1 What is the Issue? 

Tunneling is undoubtedly the most important quantum phenomena with a 
wide range of applications. Like other physical phenomena, it should be pos- 
sible to understand tunneling phenomena from a spatial and also from a 
temporal point of view (see Fig. 6.1). Of interest from a spatial point of view 
is the question as to how many particles initially on one side of the barrier 
will finally be found on the other side of the barrier at sufficiently later times. 
We are thus led to the concept of tunneling probability. Of interest from a 
temporal point of view would be the tunneling time, that is, the time a par- 
ticle ‘spends’ in the barrier region. From a spacetime point of view, it seems 
natural to characterize tunneling phenomena according to the two quantities 
tunneling probability and tunneling time. 

The lifetime is well known as a characteristic time scale associated with 
tunneling. The tunneling time must not be confused with the lifetime (see 
Fig. 6.2). When a system tunnels out of a metastable state, for example, the 
lifetime is intuitively the measure of the time during which the particle is 
trapped in the potential well (classically allowed region), while the tunneling 
time is thought to be the time that the particle is under the potential barrier 
(classically forbidden region). 

In applying tunneling phenomena, only the tunneling probability has been 
utilized. It would be safe to say that the tunneling time plays no role today in 
tools and devices based on tunneling. This raises the expectation that some- 
thing new could be built or the performance of existing tunneling devices 
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Fig. 6.1. Space-time picture of tunneling 
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Fig. 6.2. Tunneling out of a metastable state: lifetime and tunneling time 



could be improved if the tunneling time were utilized to control tunneling 
phenomena. This indicates the technological importance of tunneling time. 
Another point is that the tunneling spectroscopy as it exists today extracts 
information from the relevant system by measuring the tunneling probability 
(or tunneling current). If one could find a way to probe the system via a 
properly designed tunneling time measurement, a new area of research that 
deserved the name Tunneling-time spectroscopy’ would appear and the tra- 
ditional one would then be called Tunneling-probability spectroscopy’. This 
is where we find the importance of the study of tunneling time from an aca- 
demic point of view. Scientists have been utilizing only one aspect of tunneling 
and that the utilization of another aspect, the time of tunneling, seems to 
be a fascinating direction for tunneling-related science and technology in the 
twenty-first century. 

The barrier to overcome is high, however. There is no agreement as to 
what the tunneling time is and even as to the meaning of the concept of 
tunneling time. Regarding tunneling probability, its concept is very clear as 
described above, though the actual calculation of the probability is not always 
easy. On the other hand, the concept of tunneling time is not clear despite 
extended debate over the past 70 years. Therefore, the first thing to be done is 
to put an end to the long dispute and establish the concept of tunneling time. 
The purpose of this chapter is to present the problem of tunneling time for 
those readers who routinely deal with tunneling phenomena in their research 
but only in terms of tunneling probability. It is hoped that this chapter will 
inspire such readers to begin thinking about tunneling time. 

In Sect. 6.2, we survey a variety of ideas that have been proposed to de- 
fine and/or measure the tunneling time (see, for example, [1-8] for more 
detailed reviews of the problem). We will see that different methods of defin- 
ing/measuring the time lead to different answers. All these methods are based 
on physically plausible ideas and often give reasonable answers when applied 
to non-tunneling regimes such as free propagation or transmission over a bar- 
rier or well. This implies that such ‘physically plausible ideas’ were reached 
by implicitly imaging the dynamics in non-tunneling regimes where classical 
mechanical intuition based on a particle picture works. However, tunneling 
is a purely quantum mechanical phenomena where the validity of such an 
intuition is no longer warranted. We must therefore treat the problem in a 
way securely rooted in the foundations of quantum mechanics. 
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A pessimist, being tired of repeated arguments with no convergence in 
the results, would argue that it does not make sense to talk about tunneling 
time. It should make sense, however, simply because tunneling is a physical 
phenomena that occurs, like other physical phenomena, in space-time. In 
my opinion, an approach that is free from unwarranted ‘plausible arguments’ 
but is instead securely rooted in the foundations of quantum mechanics would 
change this controversial situation. The author reported an effort in this di- 
rection in [9]. Section 6.3 gives an overview of this approach, which is intended 
for those who were once interested in the tunneling time problem but have 
stopped thinking about it, under the negative impression that the problem 
is not meaningful. 

The question asked in Sect. 6.3 is more fundamental than the questions 
asked in most other approaches. That is, the definability of the probability 
distribution of tunneling time^ is questioned and investigated, while many 
other arguments implicitly assume the existence of the distribution and try 
to construct it in a plausible way. Such an assumption does not seem to be 
founded. What must be clarified first is the definability of the distribution. 
Until recently, however, it was not clear how one could make such a clari- 
fication in quantum mechanics. It will be argued here that the consistency 
criterion provided by the consistent (or decoherent) histories approach to 
quantum mechanics is best suited to the purpose. 

The consistent histories approach (CHA) is a general theory that deals 
with the question as to whether or not quantum probabilities are definable 
for a given set of alternatives.^ CHA was developed by Griffiths [10], Omnes 
[11,12], Gell-Mann and Hartle [13,14], and also by the present author with 
Takagi [15,16]. 

As we shall see, a negative result follows from the consistency criterion, 
i.e., the probability distribution is not definable. The implicit assumption 
behind many arguments is therefore incorrect. We should be able to convince 
ourselves of this negative result by recalling the situation we face in the two- 
slit experiment. Since the alternative that the particle goes through the upper 
slit and the other alternative that it goes through the lower slit interfere,^ 

^ Quantum mechanically, it is unlikely that the value of tunneling time is a constant 
that is uniquely determined from the initial condition of the particle and the 
barrier shape. Rather, a distribution of the values of tunneling time seems to be 
more natural to consider. 

^ Those who have studied CHA will notice that the equations that appear in 
Sect. 6.3 do not look like CHA equations. This is because, unlike most CHA 
studies where the relevant alternatives are defined by the products of projection 
operators at different times, Feynman paths are used to define alternatives in the 
present study. 

^ The two alternatives, which are classically mutually exclusive, are said to be 
interfering alternatives. 
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probabilities cannot be defined for the two alternatives in a consistent way.^ 
Likewise, the alternative that the particle spends a certain amount of time 
in the barrier region and another alternative that it spends another amount 
of time in the region interfere,^ making it impossible to define probabilities 
for the alternatives in a consistent way. 

The role of CHA in the above story is to provide a rigorous definition for 
the measure of the interference (decoherence functional). What the author 
has found using the measure is that interference never vanishes in tunneling 
regimes. This negative conclusion immediately raises the following question: 
what can then be said about tunneling time? This is the second question 
treated in Sect. 6.3. The answer given there is that the range of values of 
tunneling time is ‘speakable’. This is also understandable in analogy with the 
‘one-slit experiment’. Readers not interested in technical details can skip all 
equations in Sect. 6.3 but are instead encouraged to understand the results 
with the help of the two analogies employed there (see Fig. 6.7). The author 
does not claim to have solved the long-standing problem. Rather, the problem 
is simply redefined and a new starting point is set for constructive discussion. 
What physical results are obtained from this new starting point is mostly left 
to future studies. It is worth repeating, however, that the most important 
thing for the tunneling time problem is to have a sensible definition of the 
problem, and such a definition has, I believe, been presented. The long history 
of controversy does not mean that the tunneling time problem is meaningless. 

Throughout the text, our analysis is restricted to the case of a rectangular 
potential barrier in one-dimensional space. Since our interest is a conceptual 
issue, this would be a reasonable choice. As stressed by Stpvneng [17], of 
course, a fundamental study is not meaningful if it cannot be related to ex- 
periment, and ultimately to applications. Sadly, it is too early to speak of 
applications of tunneling time. For the moment we can only focus on con- 
ceptual issues. Once a foundation is established, applications should become 
possible. For most readers, this will be frustrating, but we should recall that 
the Esaki diode was invented in 1957 and STM in 1982, while tunneling 
phenomena were already known in the 1920s. 

6.2 Well-known Approaches 

6.2.1 Notation 

There are many theoretical ideas for defining tunneling time. Only a few of 
these will be described here. Since the derivation of all the details is not 
the purpose here, intermediate equations will be omitted as far as possible. 

^ In such a way that the superposition principle for amplitudes (wave nature) and 
the sum rule for probabilities (particle nature) do not conflict. 

^ This interference may be said to be a temporal interference, while the interference 
in the two-slit experiment is a spatial interference. 
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To avoid complication, descriptions will be made for a square barrier in one 
dimension, although some of the ideas have been analyzed in more general 
circumstances in the literature. The following descriptions are not necessarily 
faithful to the original notations. 

The square barrier of height Vq is assumed to cover the region a < x < b. 
The Hamiltonian of the particle is simply 



H = - 



2m dx‘^ 



+ ho 0ab{x) , 



( 6 . 1 ) 



where m is the mass of the particle and 



Oab{x) S I ^ 



for a < X < 5 , 
otherwise . 



(6.2) 



Some approaches use wave packets and others use plane waves. In either case, 
what is basic to the analysis is the normalized eigenfunction Uk{x) which takes 
the following form outside the barrier: 



Uk{x) 



for rc < a , 

^ v27T 

for a: >6, 

I 



(6.3) 



where k > 0^ d = b—a^ and the complex coefficients R and T are the reflection 
and transmission amplitudes, respectively, satisfying |i?p + |Tp = 1. Only T 
is important in our analysis. It is given by 



T = e 



— ikd 





\Z^2 - ky^ 




-1 

(6.4) 



with fcvb = \/2mVo/h. The decomposition of T into its modulus and phase 
is often useful: 



T= |T|e^^ = |T|e^^e~^^^. (6.5) 

For a particle described by a wave packet ?F(x, t) (initially localized in x < a 
and moving towards the barrier) , the tunneling (or transmission) probability 
Pi is given by 

Pt= lim [ dx\<P{x,t)f= [dk\ip{k)\^\T{k)f (6.6) 

Jb J 

^ \T{ko)f as \^{k)f ^ 6{k - ko) , (6.7) 

where (p{k) is the fc-space wave function, which is assumed to vanish for < 0 
(to make a right-moving wave packet) and also for k > ^/2m^/h (otherwise 
the particle can go over the barrier).^ Equation (6.7) gives the tunneling prob- 
ability for the stationary case with wave number ko. The reflection probability 
Pj. can be calculated from Pj. = 1 — P^. 

® Note that the tunneling probability is originally defined as the probability of 
finding a particle somewhere on the right side of the barrier at sufficiently later 
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Fig. 6.3. Tunneling time from the motion of wave packets 



6.2.2 Following the Motion of a Wave Packet 

This is probably the simplest idea that comes to mind from a particle point of 
view. As an incident wave packet, we consider a wave packet that is narrow in 
fc-space and thus wide in real space. This corresponds to a particle of nearly 
monochromatic kinetic energy. Explicitly, 

/ dk 

( 6 . 8 ) 

where ^ is such that \ 2 p{k — ko)\ is sharply peaked around k = ko > 0 
and satisfies Jdfc|'^P = 1. At the initial time t = 0, the incident packet is 
supposed to be localized well to the left of the barrier and is moving towards 
the barrier with its peak velocity vq = hko/m. At sufficiently later times, the 
wave function is made up of two spatially non-overlapping parts, the reflected 
and the transmitted packets. The reflected packet has a clear peak in x ^ a 
that leaves the barrier with the velocity — ug, while the transmitted packet 
has a clear peak in x 6 leaving the barrier with the velocity Vq. 

Now, let us focus on the peak position x{t) of the transmitted packet and 
compare it with the peak position Xf(t) of the free wave packet that would 
have evolved under the same initial condition without the barrier. Arguing 
that the spatial shift of the peak, Ax(t) = x(t) — Xf (t), is due to the extra time 
spent by the tunneling particle inside the barrier, we are led to the following 
expression for the tunneling time, which is called the phase time^ 

'Tphase = 'Tfree + ^x{t)/vQ , (6.9) 

where rnee = djvo is the time a classical free particle with velocity vq spends 
in the spatial region a < x <b. Now, the transmitted packet is expressed as 

J y27r 

times. In text books, the tunneling probability for the stationary case is usually 
defined as the ratio of transmitted flux to incident flux (with a velocity ratio factor 
when the barrier is not symmetrical). This procedure is, however, outside the 
standard prescription for calculating probabilities using amplitude and operators. 
The procedure presented here is completely within the standard prescription. 
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Because of the peaked nature of — ko)^ we can evaluate (6.10) using the 
stationary phase approximation. We then find that the peak position of the 
transmitted packet moves as x{t) = vot — 9\k)\k=ko^ where 9 is the phase 
defined by (6.5) and the prime denotes d/dk^ while the peak of the free wave 
packet moves as Xf(t) = vot. Using these results in (6.9), we have 



"^phase ^ 



d(j) 



E=Eo 



( 6 . 11 ) 



where E = f?k ^ = k?k\/2m, and (j) is defined by (6.5). Essentially 
the same result can be obtained by following, instead of the peak of the 
transmitted packet, a suitably defined ‘center of gravity’. 

The phase time is undoubtedly useful in a qualitative understanding of 
the motion of the transmitted packet. However, we cannot consider it to be 
the time taken by a packet peak to traverse the barrier because the peak of 
the transmitted packet is not clear near the barrier. In the above derivation, 
the clear peak motion in the region x ^ h was linearly extrapolated back to 
X — b. Near the barrier, the extrapolated peak motion, which is the basis of 
(6.11), cannot be found in the real dynamics of the wave function. The phase 
time is useful but is not the answer to the tunneling time problem. 



6.2.3 Velocity Under the Barrier 

The phase time method is based on the evaluation of (6.10) by the stationary 
phase approximation. What about more sophisticated evaluations to access 
the motion inside the barrier? Using a contour integral representation of 
the time dependent wave function, Stevens [18] analyzed how a pulse wave 
that initially has a sharp wave front (product of a plane wave and a step 
function) propagates under the step potential region of height Vq • Evaluating 
the integral by the steepest descent method, he concluded that the wave 
front propagates at a velocity ve = y^2(Vo — E)jm7 This conclusion had an 
impact because a clear velocity was assigned to the motion of the particle in 
the classically forbidden region. However, Teranishi et al. [19] reinvestigated 
the problem with more care in the analytical steps and found that Stevens’ 
conclusion is incorrect. For the case of a square barrier, this negative result 
was also confirmed numerically by Jauho and Jonson [20] and analytically by 
Ranfagni et al. [21]. Moretti [22] arrived at a conclusion similar to Stevens 
for a square barrier, but the result attributing a velocity to the particle under 
the barrier was questioned by Muga et al. [23]. 

Although Stevens’s original result is not correct, it motivated others to 
develop an exact analytical expression^ for the time-dependent wave func- 
tion. The present author considers this exact solution approach to be very 

^ This implies that the tunneling time for an opaque rectangular barrier is d/vE^ 
which is identical to the Biittiker-Landauer time to be discussed in Sect. 6.2.5. 

® In terms of the critical points (such as poles) of the transmission and reflection 
amplitudes. 
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important because it allows a non-perturbative analysis of tunneling in an 
analytical manner.^ For this approach, see [23-26]. Of course, if one aims sim- 
ply to visualize the time evolution of the wave function, one can numerically 
solve the Schrodinger equation from the beginning. There are many papers 
on this issue including the very early one by Goldberg et al. [27]. 



6.2.4 The Larmor Clock 



It may be said that the phase time is based on an equation of type x vtm 
real space, while the Larmor times discussed here are based on an analogous 
equation of type Lp utm spin space with ip being a spin rotation angle and 
u an angular velocity. There is, however, an important difference between the 
two. That is, the phase time method is based on x = vt outside the barrier, 
while the Larmor clock method is based on ip = uot inside the barrier. In 
the phase time method, x and v (to be precise. Ax and Xg), and therefore 
the resultant t, were quantities associated with the motion of wave packets 
outside the barrier. If we attempt to use an equation of type x = vt inside 
the barrier, which is preferable since we are interested in the time inside the 
barrier, we immediately face the difficulty that v is not meaningful (or it is 
imaginary) in the barrier. This kind of difficulty is absent in the Larmor clock 
method. 

We consider the tunneling of a particle with spin 1/2. The motion of 
the particle is supposed to be one-dimensional along the y-axis. A small and 
homogeneous magnetic field Bqz is applied only in the barrier region {z is 
the unit vector in the ^-direction) . The Hamiltonian now has an additional 
term 



huji, 

2 



(^zOabiy) 



on the right-hand side of (6.1) (replac xhy y). Here is a Pauli spin matrix. 
The initial spin is supposed to be polarized in the x-direction, and the initial 




Fig. 6.4. The Larmor clock method 



® Of course, the critical points must be obtained numerically. 
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orbital wave function is assumed to be monochromatic with wave number fco 
(energy Eq). 

As the particle enters the barrier, the spin starts to rotate due to the 
magnetic field. If the particle is at rest in the magnetic field, the spin rotates 
in the x^-plane through an angle cjl = per emit time, where g 

is the gyromagnetic ratio, and jib the Bohr magneton. This is a Larmor 
precession. Hence the rotation angle divided by ujb gives the time the particle 
is located in the magnetic field. In the tunneling case, it thus also seems 
natural to calculate the tunneling time by dividing the spin rotation angle 
accumulated in the barrier by (Jl- However, as pointed out by Buttiker [28], 
the motion of the spin is no longer a Larmor precession in the xy-plane in 
the tunneling case. Indeed, the spin rotates not only in the x^-plane but also 
in the xz-plane. Only the rotation in the x^z-plane was considered in earlier 
papers [29,30]. Therefore, we obtain two different tunneling times depending 
on which rotation is analyzed with an equation of type (f = out. 

That two rotations arise can be understood as follows. The initial spin 
may be written as the superposition of a spin-up and a spin-down state (‘up’ 
and ‘down’ are with respect to the 2 :-direction) . Due to Zeeman splitting, the 
spin- up component sees a modified barrier of height Vq — Uujb/2^ while the 
spin-down component sees a modified barrier of height Vq + There- 

fore, the spin-up component is preferably transmitted since the barrier is 
lower, and as a result the spin of the transmitted particle acquires a non- 
zero ^ component. Of course, the spin also acquires a y component due to 
the conventional Larmor precession. The expectation values of Sz and Sy for 
transmitted particles can be written in terms of the transmission amplitudes 
for the square barriers of height Vq i To the lowest order in they 

are found to be 





h 


^ dVo 


md 


(6.12) 


i^y) 


h 


a 

1 

11 


hko 


(6.13) 



where n = y^Vo ~ expressions in (6.12) and (6.13) are for 

opaque barriers. We thus end up with the two tunneling times Tz and Ty. 
One can of course take an average of the two tunneling times to define a 

single quantity such as tt = -h , which is Buttiker ’s traversal time 

(the subscript T in tt stands for ‘traversal’). However, there are many ways 
to define an average and, at least within the Larmor clock approach, a special 
status cannot be given to tt- The Larmor clock method does not give a unique 
answer to the tunneling time problem. 

In a narrower sense, Ty is called the Larmor time, but in a broader sense 
both Ty and are often called the Larmor times. Leavens call Ty the ‘spin- 
precession traversal time of Rybachenko’ and the ‘spin-rotation traversal 
time of Buttiker ’. 
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Fig. 6.5. Oscillating barrier of Biittiker and Landauer 



Leavens [31] generalized Biittiker ’s analysis of the Larmor clock to the 
case of arbitrary barriers by introducing an average height V of the potential, 
finding that the resultant Larmor times are also given by (6.12) and (6.13) 
with Vo replaced by V. Leavens [32] also made a generalization to the case 
of a uniform transverse magnetic field confined to an arbitrary part of an 
arbitrary barrier. Important findings are: 

• each of the resultant four Larmor times (r^ and Ty for reflection and 

transmission) is spatially additive, i.e., each satisfies (6.30) for 0 < < 

X2 < xs < dj 

• the reflection and transmission times are consistent with the local dwell 
time. For example, for reflection and for transmission satisfy (6.29). 

A series of precise measurements of Larmor times has been reported by a 
Japanese group. The results are, as they must be, in excellent agreement 
with theoretical predictions. See [33] and references therein. 



6.2.5 Time-Modulated Barriers 

The previous two approaches attempt to measure the tunneling time directly. 
The peak shift or spin rotation angles are measured and translated into time 
information through equations of the form x = vt or ^ = uot. In contrast, the 
approach to be discussed here attempts to measure the tunneling time indi- 
rectly. The tunneling particle is perturbed with a known time scale and the 
tunneling time is then estimated by finding such a time scale of perturbation 
that leads to apparent changes in the behavior of the tunneling particle. 

Biittiker and Landauer [34,35] studied an oscillatory barrier, where a 
small perturbation Vi cos(a;t) is added to the original static barrier. The 
Hamiltonian is therefore given by (6.1) with Vo replaced by Vq + Vi cos{ujt). 
Since the potential is time dependent, the energy of the transmitted parti- 
cle can differ from that of the incident particle, which is assumed to be Eq 
(monochromatic). The particle can emit or absorb energy quanta nhuj (n be- 
ing a positive integer) through the interaction with the barrier. Treating V\ 
as a perturbatio n and assuming further that huj Eo, huj Vq — Eq, and 
d^y2m{Vo — Eo)/h ^ 1, they analyzed the case of emission or absorption of 
a single quantum huo and calculated the tunneling probabilities and P_ 
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that the particle is transmitted at the sideband energies Eq^Tiuo and Eo — hcu, 
respectively. They showed for an opague barrier, with T — 

2 






exp(±27TTBL/^) - 1 



P. 



(6.14) 



7"BL 



md 

hK 



A 

Ve' 



ve = \/2(Vb - , 



(6.15) 



where P is the tunneling probability at the energy E'o, T is the modulation 
period (not the transmission amplitude), and ve is the absolute value of the 
imaginary velocity ^J2{Eo — Vo)/m under the barrier (here the subscript E 
on Ve stands for ‘Euclidean’). 

Let us consider P± as functions of the modulation period T and observe 
how the functions behave as T is varied. Equation (6.14) shows that P±/P 
is negligible for 2etblIT <C 1, while it is not for 2'ktbeIT > 1- Because of 
this crossover behavior at T ^ 27rrBL, we identify tbl as the tunneling time, 
apart from the numerical factor 2tt. For opaque barriers, we can write 



^BL 



,ain|T| 



(6.16) 



Normally, tbl is called the Biittiker-Landauer time, although Hauge and 
Stpvneng [2] called tt = ^ the Biittiker-Landauer time. For opaque 

barriers, tt ~ tbl since = tbl ^ 

As just mentioned, tbl coincides with obtained in the Larmor clock 
method. Under the name of ‘bounce time’, tbl also appears in the instanton 
technique, an imaginary-time path integral technique for calculating tunnel- 
ing probabilities. Does this popularity of tbl imply that it is the tunneling 
time? 

Before giving tbl a special status, it must be checked whether a different 
modulation of the barrier also gives tbl as a time scale characterizing some 
crossover behavior. To this end, Takagi [36] studied another time-modulated 
barrier model, where not only the height but also the width of the barrier 
changes in time (see Fig. 6.6). In particular, the tunneling through a squeezing 
potential was studied. The height V (t) and width d{t) change as 



d{t) = ii + t/T)d. 



(6.17) 




V{t) = VoHl+tlTj^ 






4 ► 

d{t)=do{l + tlT) 



Fig. 6.6. Linearly squeezing barrier of Takagi 
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A remarkable feature of this model is that the exact expression for an evolving 
wave packet can be obtained in terms of the scattering amplitudes. Since the 
potential changes with time scale T, the energy of the transmitted particle can 
differ from the incident energy. This time, however, the energy does not shift 
by an amount n27rh/T because the barrier modulation is not periodic. Thus, 
instead of P±/P^ we need a new measure that quantifies how the tunneling is 
affected by the barrier modulation. Paying attention to the barrier effect on 
the transmitted wave function in k space, rather than real space, Takagi used 
a = Apout/Apin — 1 as the measure, where Apin is the momentum uncertainty 
of the incident packet and 4ipout is that of the transmitted packet. Analyzing 
a as a function of T, he showed that for opaque barriers a is negligible for 
TTakagi/r < 1, while it is not for TTakagi/^" > 1, where 



_ 2m d‘^(f) 
rxakagi = X ^ 



k=ko 



Ti^ko 



(6.18) 



If we follow the spirit of Biittiker and Landauer, rxakagi should now be iden- 
tified as the tunneling time. 

Depending on how we modulate the barrier, we now get different an- 
swers. The difference between tbl and rxakagi is not a minor difference, 
either qualitatively or quantitatively. Like the Larmor clock approach, the 
modulated barrier approach also fails to give a unique answer to the tunnel- 
ing time problem. 

Note that the status oH/tbl as the crossover frequency in the oscillating 
barrier model has been questioned by some authors [38-40]. They numerically 
calculated the transmission probability as the function of uj but failed to find 
a crossover behavior at uo ~ ^Itbl^ 

Although Tbl is widely used as the tunneling time in the literature, its 
status must be considered as not yet established. 



6.2.6 Dwell Time and its Decomposition 

Applying the concept of collision time [41] in the scattering theory to the 
case of the one-dimensional stationary scattering problem, Biittiker [28] in- 
troduced the dwell time ra as the ratio of the number of particles within the 
barrier to the incident flux: 

Td{k) = ^ j dx \uk{x)f , (6.19) 

where Uk{x) is the eigenfunction with wave number k. As clearly stated by 
Landauer and Martin [1], the meaning of ra is simply ‘the time that the 
incident flux has to be turned on to provide the accumulated particle storage 



This is also supported by the study due to Truscott [37], who analyzed another 
time-modulated barrier but did not find tbl as the characteristic time in the 
sense of Biittiker and Landauer. 
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in the barrier’. In the literature (including [28]), is often called the ‘average’ 
dwell time. However, since it is not the case that is given as an average 
with respect to a probability distribution (of dwell time), it is not appropriate 
to regard it as an averaged quantity. It is widely believed that (6.19) is the 
(average) time spent by a particle with wave number k in the barrier region 
irrespective of whether it is finally reflected or transmitted.^^ 

The issue of dwell time is one of the major sources of controversy. A 
question arises when one tries to generalize (6.19) to the case of a wave 
packet. Leavens and Aers [42] pointed out that an interesting relation 

/ oo pb poo 17 

dt dx\W{x,t)\^= 7 ^ Td{k) (6.20) 

follows from the results of Hauge, Falck, and Fjeldly [43], where '^(fc) is 
the A:-space wave function obtained from ^{x^t) via Fourier transform. The 
right-hand side of (6.20) is a clear expression as an average. Indeed, it is the 
k average of ra(fc) with the A:-space probability distribution \ijj{k)f . Leavens 
and Aers [42] seem to have considered that the meaning of the space-time 
integral [the left hand side of (6.20)] was undoubtedly an exact expression 
for the average time spent in the barrier region a < x < 5 by a particle 
with wave function lF(x, t). They then read (6.20) from left to right to justify 
the interpretation of Td(fc) as the dwell time for a given k > 0 (|'^(fc)p ^ 0 
was assumed for k < 0). However, as pointed out by Hauge and Stpvneng 
(Appendix B of [2]), who although are in favor of the above interpretation 
of the space-time integral, the meaning of the space-time integral is not 
self-evident in conventional quantum mechanics. 

Despite this subtlety, the space-time integral is now widely believed to 
be the correct expression for the average dwell time when the particle is de- 
scribed by a wave packet. This is probably because the space-time integral 
has been reached in several different ways [the range of the time integral is 
not necessarily (— oc,oo)] in the analysis of tunneling time or related mat- 
ters. Baskin and Sokolovski [44] and Sokolovski and Baskin [45] obtained the 
space-time integral from a Feynman paths average of a traversal time op- 
erator that is non-local in time and acts on Feynman trajectories. Leavens 
and Aers [46] showed that, within Bohmian mechanics, the space-time in- 
tegral has a clear meaning as an average dwell time. Jaworski and Wardlaw 
[47,48] also postulated the space-time integral as giving the mean sojourn 

The word ‘average’ (or ‘mean’) is widely abused in the literature. We often en- 
counter phrases such as ‘average’ dwell time, ‘average’ transmission time, and 
‘average’ reflection time without finding the relevant probability distributions. 
The word ‘average’ should be used only when the relevant probability distribu- 
tion exists. Otherwise, one should write, for example, ‘a measure of dwell time’ 
instead of ‘the average dwell time’. 

It seems to the present author that by taking this for granted one goes beyond 
the original meaning of Td{k) so clearly stated by Landauer. 




156 N. Yamada 



(dwell) time in one-dimensional scattering. At the same time, they well rec- 
ognized that the origin and interpretation of the space-time integral is not 
self-evident because of the absence of real physical trajectories in quantum 
mechanics [48]. 

Is the space-time integral indeed a general expression for the average dwell 
time for the particle described by a wave packet? Or is it just a beautiful way 
of expressing the right-hand side of (6.20)? The right-hand side is indeed the 
k average of ra(fc) but it is not a quantum mechanical average since it has 
not been derived from a general prescription for calculating an average from 
amplitudes and operators. The present author considers that, in order to in- 
terpret the space-time integral as an average dwell time, one must show the 
underlying probability distribution of dwell time within conventional quan- 
tum mechanics. From this viewpoint, it is safe to say that the status of the 
space-time integral is still not clear. Apart from these conceptual questions, 
the dwell time for stationary states is widely used in the literature, especially 
in the analysis of transport properties of semiconductor heterostructures. See, 
for example, [49]. 

Another question that has attracted people is whether or not one can 
decompose the dwell time into transmission and reflection components as 

Td = PtTt + PrTr , ( 6 . 21 ) 

in which the dwell time is supposed to be given by (6.19) or (6.20), rt is 
the transmission time, the time a transmitted particle spends in the barrier 
region, and is the reflection time, the time a reflected particle spends in the 
barrier. Pt is the transmission probability and Pr the reflection probability 
(Pt Pr = !)• When discussing the above decomposition, we give a special 
status to the dwell time. That is, we simply accept (or assume) that the 
left-hand side of (6.21) is given by (6.19) or (6.20) with the interpretation 
that it is the average dwell time for a particle in a stationary state or in a 
localized state described by a wave packet. This was criticized by Olkhovsky 
and Recami [3]. Landauer and Martin [1] and Biittiker [4] also criticized 
(6.21) as a meaningless relation in quantum mechanics. On the other hand, 
Hauge and Stpvneng [2] took the ‘probabilistic identity’ (6.21) for granted 
and used it to single out the correct transmission and reflection times from 
various times that had been proposed. 

In Sect. 6.2.2, we considered the transmission time Tphase by following the 
peak motion of a transmitted packet. We can also follow the peak of the 
reflected packet to obtain the corresponding reflection time. For the time 
being let us denote these transmission and reflection times by rf and 
respectively. Hauge, Falck, and Fjeldly [43] found that the transmission phase 
time r/ and the reflection phase time rf satisfy the following equation that 
has an interference term on the right-hand side: 

Td{k) = + ^y/PrSm{0 - 2ka ) , 



( 6 . 22 ) 
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where Pt = Pv = \P\^-> and (3 is defined by P = |P|e^^ [see (6.3)]. This 
relation was generalized in [2] to the case of a wider region of space including 
the barrier region, to stress the meaningfulness of asymptotic phase times 
when they are properly averaged. 

Below, we review other approaches from the viewpoint of whether they 
yield transmission and reflection times that satisfy (6.21). 



6.2.7 Feynman Paths Approach 



If a particle follows a well-defined path in space-time, the tunneling time 
problem disappears. But there is no such real physical path in quantum me- 
chanics. Instead, quantum mechanics can be formulated with various types 
of quantum trajectories: Feynman paths, Bohm trajectories. Nelson trajecto- 
ries, and Wigner trajectories. The first three are trajectories in space-time, 
while the last one is defined in phase space. It is natural to attempt to use 
these trajectories to define and analyze the tunneling time. Here we examine 
a Feynman path approach. 

Consider Feynman paths connecting given end points (xi, ti) and (x 2 , ^ 2 ) 
[ti < ^ 2 ). We start by considering the sum over those Feynman paths that 
spend a certain amount of time r in the barrier region: 






(6.23) 



where Xlpaths stands for a sum over paths, the denominator (introduced to 
normalize ^ 2 , 1 ) is the sum over all paths and is therefore the propagator 
between the end points, Tah[x{')] is the amount of time a Feynman path x{t) 
spends in the barrier region [for its precise definition, see (6.53)], and the 
initial point x\ is assumed to be on the left side of the barrier, i.e., X\ < a. 

Now, F 2 ,i(r) is a distribution of r at the amplitude level (a complex- 
valued ‘amplitude distribution’). As pointed out by Landauer and Martin [1], 
it is not clear what procedure is needed to extract physical information about 
the time under the barrier from such an amplitude distribution, which is not 
the usual wave function. Nevertheless, since JdrF(r) = 1, let us formally 
regard F 2 ,i(r) as a ‘probability’ distribution for r and calculate its ‘aver- 
age’ according to J drrF(r). We expect to obtain a transmission time for 
X 2 > b and a reflection time for X 2 < a, which we shall denote by ti and T 2 , 
respectively. Sokolovski and Baskin [45] (see also [44]) found 



r\ 




dx 



81nT 

SV{x) 



— ih 



dluT(Vo + V,k) 
dV 



v=o 



(6.24) 



The approach presented here is mainly due to Sokolovski and coworkers and also 
to Fertig. For Feynman paths approaches by other authors, see [50-52]. For the 
Nelson trajectory approach, see [53-55]. For the Wigner trajectory approach, see 
[56,57]. The Bohm trajectory approach by Leavens is discussed in Sect. 6.2.10. A 
Bohm trajectory approach by other authors is found in [58]. 
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and a similar expression for T 2 {k) (just replace T above by R), where Vq is 
the height of the potential under consider at ion A remarkable feature of 

these times is that they satisfy [44,45,88] 

^d(^) = Ptnik) + PrT2(k) . (6.25) 

However, this is not the desired decomposition because Ti(k) and T2{k) are 
complex- valued. Note that Sokolovski and Baskin’s complex tunneling times 
are written in terms of the potential height derivative of the transmission am- 
plitude. Earlier, Poliak and Miller [59] extracted another complex interaction 
time from the flux-flux correlation function, a tool widely used in chemical 
physics to calculate rate constants and reaction rate probabilities in chemical 
reactions. The resultant complex time is similar to Sokolovski and Baskin’s 
but the derivative is with respect to the energy of the particle rather than 
the potential height. Martin and Landauer [60] also derived complex times 
written in terms of energy derivatives of the transmission amplitude. 

Baskin and Sokolovski also found a close relationship between their com- 
plex times and the Larmor times [44,45]: 

Ti{k) = Tty{k) - \nz{k ) , T2{k) = Try(fc) - \T^z{k) , (6.26) 

where rty(fc) and Tiz{k) are, respectively, Ty and in (6.13) and (6.12), and 
Trz{k) and Try{k) are the Larmor times for reflection [28]. Since the left- 
hand side of (6.25) is positive deflnite, taking the real part of (6.24) gives 
an equation of the form (6.21) with real- valued transmission and reflection 
times. We can also take the imaginary part. Equation (6.25) can then be split 
into^^ 

Td{k) = PtTty{k) + PrTj:y{k) , (6.27) 

0 = PtTtz{k) + PvTrzik) . (6.28) 

Now, both Tty and Try are real. Unfortunately, however, r^y is not positive 
definite and therefore (6.27) is not the desired decomposition.^® There are 
still other reasons to doubt the status of Tty and T^y as physically reliable 
transmission and reflection times. 

Note that Ti{k) and T2{k) do not depend on the end points. Sokolovski and Baskin 
[45] derived an approximate expression for the propagator in the presence of a 
barrier. Using this expression, it is easy to see that the denominator and the 
numerator on the right-hand side of (6.23) depend on the end points in the same 
way. Thus the final results are independent of the end points. 

For a rectangular barrier. Tty = Try because T and R have the same phase factor 
(apart from a constant). Using this property (and also Pt-\- Pr = 1) in (6.27), we 
obtain Td{k) = Tty = Try for a rectangular barrier. The second equation (6.28) 
expresses conservation of angular momentum. These were pointed out in [28]. 

It should be noted, however, that the positivity requirement on tunneling times 
is somewhat subtle. See, for example. Sect. 1 of [5] and also p. 806 of [61]. 
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• Leavens and Aers (Sect. IV of [62]) and Range and Stdvneng (Sect. VI C3 
of [2]) presented an example where taking Tty as a physically reliable 
transmission time betrays our intuition that ‘a reflected particle spends 
exactly zero time on the far side of a barrier’ [62]. 

• For an opaque barrier, Tty is almost independent of the barrier width. 
Thus, for a sufficiently thick barrier, the speed Vt under the barrier cal- 
culated from the classical equation Vt = d/Tty can be greater than the 
speed of light. Leavens and Aers [62,63] showed that this is also the case 
if the Dirac equation is used instead of the Schrodinger equation, thus ex- 
pressing a serious doubt about the status of Tty as a physically meaningful 
transmission time [5,62]. 

These arguments remind us of the status of (6.21), i.e., that it is only a 
necessary condition for tunneling times even if it is correct. Leavens and 
Aers [5] proposed the following set of conditions to be satisfled by physically 
meaningful tunneling times (for a stationary case): 

(i) (6.21) or (ii) its generalization 



Td(A:;xi,X2) = PtTt{k]Xi,X2) + Pr^r(fc; xi, X 2 ) , (6.29) 

where Pt = |T(fc)p, P^ — \R{k)\^, and Ti{k]Xi,X 2 ) is the dwell time 
(i = d), transmission time (z = t), and reflection time (i = r) for the 
interval x\ < x < X2, 

(iii) additivity 

Ti{xi,x^) =Ti{xi,X 2 ) ^Ti{x 2 ,x^) for Xi < X2 < X3 , (6.30) 

(iv) positivity 

Ti{xi,x^) > 0, (6.31) 

(v) consistency with causality for rt. 

Sokolovski and coworkers made an extensive contribution to the tunneling 
time problem by using Feynman paths. Readers are also referred to [64-68] 
for their work, and also to [69], in which the Feynman paths approach and the 
Rohm trajectory approach to be discussed later are reviewed and compared. 

6.2.8 Systematic Projector Approach 

Brouard, Sala, and Muga [70] (see also [71,72]) have developed a systematic 
method for decomposing the ‘mean’ dwell time [postulated to be given by 
the space-time integral in (6.20)] into transmission and reflection times. This 
systematic projector approach yields an inflnite number of resolutions of dwell 
time, giving, in general, equations of the form 



Td = PtTt P PrR + interference term . 



(6.32) 
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The systematic projector approach contains three projection operators P, 
1 - P), where 1 is the identity operator, and D. The projector P se- 
lects that part of associated with transmission (in the sense that 

the selected part has only the positive momentum components in the infi- 
nite future), and the projector Q selects the part associated with reflection 
(only negative momentum components in the infinite future). The expecta- 
tion values of P and Q are, respectively, the usual transmission and reflection 
probabilities defined in the limit t oc and hence independent of time. The 
projector D selects the part of ^{x,t) inside the barrier region a < x < b 
so that its expectation value (D) = {^{t)\D\^{t)) is the usual probability 
of finding the particle in the barrier region at time t. The left-hand side of 
(6.20) can thus be written as J {D) dt. The essence of the systematic projector 
approach is to rewrite D in f (D) dt using the property P + Q = 1 and the 
idempotency of a projection operator (e.g., D = DD). For instance, we can 
write D = PP + QjD, D = DP + DQ, and D = PDP-\-QDQ-\-PDQ-\-QDP. 
Using the first one in f {D) dt leads to [71] 

Td = + PvT^^ , (6.33) 

where 

rr = ^J m)\PD\nt)) dt, i Jmt)\QD\m) dt . ( 6 . 34 ) 

Equation (6.33) does not have an interference term, but and are 
complex-valued since PD and QD are not Hermitian. Muga, Brouard and 
Sala [71] found that the two complex times are nothing but Sokolovski and 
Baskin’s complex- valued transmission and reflection times, i.e., = ti 

and = T 2 [see (6.26) for ti and T 2 ]. 

Now, the second choice D = DP-\-DQ leads us to another decomposition 
similar to (6.33) with new transmission and reflection times, which we denote 
by and respectively. Since D and P do not commute, ^ 
and and therefore different decompositions of D lead to different 

expressions for rt and in general. The third choice [71] 

D - PDP -f QDQ + PDQ + QDP 

leads us to yet another decomposition of the form of (6.32) with the interfer- 
ence term given by 2Re[rint], where 

Tint = I {<I'{t)\PDQ\^{t)) dt , (6.35) 

which might be called an interference time. The real- valued times 
tQDQ ^ and Tint obtained from the present decomposition coincide with the 
transmission, reflection, and interference times introduced by van Tiggelen, 
Tip, and Lagendijk [73]. In this way, the systematic projector approach repro- 
duces many of the proposed transmission and reflection times in a systematic 
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way. Leavens [61,74] has shown, however, that the systematic projector ap- 
proach does not give positive definite transmission and refiection times that 
satisfy both the probabilistic identity (6.21) and the condition for additivity 
(6.30). The utility of the systematic approach is that it provides, completely 
within standard quantum mechanics, ‘a compact classification scheme for 
many of the proposed times and simplifies their connection and study’ [23]. 

In the systematic projector approach, quantities arise naturally that we 
are tempted to interpret as joint probabilities. In the third choice discussed 
above, we encounter a quantity {^{t)\PDP\W{t)). This has a suggestive form 
as the joint probability that The particle is now transmitted and was in 
the barrier region before. Similarly, {^{t)\DPD\^{t)) is a suggestive form 
as the joint probability that The particle is now in the barrier region and 
will eventually be transmitted’. The two joint probabilities must coincide in 
classical probability theory, but they do not in fact because D and P do not 
commute. This simply means that the concept of joint probability cannot 
be introduced in quantum mechanics when the relevant observables do not 
commute. In a word, the ‘sample space’ for the desired joint probabilities is 
not an exhaustive set of mutually exclusive events. Rather, the ‘sample space’ 
consists of interfering alternatives, for which probabilities are not definable 
in a consistent way as in the case of the two-slit experiment. Even so, the use 
of such ‘generalized probabilities’ which are not actual probabilities is still 
interesting if one can extract useful information and/or if it throws new light 
on known results. 



6.2.9 Weak Measurement Approach 

Steinberg [75,76] and lannaccone [77] have shown that the complex- valued 
conditional probabilities motivated by the weak measurement theory of 
Aharonov, Albert, and Vaidman [78,79] yield some of the known results about 
tunneling time. Note that (6.21) is immediate from 

= Pi P{x,t\trsins) PrP{x,t\ie&) , (6.36) 

provided that the conditional probabilities P(a;, tjtrans) and P(x,t|refl) are 
definable, where the former is supposed to be the probability distribution 
of the particle’s position at time t conditioned on the information that the 
particle is found to have been transmitted at sufficiently later times, and 
similarly for P(x, tjrefi). An equation of the form of (6.36) is immediate from 
the identity 

{i\A\i) = J2\{n\i)f where A ^ ^ . (6.37) 

Let us put, with h = 1, 



(6.38) 
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/ Ax\x){x\, (6.39) 

J a 

(n| = (out„(i')|e“‘lt (6.4Q) 

where |^(0)) is the initial state of the particle (far to the left of the barrier) 
and {(outn(t')l \n = t,r} = {(trans(t')|, (refl(t')|} is a set of transmitted and 
reflected states at a later time (0 < t < t'). Substituting (6.38)-(6.40) into 
(6.37) and taking the limit t' oo, we obtain (6.36) integrated over x from 
a to b. Identifying the resultant equation multiplied hy m/hk with (6.21), we 
have^^ 



* hk^'^ 



m 



m {tvdins{t)\A\^{t)) 
hk (trans(t)|l^(t)) 
b 



1 



hk (trans(t)|lZ^(t)) 



J a 



(6.41) 



where (trans(t)| = (trans(t')|e~' = (x|trans(t)), andlZ^(x,t) = 
(x\^{t))}^ We also obtain a similar expression for Let us consider a func- 
tion a{x) that is unity if x is inside the barrier and zero otherwise. Equa- 
tion (6.41) shows that rt can be viewed as the conditional average of a{x) 
calculated with the conditional probability 



P(x, tjtrans) 



^^{x,t)^{x,t) 

(trans(t)llZ^(t)) 



(6.42) 



Note that this is a spatial average, not a temporal average. 

We saw that it is actually possible to construct an equation of the form 
(6.21) from the probabilistic statement (6.36) with explicit expressions for the 
conditional probabilities P(x,t|trans) and P(x, tjtrans). However, the condi- 
tional probabilities and the resultant times (rt and r^) are in general complex- 
valued, so that we cannot take those quantities literally. From the viewpoint 
of dwell time decomposition, the status of the present method is similar to 
that of the systematic projector approach. 

Steinberg [75] has shown that rt given by (6.41) agrees with Sokolovski 
and Baskin’s complex tunneling time ri given by (6.24). See also [77] for 
the connection between the times obtained by the present method and those 
obtained by other methods. In the weak measurement theory, is called 
the weak value of A in the time interval between pre-selection of the state 
|?Z^(0)) and post-selection of the state |outn(t^)). In the weak measurement 
theory, the system of interest interacts with an apparatus (measuring device) 
in a certain manner and what we measure is the position of the pointer of the 

The superscript t attached to Aw, which simply means Transmitted’, must not 
be confused with the lower bound t of the s integration. 

(trans(t)|lZ^(t)) is in fact time independent. See [75] for details. 
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apparatus. A weak value then has the following meaning: when the pointer 
position is measured, the real part of the weak value gives the mean shift 
in the pointer position, while the imaginary part gives a shift in the pointer 
momentum. It thus follows that in the weak measurement theory the real 
part and the imaginary part of the complex tunneling times are related to 
such physical shifts in apparatus variables. 



6.2.10 Bohm Trajectory Approach 



Leavens [46] has shown that it is indeed possible within Bohmian mechanics to 
find positive definite Tt and Tj- that satisfy (6.21). It is a direct consequence 
of the basic property of Bohm trajectories that they never intersect after 
starting from different initial positions. In Bohmian mechanics, a quantum 
mechanical particle has a well-defined position x{t) and velocity v{t) at each 
time. The velocity field is given by 



v{x, t) 



J{x, t) 



x=x{t) 



(6.43) 



where J is the usual probability current density. Solving the particle’s equa- 
tion of motion dx{t)/dt — v{x{t)^t) yields a causally determined trajectory 
x{x^^\t)^ where x(0) = is the particle’s initial position. The dwell 

time for a Bohm trajectory x{x^^\t), the time the trajectory spends in the 
spatial region of interest, is defined by 

poo pX 2 

rd(xi,X 2 ;x^°^) = / dt dx 6[x — x{x^^\t)) . (6.44) 

Jo Jxi 



Since the initial positions of the trajectories distribute according to the prob- 
ability distribution 0)p, the average dwell time is given by 

/ CX) 

. (6.45) 

-CX) 



Substituting (6.44) into (6.45) and carrying out the integral first, we 
obtain 

pOC pX2 

Td(xi,X2)= / dt dx\^{x,t)\‘^ . (6.46) 

Jo Jxi 

The left-hand side of (6.20) is indeed, within Bohmian mechanics, an average 
dwell time. To decompose rd(xi,X 2 ) into transmitted and reflected compo- 
nents, we begin by noting that there exists a special trajectory Xc{t) that 
separates transmitted and reflected trajectories at any time, as implied by 

Let us consider here the dwell time in a general spatial region x\ < x < X2- The 
dwell time in the barrier region is immediate by putting x\ — a and X2 = b. 
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the non-intersecting property of Bohmian trajectories. The trajectory Xc{t) 
is defined by 

pOO 

/ dx \'F{x,t)f — transmission probability (Pt) ■ (6-47) 

Jxc{t) 

We can decompose the probability density as 

\^x, t)f = 0{x - x^mnx, t)\^ + 0{xS) - x)mx, t)f . (6.48) 

Since Xc{t) separates transmitted and reflected trajectories at any time t, the 
first term on the right-hand side of (6.48) is the to-be-transmitted compo- 
nent of |^'(a;,t)p, while the second term is the to-be-reflected component. 
Integrating both sides of (6.48) with respect to x and t, we immediately 
obtain 



Td{xi,X2) = PtTt{Xl,X2) + PvTr{Xi,X2) , (6.49) 

-j pOO pX 2 

Tt(a;i,a; 2 ) = — / dt dxO{x - Xc{t))\^{x,t)\^ , (6.50) 

Jo Jxi 

1 poo nX2 

Tr(xi,X 2 ) = — / dt dxO{xc{t)-x)\^{x,t)\^ . (6.51) 

Tr Jo Jxi 

Equation (6.49) has the same form as (6.21) with positive definite times on its 
right-hand side. Moreover, ra, rt, and above satisfy the additivity condition 
(6.30). For now, the use of Bohmian mechanics is the only way to arrive at 
this remarkable result. Leavens, at the same time, carefully points out the 
price one has to pay. That is, the connection between the case of a wave 
packet and the case of a stationary state is not simple: 

n{xuX2) ^\mf\mW{x,,X2-,k), (6.52) 

where rt(xi, X 2 ; k) is the transmission time for a stationary case, and similarly 
reflection time. The lack of the simple relation is due to the fact that rt(xi, X 2 ) 
and Tr(xi,X 2 ) derived above are not bilinear in because they depend 

on Xc{t), and Xc{t) depends on ^{x,t) through (6.47). 

The separability of Bohm trajectories into to-be-transmitted and to-be- 
reflected parts is at the heart of the derivation of (6.49)-(6.51). This prop- 
erty is often criticised as unphysical. However, the possibility remains that 
the results obtained are correct and justiflable even if the separation itself 
is unphysical. At the same time, the status of (6.49)-(6.51) in conventional 
quantum mechanics needs to be clarified. More fundamentally, we must keep 
in mind, before arguing the physical meaning of (6.49)-(6.51), that the de- 
composition (6.21) itself could be meaningless in quantum mechanics. We 
should not take (6.21) for granted, simply because our classical intuition be- 
hind it could be wrong for purely quantum mechanical phenomena. 
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For Leavens’ earlier contributions regarding Bohm trajectories, see [80,81]. 
One of the conclusions is that a Bohm- trajectory calculation of tunneling time 
gives answers that are very different from those obtained from the Larmor 
clock or the time-modulated barrier method by Biittiker and Landauer. 

6.3 Tunneling Times as Interfering Alternatives: 

A Feynman Paths Approach 

6.3.1 Lesson from Two-Slit Experiments 

I claimed in [9] that we should understand the nature of tunneling time from 
the viewpoint of interfering alternatives and argued that what makes sense 
in the tunneling time problem is the range of tunneling times. An important 
point is that probabilities cannot be considered in the range. We cannot say 
that the tunneling time takes a certain value in the range with such and such 
a probability. This is because different values of tunneling time are interfering 
alternatives, i.e., they are not exclusive. Recall that probabilities are defined 
only for an exhaustive set of mutually exclusive alternatives. An analogy 
with the two-slit experiment would be helpful to understand this point. In 
the two-slit experiment, due to the interference between the two alternatives 
of a particle going through the upper slit and going through the lower slit, 
probabilities cannot be defined for the two alternatives without disturbing 
the motion of the particle (see Fig. 6.7). 

In the tunneling time problem, the y axis (at x = 0) on which we con- 
sidered position alternatives corresponds to the time axis on which we con- 
sider tunneling time alternatives (alternatives correspond to different values 
of tunneling time). To discuss quantum interference between tunneling time 
alternatives, we need a measure of the interference. We can find it by work- 
ing in the real time formulation of Feynman path integrals. The starting 




Fig. 6.7. Left: we cannot say that ‘a particle passes through slit 1 or slit 2 with 
some probability’ in the two-slit experiment. Right: we can say, however, that ‘a 
particle passes through the slit of width ISy with probability approximately unity’ 
in the one-slit experiment. Therefore, although the precise point the particle passes 
through cannot be considered, the range of space (At/) it passes through can be 
considered 
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Fig. 6.8. The time a Feynman path spends in the barrier 

point is the amount of time t[x(-)] which a Feynman path x(-) traversing 
the barrier spends inside the barrier region. Since Feynman paths can cross 
the two points x = a and x = b many times (in Fig. 6.8, the single crossing 
case is shown), several definitions are possible for r[x(*)]. In [9], I used the 
‘resident-time’ definition widely used in the literature 

'r[a^(-)] = [ , (6.53) 

Jto 

where 0ab{^) is unity for a < a; < 6 and zero otherwise. 

Since there is no classical solution to the equation of motion in the tunnel- 
ing regime, many paths contribute to the tunneling propagator and the value 
of r[x{')] is not unique. We can classify the paths according to the values of 
r[x(-)]. The sum over paths defining the tunneling propagator can then be 
written as a time integral over all values of r[a:(-)], and as a result we can 
write the transmitted packet as 

^{x^t) = j dr^(a),t|r), (6.54) 

where lZ^(x,t|r) is the product of the ‘partial propagator’ K(x, t; xq, to|r) to 
which only the paths satisfying r[x(-)] = r contribute and the initial state 
?F(xo,to) that is nonvanishing only on the left-hand side of the barrier, fol- 
lowed by an xq integral. What corresponds to (6.54) in the two-slit experiment 
is 



^ ^ (6-55) 

i=i,2 

where ^{y,t) (omitting x dependence) is the total wave function after the 
slits and ^j{y,t) is the partial wave function from slit j {j = 1,2). The 
probability that a particle is found at point y on the final screen at time t is 
given by \^{y,t)\‘^ = \^i{y^t) -h ^ 2 {y^t)\‘^j where x implicit in ^ is supposed 
to correspond to the screen position. The alternative that a particle goes 

In [82] I used a ‘passage-time’ definition and applied the CHA criterion to the 
tunneling time problem for the first time. The passage-time definition was first 
employed by Schulman and Ziolkowski [52]. 
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through slit 1 and arrives at y and another alternative that it goes through 
slit 2 and arrives at y are not exclusive because 

+ + ■ (6.56) 

The cross term 

Re'^l{y,t)'p2{y,t) (6.57) 

is the measure of interference between the two alternatives. 



6.3.2 Tunneling Time Probability Distribution is not Definable 



Now, for the tunneling time problem, we follow the same line as above. The 
tunneling probability, being the probability that a particle is found some- 
where on the other side of the barrier at sufficiently later times, is given 

by 



poo 

P= lim / dx|^(x,/:)p. 


(6.58) 


Substituting (6.54) into (6.58), we can write 




P = J dr J dr'Z?[r;r'], 


(6.59) 


where 




pOO 

D[t;t']= lim / dx (x,t\T)'I'{x,t\T') . 


(6.60) 



Note that only the real part of D contributes to the double r integrals since 
jD*[r;r'] = D[t';t]. Now, ReD[r;r'] gives the measure of interference be- 
tween the two alternatives of a particle spending time r and spending time 
r' in the barrier region. It is, of course, a subject for investigation whether 
the interference vanishes or not, i.e., whether an equation of the form 



Re D[r; r] = — t')P{t) 



(6.61) 



holds or not. Equation (6.61) reads as follows: ReD[r;r'] is proportional 
to (5(r — t') with a proportionality coefficient Pir). To clarify the physical 
meaning of P(r), let us see what happens if we actually assume (6.61). Under 
this assumption, one finds 



P(r) > 0 and 




(6.62) 



which are exactly the properties that a probability distribution of tunnel- 
ing times must have. Based on this observation, I consider (6.61) to be the 
criterion to judge the possibility of a probabilistic description of tunneling 
time: If (6.61) holds, a probability distribution of tunneling time is definable 
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and it is given by the proportionality coefficient P{r); otherwise, a probability 
distribution is not definable. 

To investigate whether (6.61) holds or not, we must calculate the right- 
hand side of (6.60) using a concrete expression for t\r). This was done in 
[9] using (3.3) in [66], and I shall only repeat the consequence here. D[t;t'] 
was found to be 

D[r; r'] = , (6.63) 

where 

%{t) = I (6.64) 

and (•••) = J dfc • • • with '0(fc) the fc-space wave function of the in- 

cident particle and T{V, k) the transmission amplitude when the potential 
height is V. The range of the V integration is (— 00 , 00 ). To guarantee that the 
incident particle moves towards the barrier with energy less than the barrier 
height, fj{k) is so restricted that it vanishes for fc < 0 and for k > y/2mVo/h, 
and satisfies the normalization J dfc|'0(fc)p = 1. I found that (6.61) with 
(6.63) leads to 

Pvo-v + Pvo-ev = 2 Pvo 5 (6.65) 

for arbitrary V and then showed that (6.65) leads to a contradiction, namely 
that the tunneling probability is higher for higher barriers. This means that 
(6.61) does not hold. The physical conclusion is, therefore, that a probability 
distribution for tunneling times is not definable since different values of r 
interfere. 

6.3.3 Relation to Gell-Mann and Hartle’s Formulation 
of Quantum Mechanics 

In terms of Gell-Mann and Hartle’s generalized quantum mechanics (GQM) 
(or, the decoherent histories approach to quantum mechanics) [13,14], D[r] r'] 
and (6.61) are, respectively, the decoherence functional and the weak deco- 
herence condition (WDC) [83]. GQM provides a general scheme for the study 
of quantum probabilities for ‘histories’, which are alternatives not restricted 
to a moment of time. The negative result reached above may be regarded 
as a result of GQM with continuous histories (Feynman paths), rather than 
discrete histories (the products of projection operators at different times). Al- 
though GQM with continuous histories has received less attention, our study 
demonstrates how it can be applied to concrete problems. 

6.3.4 Td = PtTt + Is Doubtful 

In situations where discussions on the tunneling time problem do not con- 
verge, some authors proposed criteria for ‘correct’ tunneling times. The most 
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widely used among them is the following ‘probabilistic identity’ introduced 
by Hauge and Stpvneng in [2]: 



"Td = Ptn + -PrTr , (6.66) 

where ra is the dwell time, the total average time a particle spends in the 
barrier region irrespective of whether it is finally reflected or transmitted, 
Tt is the transmission time, the average time a transmitted particle spends 
in the barrier, and is the reflection time, the average time a reflected 
particle spends in the barrier. Pt is the transmission probability and the 
reflection probability, satisfying Pt + Pr = 1. Here, as well as in [9], we 
concentrate on the transmission time, calling it the tunneling time (our P 
is Pt). In many arguments, ra, R, and in Eq. (6.66) are assumed (often 
implicitly) to be averaged quantities. In other words, the existence of the 
probability distributions of these times are often assumed behind Eq. (6.66). 
In fact, Eq. (6.66) seems most rigorously founded if such distributions exist. 
As concluded above, however, the probability distribution is not definable, at 
least, for Tt. (It is also unlikely that and ra allow probabilistic descriptions 
because WDC is a very severe condition.) Thus, the validity of Eq. (6.66) is 
doubtful. In my opinion, many arguments in which Eq. (6.66) plays essential 
roles should be reconsidered. 

At the same time, it should be pointed out that our negative result does 
not exclude the possibility of a distribution of r that is not fully inherent 
in tunneling but is valid only with a certain measurement setup. Although 
such non-intrinsic distributions with limited validity are not the subject here, 
the probabilistic identity could find its meaning with such non-intrinsic dis- 
tributions. Sokolovski and coworkers [84,85] have proposed distributions of 
tunneling times by using a Larmor clock as the measuring device, where the 
angle of spin rotation due to a magnetic field in the barrier is measured and 
compared with theory. 

6.3.5 Range of Tunneling Times 

From a classical point of view, it is hard to understand in the two-slit exper- 
iment that going through the upper slit and going through the lower slit are 
not exclusive, but this is how nature behaves and we just accept it. In the 
same way, we must accept the non-exclusiveness of tunneling times and give 
up a probabilistic description of tunneling time. This should not be regarded 
as a ‘discouraging’ result. It just shows how nature behaves. What we need to 
do is to consider what is meaningful when nature behaves that way. An idea 
arises from the one-slit experiment (see Fig. 6.7): just as the spatial range 
Ay the particle goes through can be meaningfully considered, we should be 
able to consider the range of tunneling times. Let us denote by Cr the class 
of Feynman paths for which r[x(-)j = r. Roughly speaking, the range would 
be such that the sum over paths for the tunneling process is approximated 
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r< Ty 



Fig. 6.9. How to estimate the range of tunneling times using the function G{r) 



with sufficient accuracy by the sum over those paths that belong to Cr with 
T exhausting all values in the range. More precisely, when 




dr' D[t\ t] 



s 



p 



(6.67) 



holds with sufficient accuracy for Ms < r< and for Vt > r>, I propose to 
regard (r<,r>) as the range of tunneling times. Rather than dealing with 
dr dr' D[t; r'], which is a function of two variables s and t, let us work 
with a function of one variable 



G(t) 




dr"a[r';r"], 



0 



(6.68) 



which satisfies G(0) = 0 and G(oc) = 1. The range of tunneling times can 
then be estimated by finding r< and r> such that, for e <C 1, 



G(t) < e for Vr < r<; |1 — G(r)| < e for Mr > Ty . (6.69) 



See Fig. 6.9). 

It is interesting to see that, if WDC [(6.61)] held (though it does not in 
fact hold), the range of tunneling times would simply be the support of the 
probability distribution of tunneling time. This is because, assuming (6.61), 
(6.67) reduces to 



/ drP(r)wP. 



(6.70) 



It would be natural to call the support (s, t) of P{r) the range of tunneling 
times if P{r) were definable. The concept of the range of tunneling times can 
be introduced, however, regardless of the existence of P(r). In (6.67), the 
range of tunneling times is defined at the level of amplitudes, not at the level 
of probabilities. 

The function G(r) can only be obtained analytically for delta function 
potentials, so we must rely on numerical calculations. From (6.63) and (6.68), 




J — oo 
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P = jdk\^{k)f\T{Vo,k)f. (6.72) 

It is clear that the initial state of the incident particle [or, ^{k)] and the 
barrier shape [or, T(R, /c)] determine the function G(r) and thus determine 
the range of tunneling times. More fundamentally, the decoherence functional 
D[r-y] is written in terms of the initial state and the barrier shape [see 
(6.63)]. It is quite reasonable that the two elements (the particle and the 
barrier) finally determine speakables and unspeakables in the tunneling time 
problem. Equation (6.71) is not only of conceptual importance but also of 
practical importance as a basis for numerical calculations of G(r). 



6.3.6 An Example: The Monochromatic Case 



Here, following [9], let us consider the case of an incident particle having a 
constant energy. We thus put |'0(fc)p = 5(/c — fco), which leads to 



G(r) = 



T{Vo,k 



Oj Jinfty 



dV sm{VT/h) 

E 



T(Eo + E,fco) 



(6.73) 



where T(E, k) is given by (6.4) with repacement Vq V. Inside the square is 
the function C(r) introduced by Fertig [86,87]. To be precise, our formulation 
loses its validity when the incident energy is constant because the assumption 
of well-defined refiected and transmitted packets, which is used in the deriva- 
tion of (6.63), is not satisfied in such a case. Therefore, |'0(fc)p = 6{k — k^) 
should be understood as \'ip{k)\‘^ S{k — ko). Figure 6.10 shows a typical 
example of G(r) for opaque barriers in the monochromatic limit, where the 
free passage time of the region, rf = md/hko^ is used as the time unit. The 
parameter values are kod =1 {d = b — a) and Vq/Eq = 5 [Eq = {hko)‘^ /2m], 
for which P = \T{Vo,ko)\‘^ = 4.64 x 10"^. From Fig. 6.10, we can read, by 
choosing e = 0.01 for instance, r</rf 0.0335 and T>/rf 5.95. 

Among various tunneling times proposed in the literature, the most fa- 
mous ones are the Larmor time tlm — ^/Eoy/Vb/Fo — 1 [29,30] and the 
Buttiker-Landauer time tbl = d^/m/2{Vo — Eq) [34]. (Another famous one 
is the phase time.) In Fig. 6.10, LM and BL indicate the Larmor time and the 
Buttiker-Landauer time, respectively. Note that r< ^ tlm and r> ^ tbl- 
This is true for a wide range of parameters corresponding to opaque barriers. 
If both T< and r> were close to tbl (or Ihon we would be able to regard 
r’BL (or Tlm) as The tunneling time’ in a loose sense, but this is not the case. 
In [52,65], Tbl was considered as a measure of the range of tunneling times, 
but the range obtained here is much wider than tbl- 

In addition to r< and r> , we can find (though only loosely) another char- 
acteristic time Tc. This is the time at which G(r) takes its maximum when its 
noisy behavior is ignored. Figure 6.10 shows Tc ~ tbl, which is also widely 
observed in our numerical results for opaque barriers. Although G(r) takes 
significantly large values around Tc, we cannot regard Cr^ as the dominantly 
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Fig. 6.10. G{t) for the monochromatic case 



contributing class of paths, since, as implied by the decreasing behavior of 
G{r) for larger r, the contribution from the class is mostly canceled by con- 
tributions from other classes of paths. The idea of identifying ‘the tunneling 
time’ by finding a sharply defined ‘dominantly contributing class of paths’ 
does not work. There is no such class. 

We have found three characteristic times associated with tunneling: r<, 
r>, and Tc- The physical meaning of these times would be clarified in situa- 
tions where the particle interacts with other degrees of freedom in the barrier. 
I would expect that and ^ determine the parts of the environment (in 
terms of its frequencies) that can affect the motion of the tunneling particle. 
The large values of G{r) around r = might be an implication that the 
tunneling probability would be greatly affected if the Feynman paths con- 
tributions from around Cr^ are altered. If so, a time-dependent interaction 
whose time scale Tint is about would cause a significant change in the tun- 
neling probability. The significant change would occur not at a specific value 
of Tint but over a range of Tint around Tint = because G{t) takes large 
values around r = 

The width of the range of tunneling times is essentially determined by 
r>. In [9], I derived an asymptotic expression for G{r) to give a formula for 
T> for opaque barriers: 






Im(Zo) 



Wu — 1 + In 



— j , u=^, l = kod, (6.74) 

2e{zo -u) \ Eo 



where zq = ZojEo, r = R/Eq = 4(1 — zq)/{2 — ilzo). Among the poles of 
T{V^k) in the complex V plane, Zq is the pole that is closest to the real 
V axis, and R is the residue of T(V, k) at the pole. A graphical method for 
obtaining zq is given in [9]. 
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6.4 Conclusion 

In my opinion, the main reason for the long-lasting controversy over the 
tunneling time problem is that it has not been noticed^^ that different values 
of tunneling time are not exclusive. Just as probabilities cannot be considered 
for the non-exclusive alternative of which slit a particle passes through in the 
two-slit experiment, probabilities cannot be considered for different values of 
tunneling time since they are not exclusive. What can be meaningfully defined 
is the range of space the particle passes through in the one-slit experiment 
and the range of tunneling times in the tunneling time problem. The range 
of tunneling times is so broad for opaque barriers that a unique tunneling 
time cannot be considered even in a loose sense. 

A significant of the range of tunneling times is suggested even from the 
preliminary investigations of Sect. 6.3; various tunneling times that were pro- 
posed in the past are found to fall within the range, implying that the in- 
formation of important time scales of tunneling is ‘embedded’ in the range. 
Moreover, I expect that the further studies of the range in the presence of 
other degrees of freedom will provide a new insight into how the coupling 
between the particle and the environments affects the dynamics of tunnel- 
ing and hence alters the tunneling probabilities or the tunneling currents. 
Through these studies, we would be able to learn how the tunneling time 
range can be measured experimentally and how the system information can 
be extracted from the range. The range of tunneling times under the influence 
of environments thus seems to be an exciting direction to explore toward the 
‘tunneling time spectroscopy’. 
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Exceptions are [84] (see Sect. X there) and [87], in which completely different 
approaches from our own were pursued for probability distributions of tunnel- 
ing time, although they also worked with real-time Feynman path integrals. A 
fundamental difference between their approaches and the present one is that we 
have a solid criterion (WDC) for judging the definability of a probability dis- 
tribution. Definability thus becomes a subject for investigation, rather than an 
assumption. Another difference is that the present approach introduces the con- 
cept of the range of tunneling times in a precise manner based on the decoherence 
functional. 
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